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ABSTRACT 


Synchronous, permanent magnet, and induction machines 
are modeled using computer programs. The computer programs 
incorporate an optimization algorithm which converges on 
lowest weight, volume, and inefficiency. Machine designs 
for high and low rpms are performed, with a varying number 
of pole-pairs. The machine designs are analyzed to find the 
optimum combination of generator and motor for inclusion in 
a naval ship propulsion system. 

Three ships are used for the system study: a baseline 
mechanical transmission ship, a ship retaining the same sub- 
division as the baseline but with the electric machinery, 
and an electric transmission ship with subdivision and 
machinery box arrangement chosen to benefit from the in- 
herent arrangeability of electric transmissions. 

Two generator/motor combinations are used in the final 
ship analysis. Both employ a 3600 rpm, six-pole synchronous 
generator, which turns at the shaft speed of the prime 
mover. One combination uses a 180 rpm, direct-drive, 16- 
pole synchronous motor, and the other uses an 1800 rpm, 
geared, 8-pole synchronous motor. Power converters are used 
in both combinations to control motor speed. 

The geared combination in the rearranged ship 
demonstrated the best endurance speed efficiency, reducing 
the endurance fuel load by 18%, while maintaining the maxi- 
mum and sustained speed of the baseline ship. The savings 
in ship volume translated to an additional twenty Tomahawk 
missile cells in the rearranged ship. When the fuel load 
was held at the tonnage of the baseline ship, endurance 
range increased as much as 25%. 

Permanent magnet machines were not competitive in this 
study due to their high weight and volume, even though their 
individual machine efficiency was the highest of all types. 
Induction machines were not used as propulsion generators 
because of the inherent difficulties in control. The induc- 
tion machine motor candidates were not competitive because 
of off-design-point inefficiency. 
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Chapter One: Introduction 


The use of electric drive as a propulsion method for 
naval ships brings to the ship design process improved ar- 
rangeability and efficiency, though electric machines may 
increase the weight of the plant. Water-cooled electric 
machines are being studied today for ship transmissions! ; 
these are smaller and lighter, for the same power output, 
than air-cooled machines. They promise reduced overall 
weight for the ship through more economic prime mover load- 
ing, as less fuel will be needed on board. A review of the 
literature has found no work comparing various types of con- 
ventional motors and the effect of each type on the overall 
ship system when used as a propulsion method. St. John [1] 
showed the effects of a superconducting generator/motor 
transmission on the design of a DD963 destroyer hull. Many 
Simulations of electric motors and their transients have 
been done. There have been several papers written on naval 
ship integrated electric propulsion systems2. Also, much 
effort has been expended in the area of electric motor 
design and optimization. Herein, various kinds of conven- 
tional electric machines are modeled. Those machines were 
used in ship designs to find the sensitivity of the designs 


to their use. 


1. Greene, Mole, Welch, and Seng, “Analysis of a High-Power 
Water-Cooled Electric Propulsion System," SNAME Trans. , 
Vol 86, 1978, pp 140-162. 


2. Ames, "Marine AC Generation Systems," LSE Journal, 
Yol 12(1), pp 13-29. 





1.1. Review of electric drive 


Navy ships operate in almost every salt water location 
in the world, including the Black Sea and the Indian Ocean. 
Regular deployments are made to the Mediterranean Sea, North 
and South Atlantic Ocean, and all areas of the Pacific. 
Ships transit the Suez and Panama Canals, the Saint Lawrence 
Seaway, and operate in the Arctic and Antarctic Oceans. 
These operations are made under greatly varying environmen- 
tal conditions, ranging from the sub-freezing temperatures 
of the high latitudes to the hot, dusty conditions of the 
Middle East. All this variety requires ships (and men) 
capable of sustained and efficient operation under any known 
condition. To that end, naval ships are tremendously com- 
plicated systems. 

It is impractical to equip naval ships for every con- 
tingency, but ship designers try to include as much as pos- 
sible when deciding what systems, equipment, and spares to 
put aboard a ship. Once a ship is designed, or as part of 
the design process, political acceptance of the product and 
its purchase is required. Since ships cost tax dollars, 
there are usually limitations on the size and complexity of 
the design. Still, the designers try to work within the 
given constraints and produce an acceptable and survivable 
(in both the battle and political senses) design. Usually 
this results in a ship that has very small margins in avail- 
able weight and volume.3 

Weight (displacement) is a semi-direct measure of ship 
cost. Volume is required to place desired systems aboard a 
ship. Therefore, any design change that results in less 
required weight or volume with no decrease in ship 


3. To have small margins is to be limited in the quantity of 
additional system weight and volume that can be added to the 
ship over its lifetime. If a ship is limited in this 
fashion, the flexibility one has in backfitting new systems 
as the ship ages is significantly decreased. 


10 





effectiveness, is usually a welcome one. 

There are several different kinds of ship propulsion 
systems now in use. They include steam, nuclear, diesel, 
and gas turbine, and there are two principle drive systems: 
mechanical and electric drive. Below is a crude comparison 
of the various propulsion systems, for the purpose of plac- 
ing the thesis in perspective. 

Steam plants burning coal or oil have been in use for 
over a hundred years. They require large amounts of prime 
ship volume, in the center part of the ship.4 Steam is 
produced in boilers and used to power turbines that rotate 
the shafts and propellers through reduction gears. Steam 
plants are very reliable mechanically, but are not terribly 
economic. The large size of the system components demands 
that the boilers and turbine machinery be placed in the cen- 
ter of the ship. This necessitates long runs of shafting 
from the center to the stern of a ship. Shafts typically 
are 18°-24" hollow steel cylinders of two to four inch 
thickness; they are heavy and their required placement and 
length makes valuable volume unavailable for other use. All 
propulsion plants except for electric drive have this ar- 
rangement and shafting disadvantage. Steam plants are used 
in all sizes of ships, from the 3000-ton displacement 
frigates to the 50,000-ton battleships to the 80,000-ton 
aircraft carriers. See Table 1 for a list of ship types and 
principle dimensions. 

Nuclear plants are steam plants with a different heat 
source. They also produce steam to power turbines and suf- 
fer the same volume disadvantages as conventional steam 
plants. They are also very heavy and very costly. Manning 
requirements are more stringent, since personnel levels are 


4. The center part of the ship is the most useful ship 
region to place and arrange systems. Ship designers try to 
keep free as much center ship volume as possible. This al- 
lows much greater flexibility in arranging systems that have 
large objects, such as boilers and turbines. 
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Table 1. Typical principle dimensions for various ship types 


Ship types Total shaft horsepower Length Tonnage 
Patrol Hydrofoil 18,000 shp 145 ft. 240 
Frigate 37,417 441 3880 
Destroyer 78,555 510 6395 
Cruiser (non-nuclear) 82,462 556 8872 
Cruiser (nuclear) 77,460 590 9487 
Carrier (non-nuclear) 280,000 1050 80,000 
Carrier (nuclear) 270,000 1100 73,000 
Battleship 212,000 887 58,000 
Submarine (nuclear, attack) 

15,000 292 4640 
Submarine (nuclear, strategic missiles) 

16,000 410 7880 


—=enmenmwremwremreamw emer eamwremr ee ewe ewe ee ss oS ee es awe ess as eS aw eS aS a eS SP eS a ee ew ewes ew ew ess eB awe eawrawrear saws awe HY§ aweas asa away =] ay =] = 


Source: Jane’s Fighting Ships 1985-86, Jane’s Publishing 
Company, Ltd., London, edited by John Moore. 


These numbers represent the geometric mean of several 
classes of ships within a type and should not be taken to be 
those of a particular ship. Their value lies in the ap- 
preciation of the differences between various ship types. 


An apt weight comparison would be that a typical forty foot 
Sailboat might displace fifteen tons. 


The U. S. Navy has other ship types besides those listed 


above. They include amphibious warfare ships, submarine and 
destroyer tenders, and fleet oilers, and supply ships. 
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rigidly controlled and the operators of the propulsion plant 
must be nuclear trained. Nuclear plants are used on 
cruisers, submarines, and aircraft carriers. 

Diesel propulsion plants have high weight to volume 
ratios, making their use costly in weight dollars, but 
cheaper in volume dollars.5 They are very noisy, ruling out 
their use in antisubmarine warfare ships. Although very 
reliable mechanically, they require much tinkering and 
tuning. Their specific fuel consumption® is among the 
lowest of all the plants. Medium speed diesel engines are 
not commonly manufactured in the 25,000 hp range, which 
means diesel plants cannot be used in ships requiring high 
shaft horsepower. They are typically used in smaller ships 
as cruise engines and in amphibious warfare ships as main 
propulsion. (Amphibious warfare ships typically have lower 
top speeds than frigates or destroyers. ) 

Gas turbine propulsion plants seem to have many good 
points. They are reliable, quiet, relatively low weight, 
and come in power ranges that are useful in ships ranging 
from 300 ton hydrofoils to 8000 ton destroyers. They 
require large amounts of volume for intake and exhaust duct- 
ing, but this is not a great disadvantage. Their fuel 
economy is not as good as other plants, but this is not 
intrinsic to the gas turbine engine. It is a fault of the 
operating method; gas turbine plants have mostly been built 
with mechanical transmissions. Usually one or two engines 
are coupled to each shaft. If the ship is proceeding at 
high speed, the gas turbines are operating at their full 


5. When designing a ship, total ship cost is monitored by 
the use of marginal cost factors. Every additional cubic 
foot or ton of weight added has a marginal cost associated 
with it. When a ship’s total cost is constrained, design 
changes that add cost are discouraged or must be offset by 
the reduction of other systems’ weight or volume. 


6. Specific fuel consumption is the ratio of pounds of fuel 
burned per horsepower-hour. SFC=lbf /hp-hr. 
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load design point and are relatively fuel economic. Good 
fuel economy is not usually realized, however. In the main, 
the ship proceeds at a cruise speed, using one gas turbine 
for each shaft, and the gas turbines operate at about half 
power. Specific fuel consumption rises rapidly as gas tur- 
bine power level drops, which makes for inefficient opera- 
tion. 

Usually, the above propulsion plants have mechanical 
transmissions. This means the main engines, whether they 
are diesels, gas turbines or steam turbines, are mechani- 
cally connected to the shafts and propellers. There is 
usually a reduction gear between the engine and shafting. 
These gears are large, very heavy, and expensive. To 
provide an idea of size, the largest, or “bull” gear ina 
typical locked-train double-reduction gear is about seven 
feet in diameter. The reduction gear must be placed in-line 
with the shafting, thereby using more of that prime ship 
volume. Some mechanical transmissions have cross- 
connections between shafts, but this is not common. 

Electrical transmissions are characterized by prime 
movers of any type providing power to generators. The out- 
put electricity is conditioned and sent to propulsion motors 
via a distribution network. Cross-connection is done with 
switches and breakers. There can be a mechanical reduction 
gear if it is desired to operate the propulsion motors at 
higher than propeller rotational speeds. Direct-drive 
motors may also provide the desired propeller rpm, e.g., by 
controlling the field current ina DC motor. The propulsion 
motors can be very near the propeller, i.e. aft, eliminating 
the long runs of shafting associated with a mechanical 
transmission. 

Naval ship propellers are of two types, controllable or 
fixed pitch. The pitch of a propeller is the distance the 
ship moves forward in the water for one turn of the propel- 
ler. <A fixed pitch propeller has this characteristic dis- 
tance the same at all times. A controllable pitch propeller 
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can vary this distance by changing the angle of attack (the 
angle at which the blade slices through the water) of the 
propeller blades, including reversing the blade so that the 
ship moves astern. Controllable pitch propellers are prac- 
tically required for propulsion plants that have non- 
reversing main engines, such as diesel and gas turbine 
plants.? Steam plants can reverse their propellers and 
shafts by use of an astern turbine, albeit much more slowly 
than a ship with a controllable pitch propeller system. 
Fixed pitch propellers have a slight advantage in efficiency 
(1-3%) over the controllable ones. This is mostly due to 
the large propeller hub required for varying the blade angle 
of a controllable pitch propeller. Quick reversal of shaft 
direction or propeller pitch means quick ship braking and/or 
ship reversal. This ship quickness is mandatory for an- 
tisubmarine operations and safe navigation. For example, 
the ability to stop "on a dime" may be important in a 
crowded sea lane, where a small wooden sailboat has the 
right of way over a powered naval vessel. 

Electric drive seems to combine the best of all the 
propulsion plants. It has all the advantages of a conven- 
tional prime mover plus the advantage of electrical cross 
connection and better arrangements. In the cruise scenario 
above, the electric drive ship could have both shafts 
operating from one gas turbine engine. That engine would be 
coupled to an electric generator which would produce enough 
power to run the motors that turn each shaft. The drive 
motors would be placed at the stern of the ship, near the 
propellers, on the same level. The long runs of shafting 
would be replaced by electric cable, which is smaller, 
weighs less, and can be placed in non-prime real estate. 


7. Ships with non-reversing prime movers can also have a re- 
versible reduction gear with a fixed pitch propeller instead 
of a controllable reversible pitch propeller. This is new 
technology for the United States and only the latest naval 
ship design, the DDG51, has a reversible reduction gear. 
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Cable is also in many cases cheaper than shafting, espe- 
cially to repair. Since one gas turbine would provide power 
to both shafts, it would operate at a higher power level and 
would be therefore more fuel economic. A typical propulsion 
plant might consist of three gas turbines with three gener- 
ators. Most ships have an even number of prime movers be- 
cause mechanical shaft cross connection between shafts is 
not often used and each shaft in a mechanical transmission 
ship requires the same number of prime movers to balance 
loading at high power levels. The extra prime mover 
requires a lot of weight and volume. An advantage of 
electric drive is that it becomes possible and perhaps 
desirable to use an odd number of prime movers. Each of the 
two shafts would have one propulsion motor. The heavy 
reduction gears could be replaced by the motors, which would 
have an infinitely variable reduction ratio. The control- 
lable pitch propeller system so far required by this gas 
turbine ship would be replaced by the cheaper, slightly 
smaller, and far less complicated fixed pitch propeller, ad- 
ding a small efficiency gain. The hydraulic system used to 
vary blade angle would be eliminated. Ship braking and 
reversal would be accomplished by electrically controlling 
the motor rotation direction, combined with energy dissipa- 
tion through the use of resistor banks. 

A disadvantage of this arrangement would be the high 
weight of the propulsion motors. They would be special 
designs and have a high capital cost. Hopefully, the high 
weight of the motors would be offset by the reduction in 
shafting and fuel weight and the possible elimination of the 
reduction gears. The high cost would be made palatable by 
the savings in fuel over the life of the ship. A Life Cycle 
Cost comparison of various propulsion plants, including 
electric drive, is available in reference two. 

The change to electric drive would likely be accom- 
panied by an overall decrease in propulsion plant weight and 
volume. The ship could be smaller and lighter, and would 
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require less onboard fuel for the same endurance range (the 
distance the ship can travel without refueling). Since less 
volume would be required for the fuel, the ship could be 
smaller and lighter. Since the ship would then be smaller 
and lighter, less horsepower would be required to achieve 
the same top speed. Since less horsepower would be 
required, the ship could be smaller and lighter. This is an 
example of the design spiral that would result ina smaller, 
lighter, cheaper, more risk-free ship. An example of this 
type of ship improvement is given in reference one. There 
is a limit on ship improvement, usually due to the non- 
propulsion systems or payload. One cannot make an ocean 
crossing missile ship the size of a small yacht. 

So why are not all Navy ships electric drive? They are 
not largely because the technology has not existed ina 
usable, fully developed, and manageable form. Because of 
the high cost of naval ships (a small one may cost $350 
million) and the lack of experience with current electric 
drive technology, the Defense Department is reluctant to 
build large electric drive ships. There have been electric 
drive ships, including five battleships with 21 MW shaft 
output and two aircraft carriers with 135 MW shaft output. 
Over 160 escort vessels were built during the Second World 
War with turboelectric or diesel-electric drives ranging 
from about 4.5 to 9 MW.8 A new class of ocean surveillance 
ship, the T-AGOS 19, is being built with diesel-electric 
propulsion, but it is only a 3500 shaft horsepower (shp) 
ship. 

Electric drive was replaced by conventional mechanical 
transmission plants after World War Two because of the com- 
petition afforded by improved gear cutting methods. Double- 
reduction locked-train gear transmissions became the 
standard. Since the electric drive ships all had non- 


8. Doyle, T. J. and Harrison, J. H., “Navy Superconductive 
Machinery Program," Trans. SNAME, 1978, p. 20-1. 
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superconducting, air-cooled motors and generators, they had 
higher weight and increased space requirements and suffered 
in comparison with the mechanical transmission ships.3 The 
importance of the improvement in power electronics must be 
mentioned. World War Two ships did not have the advantages 
afforded by those electronics. 

Integrated electric drive propulsion must also be men- 
tioned. This propulsion plant is the same as any of those 
discussed above, except that ship service electrical power 
is derived from the main propulsion plant, usually by taking 
power off the reduction gears or main propulsion generator. 
Power conditioning equipment, such as a cycloconverter, is 
needed to "clean up” the power and change it to fixed 
frequency for use in other equipment. Variable speed con- 
stant frequency equipment and concepts embody the integrated 
electric drive concept. The U. S. Navy has investigated 
this in some detail.1° 

Some requirements of electric drive may be viewed as 
disadvantages. The power from the electric generators has 
to be conditioned to provide frequency control of the 
propulsion motors. The power conditioners add weight and 
volume to the overall system, as well as reducing the 
efficiency of the transmission. Braking resistors, used to 
dynamically and quickly slow the propulsion motors, add more 
weight and volume to the system. There also may be a high- 
frequency radiated noise signature associated with alternat-— 
ing current systems that may be deleterious to the mission 
of the ship. 

The research done to date has not explored specific 
motor types in detail. How can it be decided whether to put 


a synchronous, inductive, permanent magnet or other motor in 


10. Robey, Stevens, and Page, “Application of Variable Speed 
Constant Frequency Generators to Propulsion-Derived Ship 
Service," Naval Engineers Journal, May 1985. 
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the electric drive system? What makes the ship system 
"best"? The effect of each motor type on the ship system 
has not been analyzed. The electric transmissions used in 
the current version of the Advanced Surface Ship Evaluation 
Tool (ASSET), a ship design computer program written for the 
United States Navy by Boeing Computer Systems, Inc., are 
generic combinations of AC and DC motors and generators, 
using rough estimates of weight, volume, and efficiency. 
Ship designs more involved than the feasibility level need 
detail on just those items. 

Motor design is a well known subject and there are many 
texts on the subject. The use of motors and pertinent tech- 
nologies in a ship as a propulsion method is discussed in 
Greene, Powell, and Gripp [3]. The advantages of electric 
drive include flexibility of arrangement, controllability, 
variable reduction ratio, reliability, and provision of ship 
service power from the main bus. Jolliff and Greene [4] go 
on to propose a specific water-cooled Advanced Integrated 
Electric Propulsion Plant (AIEPP) for a frigate/destroyer- 
sized ship. They discuss the essential characteristics of 
such a plant, establish the feasibility of the drive system 
and identify the method to technically demonstrate the sys- 
tem. Acker, Greene, and Jolliff [5] present several model- 
ing techniques and scaling relationships that allow estima- 
tion of volumes and weights of propulsion motors and gener- 
ators, solid state power conditioners, electrical 
switchgear, and associated electrical propulsion systems 
components as functions of propeller shaft power. A case 
study of AIEPP is given in the paper by Kastner, Davidson, 
and Hills [6]. 

Simulation of electric motors and associated systems is 
a popular topic. Many persons have done work in this area, 
from the micro-consideration of high frequency inductance 
changes to the more macro-consideration of hunting tran- 
sients, etc. Smith, Stronach, and Tsao [7] model a complete 


electromechanical marine drive system while Smith, Stronach, 
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Tsao, and Goodman [8] concern themselves more with a marine 
power system, including pump drives. Nonlinearities and 


operational transients are addressed. 


1.2. Optimization 


Optimization is the process of making a system, sub- 
system, or idea the best it can be. "Best" is defined by an 
“objective function," a measure of what is optimum. For ex- 
ample, an optimal manufacturing process may produce the max- 
imum number of units at the lowest cost. The objective 
function would combine units-produced with cost in an equa- 
tion that could be analyzed to find the proper production 
level. The output of the objective function is a single 
scalar measure of “goodness.” It may be difficult to repre- 
sent complicated processes with only one number. 

Optimization can be performed on a global or subor- 
dinate basis. The optimum motor might be the one that has 
the highest efficiency, even if that efficiency was achieved 
by designing a very large and heavy motor. The sub-system 
(the motor) has efficiency as its objective function. The 
ship in which the motor is to be placed may be optimum when 
its overall weight and volume are the lowest (ignoring cost, 
for example). A large, heavy motor, then, may not be op- 
timum for the ship, even if it is very efficient. A good 
case study of motor optimization is the EPRI report authored 
by Fuchs, et al. [9] 

“Optimization” can be an ill-defined term but there are 
fairly well defined methods of achieving it. Linear 
programming, Markov modeling, and Monte Carlo schemes are 
examples of these methods. The accessibility of high-speed 
digital computers has made multiple random excursions in a 
multi-variable space a much easier way of finding the 
“optimum” solution, provided an objective function and con- 
straints can be devised to describe the problem. This 


method of random excursions (Monte Carlo scheme) and ex- 
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amples of it are among the methods described in references 
ten, eleven, and twelve. 

Monte Carlo schemes take their name from the action of 
the roulette wheel in the gambling casinos of Monte Carlo. 
Around and around the wheel goes, stopping on random nunm- 
bers. If the computing power is available, this is an ac- 
ceptable method of exploring a large variable space. It can 
be much quicker than looking at every possible permutation 
of all variables. 

The steepest-descent scheme is so named because op- 
timization moves down the sharpest gradient of the objective 
function. From a valid design point, random steps are taken 
in every variable and the design point is moved "downhill" 
toward the objective function over the steepest slope. This 
is different from the “drunkard’s walk," where the random 
steps are only evaluated on whether or not the objective 
function’s output has improved, not if it is improving at 
the fastest possible rate. 

Optimization is almost always subject to constraints. 
In the previous manufacturing process, warehouse space may 
be limited, so only a certain number of units may be stored. 
This could act to limit production. For motor design, con- 
straints include maximum rotor tip speed, maximum current 
density, minimum power output, etc. All constraints should 
be combined with the objective function to yield a 


“constrained optimization.” 


1.3. The objective function 


Optimization is not possible without an objective func- 
tion. It may be very difficult to devise a good objective 
function for a complicated system such as a ship. It may be 
even harder to find one for a sub-system of that complicated 
system. There are very many characteristics that could be 
optimizing variables, and assembling them into one objective 


function with all the constraints is not easy. 
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Even with a properly defined objective function, it may 
be difficult to choose among designs that result from the 
constrained optimization. For example, a low-volume and 
low-weight motor may have poor efficiency. A very efficient 
motor may also be large and heavy. If all three elements 
are important, which is the best motor? Deciding between 
competing designs has been the subject of various papers, 
one of which is by Schweppe and Merrill [13]. In that 
paper, the authors suggest the use of "knee curves," saying 
that the essential characteristics of a multiple attribute 
tradeoff can be plotted on a series of x-y graphs. Uncom- 
petitive designs are easily discerned and discarded. The 
decision process can be limited to only those designs that 


are competitive. 


Table 2. Optimizing characteristics for ship and transmis- 


sion 


Weight 

Volume 
Efficiency 

Cost 
Reliability 
Maintainability 
Commonality 
Manning 


The above table lists many of the possible optimizing 
characteristics for the ship and its propulsion sub-system. 
Manning estimates are typically based on historical data and 
do not indicate that the baseline and variant ships will 
require a significantly different number of personnel. Com- 
monality measures the use of the same equipment in other 
ships. Since there are no other electric drive ships at the 
power levels used in this thesis, commonality is not an 
issue. It is very difficult to discern maintainability and 
reliability differences between designs that are as close as 
the machine designs of this thesis, so these two charac- 
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teristics were not made part of the objective function. 

Cost is a measure that should be part of this sort of op- 
timization. The cost of the ship system, quantified in the 
ship displacement, was used in the final recommendations for 
the transmission sub-system. In the case of permanent mag- 
net machines, the relative costs of magnet material and mag- 
net steel were included in the objective function. 

Weight, volume, and efficiency were made an explicit 
part of the objective function for the computer design of 
the machines. An Effective Weight was calculated for every 
machine design. The design with the lowest effective weight 
was the "best" within its class. The generic objective 


function is 


Effective Weight = weight + kex(i-efficiency) + kv*volume 


where ke and kv are weighting factors for efficiency and 
volume, respectively. The weighting factors were obtained 
from changes in ship displacement for marginal changes in 
efficiency and volume of the transmission. They were 
modified to reflect the actual designs resulting from the 


process. 
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Chapter Two. General Considerations 


Only steady-state behavior of electric machines was 
modeled. The modeling of dynamic behavior is very difficult 
and was not viewed as being within the design problems posed 
by this thesis. The changes in machine design necessary to 
solve dynamic instabilities, etc., are much smaller than the 
approximate nature of the algorithms used here. 

All derivation work was performed without specifying 
the number of winding turns or the number of rotor or stator 
slots. The only exception to this was the case of induction 
machines, where an arbitrary number of rotor slots was 
selected. This selection was necessary for the calculation 
of the equivalent circuit components. The number of rotor 
slots chosen, 71, was a number designed not to induce pole 
harmonics. Since no turn numbers were specified, units in- 
clude volts-per-turn, ampere-turns, and impedance-per-turns- 
squared. Power is measured in watts. 

All machines used as their synchronous frequency the 
maximum allowed by the particular combination of pole pairs 
and shaft rpm. Developmental work showed that the optimiza- 
tion algorithm converged to the highest frequencies, so the 
algorithm now starts at the highest possible frequencies. 

Up to six pole pairs were used in the higher rpm 
machines and up to 25 pole pairs in the 180 rpm machines. 
Diminishing improvements in volume, weight, and efficiency 
show up at half these limits. 

The random number generator was taken from Kelley and 
Pohl [14], with one change. After every run of each 
program, the random number generator seed is stored. This 
means the sequence of pseudo random numbers is not repeated 
until the full range possible has been used. For the 
machines, it gives differences at every run and means the 


multidimensional variable space is more fully explored. 


24 





2.1. Optimization method 


The chosen optimization method is a combination of the 
Monte Carlo and steepest-descent schemes. A design point is 
established by randomly selecting machine geometric 
parameters, subject to constraints. Ten random steps are 
taken around the design point, in all variable directions. 
The effective weight of each random step is evaluated and 
compared with that of the design point. The best of the 
eleven is designated the new design point. More random 
steps are taken, and the process continues until no more im- 
provement is seen in effective weight. At that point, the 
size of the random steps is halved, and the process repeats 
itself, with the step size continually halved (up to ten 
times). The best effective weight is the index to the best 
design. The number of original design points used in any 
particular run of a program is under user control. If ten 
original design points are desired, the algorithm will look 
at over a thousand designs. 

The purpose of having original design points is to 
start the optimization process in different sectors of the 
multidimensional variable space. In this fashion, the op- 
timization process zeroes in on several local "best" points. 

The variables that are randomly selected include stator 
current density, rotor radius, air gap dimension, stator 
slot space factor, and rotor slot space factor. The back 
iron dimension (the iron behind the stator teeth) is sized 
to handle a saturation level of flux. The stator slot depth 
is originally sized as a random fraction of the back iron 
dimension, and the rotor slot depth is originally a random 
fraction of the rotor radius. 

Only steady-state behavior was modeled in this algo- 
rithm. Dynamic modeling may or may not show different op- 


timum configurations for machine types. 
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2.2. Constraints 


The constraints placed on the optimization process are 
listed in the following table. The most difficult con- 
straint to satisfy while still achieving a valid design was 
the rotor current density constraint in induction machines. 
Only a few valid designs were achieved in induction machines 
using the above algorithm, leaving some question about the 
application of the algorithm in the case of induction 
machines. Only those induction machine designs in which 
there was reasonable confidence were included in the thesis 


analysis. 


Table 3. Optimization constraints during machine design 


Minimum air gap flux density 1.05 tesla rms 
Maximum (saturation) flux density 1.5 tesla rms 
Maximum rotor radius 2.0 meters 
Maximum rotor tip speed 200.0 meters/sec 
Stator and rotor space factor 0.35 
Maximum rotor slot depth 33% of rotor radius 
Maximum synchronous reactance: 
synchronous machines 2.0 per unit 
permanent magnet machines 3.0 per unit 
Power factor 0.8 


The magnet steel chosen was 26 gauge M19. Its magnetic 
properties were found in USX technical data [15]. It has 
been observed that saturation flux levels in electric 
machines occur first where the area perpendicular to the 
flux path is the smallest. If the back iron dimension is 
made appropriately large, this saturation will first occur 
in the teeth, as is desirable. Accordingly, the back iron 


dimension was set to 


where Br is the radial air gap flux density, r is the rotor 
radius, Bsat is the saturation flux density, and p is the 
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number of pole pairs in the machine. This equation is 
derived from Gauss’ Law. 

If electric machines are to be installed in a ship, th 
78 ey obviously will need to fit into the space designated 
for them. An electric machine with a two meter rotor radius 
is at least thirteen feet in envelope diameter. This is a 
very large machine to install in a machinery space where 
volume is already at a premium. The rotor radius limit de- 
scends from the physical ability to fit an electric machine 
in a ship. 

The tip speed limitation represents the physical limit 
on material strength with regard to the rotor conductors. 
Rotor conductors may break free from the rotor at higher 
tangential velocities than this limit. The magnitude of the 
limit was taken from a tip-speed-limited, 3600 rpm, two-pole 
turbogenerator, and was verified against standard Navy 
design practices. Several as-built electric machines were 
analyzed and this number seemed to fit their characteristics 
well. The tip speed limit arises when choosing a rotor 
radius (given a particular frequency and number of poles), 
and is less stringent a constraint than maximum rotor 
radius. 

Thermal considerations are often extremely important in 
machine design. The heat build-up in electric machines has 
led to many cooling schemes over the years, including 
natural convection, forced air cooling, and hydrogen cool- 
ing. One of the latest methods is liquid cooling of the 
stator and rotor conductors through cooling passages through 
the copper itself. This has been made possible by better 
de-ionizing methods for cooling fluid and better rotating 
seals for the rotor. Naturally, the cooling passages and 
insulation limit the amount of copper area in a slot cross 
section. The copper area in a typical conductor bar was 
measured and found to be about thirty-five percent of the 
bar cross section. This number was used for the stator and 


rotor slot space factors. 
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Other thermal considerations must be made for permanent 
magnet materials, which suffer a degradation in flux as tem- 
peratures rise. Flux loss rises slowly with increasing tem- 
perature until about 100°C. Above 100°C, flux loss is more 
rapid. An assumption in the design of these machines is 
that there will be sufficient cooling in the operating space 
to limit ambient temperature to about 80°C. This, combined 
with the machine liquid cooling, should keep flux toss to a 
minimum. Transmission lines were assumed to function satis- 
factorily at the same temperature. 

Insulation also has a thermal rating. No insulation 
class was specified in this thesis but a typical insulation 
used in electric generators by the Navy is Class F. For 
this class, a permissible rise of 100°C over an ambient tem- 
perature of 50°C is standard, but lesser insulation classes 
must run cooler. If a machine must be designed with a les- 
ser insulation class, the consequent lower temperatures may 
result in a quieter machine and longer machine and insula- 
tion life. It probably will be larger than a machine with a 
greater class of insulation. The temperatures quoted above 
are at hot spots, not in the bulk of the machine.!1! 

Along with the reduction in copper area for liquid 
cooling, a maximum current density was imposed. The copper 
losses, in the form of heat, have to be removed by the cool- 
ing fluid. There is a tradeoff between the size of the 
cooling passages, the allowable current density, and the 
rating of the machine. Twelve million amperes per square 
meter equates to forty amperes in a twelve gauge copper 
wire. 

Rotor slots were constrained to no more than one third 
of the rotor radius. Some reasonable shaft diameter is 
required to transmit the mechanical torque. Stator slots 
were allowed to grow as needed to meet the stator current 


11. Private communication, D.F. Schmucker, Naval Sea Systems 
Command, Code 56231, 11 December 1986. 
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density limit. 

Synchronous reactance limits were taken from as-built 
machines. A power factor of 0.8 was used for all machines, 
though one researcher has indicated a power factor of 1.0 


might be best for permanent magnet machines.!2 


2.3. Geometric considerations 


End turns were modeled as described in Appendix B. 
While not exactly as machines are constructed, this model 
gives reasonable results. A length allowance equal to one 
rotor diameter on either end of the active length of the 
machine was made to allow for containment of the end turns. 

Fractional slot pitches were not considered. A stator 
winding pitch of 0.8 was assumed, resulting in the elimina- 
tion of the fifth harmonic from the steady-state output 
waveforms. 

Once weight and volume were calculated, an extra ten 
percent was added to allow for the frame and foundation of 
the machine. The calculated weight included an additional 
three percent of the rotor weight to allow for bearings. It 
is the final envelope weight and volume that were used in 


the decision process. 


2.4. Efficiency and losses 


The general equation for efficiency is 


minpwr 
Pestle @eN CY a= ieee = = — o oa > So ---- se 
minpwr + ph + pe + i2dr + i2rr 
where minpwr is the minimum mechanical power expected of the 


machine, ph is the hysteresis loss, pe is the eddy current 


12. Robey, H.N., “Permanent Magnet Machine Technology 
Assessment," DTNSRDC Report TM-27-80-87, September 1980. 
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loss, i2r is the stator copper loss, and i2grr is the rotor 
copper loss. This formulation is for a motor, but the ef- 
ficiency calculated will not be significantly different if 
the machine is a generator. 

Hysteresis and eddy current losses arise from currents 
circulating within the magnet steel that forms the rotor and 
Stator. They are two different mechanisms and depend on the 
metallurgy of the steel. 

Eddy currents are a result of the time-varying magnetic 
fields within the machine, and they oppose the change in 
flux density within the machine. Eddy current losses in- 
crease as the square of the electrical frequency of the 
machine and also as the square of the peak flux density. 

One method of lessening eddy current losses is to use thin 
laminations to build up the rotor and stator. If the var- 
nish used on the laminations is sufficiently insulating, the 
eddy currents are limited to azimuthal circulation. Axial 
circulation is practically zero because of the small lamina- 
tion thickness. 

Hysteresis losses are inherent to magnetic materials, 
and are proportional to the total volume of the material, 
the area of the hysteresis loop, and the machine electrical 
frequency. 

USX has developed equations to calculate eddy current 


and hysteresis losses. They are 


0.01445 6 £ Br Hc 0.4818 fl Bre t? f2 


In these equations, @ is the hysteresis loss factor (the 
ratio of the actual hysteresis losses to the area of a 
square hysteresis loop passing through Br and Hc), f is the 
frequency in Hertz, Br is the residual induction in 
kilogauss, Hc is the coercive force in oersteds, D is the 
density in grams per cubic centimeters, rho is the electri- 
cal resistivity of M19 in microhm-cm, fl is the anomalous 
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loss factor. The losses are in watts per pound of material. 
The numerical factors at the beginning of both equations 
were altered to reflect the use of SI units. The factors @, 
NM, Br, He, rho, and D change with the type of magnetic 
material used. The thickness of the laminations, t, is 
0.014 inches for 26 gauge steel. These equations do not 
reflect some variations caused by differences in silicon 
content and differences in processing treatments leading to 


variations in grain size and crystallographic texture. 
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Chapter Three. Synchronous machines 


Synchronous machines operate because of an interaction 
between stator and rotor flux waves. The rotor flux wave is 
developed by a field winding. The stator flux wave is 
developed by an armature winding. These two waves try to 
align themselves, which is how the machine action is 
produced. In the case of a motor, the armature wave is 
“rotating” around the periphery of the stator bore because 
of the 120° separation between the three phases. A rotating 
action ensues. In the case of a generator, the rotation is 
provided by a prime mover, such as a gas turbine, and volt- 
age is induced in the stator phases. 

Synchronous machines operate at a steady-state shaft 
speed specified by the number of poles and the electrical, 
or synchronous, frequency. This synchronous speed is main- 
tained despite changes in load. This feature makes 
synchronous machines attractive for applications where speed 
control is important. Shaft rpm = (120°frequency)/(poles). 

A derivation of the equations of synchronous machines 
and the computer modeling program are presented in 


Appendix B. 


3.1. Assumptions 


The rotors of synchronous machines may exhibit 
saliency, or may be smooth cylindrical rotors. The dif- 
ferences in properties and parameters among salient and 
round-rotor machines amount to only a few percenti3. The 
approximate nature of the modeling means the saliency ef- 
fects will not be important. Therefore, no special provi- 
sion for salient rotors were made. 

The ships used in this thesis have displacements of 


13. Fitzgerald, Kingsley, and Umans, Electric Machinery, 
McGraw-Hill, New York, 1983. 
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about 5000 LT ( 1 LT = 2240 lbs). Their fuel load is deter- 
mined by the required range and transmission efficiency at 
endurance speed. It was first estimated that a one-percent 
increase in transmission efficiency would reduce total ship 
displacement by about 89 LT. The efficiency factor (ke = 
90,000 kg/percent) corresponds to this 89 LT change in ship 
full load displacement. (89 LT x 2240 lbs/LT x 2.205 kg/lb) 
When this produced machines with efficiencies about 95%, it 
was doubled to 180,000. Obviously, this factor may be ad- 
justed to any level. The volume efficiency factor (Kv = 
1286.1 kg/m?) corresponds to the density of a LM-2500 gas 
turbine module. These factors were used throughout the 
thesis. 


3.2. Machine description 


Machines with shaft speeds of 1800, 2400, 3000, 3600, 
and 7200 rpm, with the number of pole pairs varying from one 
to six, were modeled. Also, 180 rpm machines using from one 
to twenty-five pole pairs were modeled. This provided a 


good coverage of the variable space. 


3.2.1. Efficiency 


Synchronous machine efficiency at full load was about 
98.5% for the higher rpm machines, while the 180 rpm 
machines hovered around 93% to 94% efficiency. The number 
of pole-pairs seemed to have little effect in the high rpm 
machines, but there was an "arch" in the efficiency curve of 
the 180 rpm machines, peaking at 97% with 36 poles. Though 
not fully understood, the 24- and 26-pole 180 rpm machines 
had very low efficiencies. Generally speaking, machine ef- 
ficiency was higher when rpm was higher. This was an ex- 
pected result. Off-design-point efficiency was good for the 
higher rpm machines, but bad for the 180 rpm machines. (See 


the discussion in Chapter Seven. ) 
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3.2.2. Weight and volume 


Weight and volume increased almost linearly with the 
number of pole-pairs in the 180 rpm machine. This is almost 
certainly a function of tip-speed limitations, as the maxi- 
mum tip speed in a machine is a function of the number of 
poles in the machine. When the rotor radius is limited, the 
machine must grow in length to develop enough torque. The 
weight and volume of these machines were much higher than 
for the higher rpm machines. 

The higher rpm machines saw significant decreases in 
weight and volume when the number of pole-pairs increased 
from one to two. There also was an observable increase in 
weight and volume as the number of pole-pairs further in- 
creased. As rpm increased, the machines became smaller and 
lighter. 
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Figure 6. 180 rpm synchronous machine weight 
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3.3. Machine characteristics 


The following tables give machine characteristics for 
many of the higher rpm and the 180 rpm synchronous machines 
designed for this thesis. 

The stator slot factor tended to increase to the limit 
of 0.75, while the rotor slot factor moved around the value 
0.58 quite a bit. This demonstrates a partial limit on the 
depth of the stator slots (to the same dimension as the back 
iron depth). The overall diameter was limited and more 
stator slot area was needed to develop the required power. 
Similarly, the stator current density converged to its cool- 
ing limit. 

The longest higher rpm machine was 5.42 meters, while 
the largest overall diameter was 0.85 meters. Machines of 
this size will cause no difficulties when placed in the 
machinery spaces of most ships. The 180 rpm machines are 
larger, typically less than 6 meters long (discounting one 
16 meter machine) and 1.8 meters in diameter. They are also 


good candidates for ship systems. 
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Table 4. 


nuaber of pole pairs 
power, hp 

efficiency factor 
voluae factor 

shaft rpae 

stator current density 
synchronous frequency 
rotor radius 

Qap dimension 

back iron depth 

stator slot depth 

rotor slot depth 

Stator slot factor 
rotor slot factor 
envelope voluae 
envelope weight 
hysteresis loss 

eddy current loss 
stator copper loss 

full load efficiency 
active length 
full load current density 
no toad current density 
xs/turns-squared, p.u. 
internal volts/turn, p.u. 
overall! length 


overall! diaaeter 


Characteristics 


I 

29779 
180000 
1286. 1 
1800 
8.00E+06 
30 
0.1720 
0.0158 
0.1204 
0.0772 
0.0401 
0.693 
0.690 
2.781 
20835. 40 
16926. 49 
3033.2 


of 1800 rpm synchronous machines 


2 

oa713 
180000 
1286.1 
1800 
1.20E+0/ 
40 
0.2077 
0.0440 
0.0727 
0.0722 
0.1859 
0.750 
0.447 
1.522 
9827.96 
10677.39 
3826.7 


279834.4 374646.9 


0.985 
4,670 
1.19€&+07 
4, 40E+04 
Neth 
2.71 
9.42 
0.77 


0.980 
1.794 
1.06E+07 
6. SPE+04 
0.86 
1.66 
2.80 
0.79 


3 

25779 
180000 
1286.1 
1800 

1. 20E+07 
70 
0.2603 
0.0250 
0.0607 
0.0607 
0.0747 
0.750 
0.687 
1.509 
8766.09 
17111.27 
9198.8 
293339.9 
0.984 
1.498 


: rn) 

29775 = 23773 
180000 180000 
1286.1 1286.1 
1800 1800 
1.20E+07 1. 20E+07 
120 150 
0.2905 0.3267 
0.0124 0.0090 
0.0908 06,0457 
0.0508 0.0457 
0.1333 0.0907 
0.7479 0.750 
0.699 0.449 
1.541 1.697 
9186.03 8661.17 
19469.83 $1370.98 
13955.6 28107.7 
164289.8 173116. 1 
0.990 0.988 
1.513 1.348 


L.17E+07 4.58E+06 8. 34E+06 


7.00E+06 
0.87 
1.67 
2.64 
0.81 


2.43E+06 4. 24E+06 
1.11 1.20 
1.89 1.97 
2.72 2.69 
0.81 0.85 


) 

20775 
180000 
1286.1 
1800 
1.20E+07 
180 
0.3189 
0.0092 
0.0372 
0.0372 
0.0289 
0.748 
0.428 
1.709 
8853. 21 
43892.02 
47191.0 
269901 .7 
0.982 
1.741 
1.B9E+07 
1.19€+07 
0.78 
1.60 
3.05 
0.80 





Table 5. 


nuaber of pole pairs 
power, hp 
efficiency factor 
voluae factor 
shaft rpa 
stator current density 
synchronous frequency 
rotor radius 
gap dimension 
back iron depth 
stator slot depth 
rotor slot depth 
stator slot factor 
rotor slot factor 
envelope voluse 
envelope weight 
hysteresis loss 
eddy current loss 
stator copper loss 
full load efficiency 
active length 
full load current density 
no load current density 
xs/turns-squared, p.u. 
internal volts/turn, p.u. 
overall length 


overall diaseter 


Characteristics 


I 

29775 
180000 
1286.1 
2400 
1.19&+07 
40 
0.1563 
0.0273 
0.1094 
0.0929 
0.0416 
0.621 
0.676 
1.411 
11619.57 
11981.17 
2862.6 
444605.4 
0.977 
Deon 


of 2400 rpm synchronous machines 


2 

25779 
180000 
1286.1 
2400 
1.20€+07 
B0 
0.1984 
0.0164 
0.0695 
0.0695 
0.0215 
0.750 
0.749 
1.112 
7204.42 


3 

29775 
180000 
1286. | 
2400 
1.19€+07 
120 
0.2388 
0.0089 
0.0557 
0.0557 
0.0088 
0.750 
0.790 
1,133 
6477.15 


14511.75 21941.41 


6934.5 


13727.2 


427167.9 482231.4 


0.977 
1,709 


0.974 
1.992 


4 
29779 
180000 
1284.1 
2400 
1.20E+07 
140 
0.2890 
0.0135 
0.0506 
0.0505 
0.0563 
0.748 
0.432 
1.329 
6733.44 
27615.01 
26391.9 
232958. 2 
0.985 
Velaz 


Z.11E+O7 S.L1E+07 5. S9E+07 1.73E+07 
1,2BE+07 2. 11E+07 9.646406 


7,.75E+06 
2.00 
Za 
Sale 
0.77 


1.49 
2.43 
Zoae 
0.71 


1.83 
2.06 
2.04 
0.72 


1.01 
1.80 
2.36 
0.81 


s) 

25779 
180000 
1286.1 
2400 

1, 20E+07 
200 
0.3141 
0.0109 
0.0440 
0.0454 
0.1047 
0.750 
0.453 
1.436 
7106.43 


6 

25775 
180000 
1286.1 
2400 
1.07E+07 
240 
0.33597 
0.0064 
0.0396 
0.0396 
0.1040 
0.750 
0.327 
1,433 
7727.01 


$2769.75 47103.39 


391479 
143185.6 
0.989 
1.144 


67525.7 
114790.8 
0.988 
1,225 


7,83E+06 7.99&+06 
4, S5E+06 4, 32E+06 


0.92 
1.72 
2.44 
0.82 


1.07 
1.85 
2.41 
0.85 





Table 6. Characteristics of 3000 rpm synchronous machines 


nuaber of pole pairs 
power, hp 

efficiency factor 
voluae factor 

shaft rpe 

Stator current density 
synchronous frequency 
rotor radius 

gap diaension 

back iron depth 

Stator slot depth 

rotor slot depth 

stator slot factor 
rotor slot factor 
envelope voluse 
envelope weight 
hysteresis loss 

eddy current loss 
stator copper loss 

full load efficiency 
active length 

full load current density 
no load current density 
xs/turns-squared, p.u. 
internal volts/turn, p.u. 
overall length 


overall diaseter 


J 

Zale 
190000 
1286.1 
3000 
1,20E+07 
20 
0.1494 
0.0390 
0.1046 
0.0862 
0.0904 
0.730 
0.098 
1.205 
9815.49 
9074.79 
2999.6 
380289. 7 
0.980 
1.670 
1.53E+07 
6. 03E+06 
1,980 
2.93 
2.42 
0.76 


2 5 4 a) 
Ziiveeeeniia  Zalia) 920775 
180000 §=180000 180000 180000 
1286.1 1286.1 1286.1 1286.1 
3000 3000 3000 5000 

1, 20E+07 1.20E+07 1.20E+07 1. 10E+07 
100 130 200 220 
0.1890 0.2304 0.2432 0.2872 
0.0172 0.0101 0.0052 0.0068 
0.0662 0.0538 0.0426 0.0402 
0.0628 0.0537 0.0422 0.0402 
0.0879 0.0372 0.0101 0.0784 
0.750 0.730 0.790 0.790 
0.494 0.724 0.730 0.255 
0.963 0.972 0.991 1.205 
6499.77 3801.99 9725.40 6176.39 
14072.41 21089.10 31331.61 41076.93 
8405.7 19895.4 37429.9 61340.9 
212641.9 183124.5 204027.7 140446.5 
0.998 0.989 0.986 0.988 
1.657 1.365 1.990 1.306 
L.13E+07 1.386407 2.42E+07 1.40E+07 
J. 20E+06 6.09E+06 9. 63E+06 7. 44E+06 
1.41 os 1,78 1.11 
2.16 2.2] 2.01 1,89 
2.48 2.32 2.08 2.48 
0.67 0.70 0.67 0.79 


44 


6 

yeny il 
180000 
1286.1 
3000 

1. 18E+07 
300 
0.2974 
0.0088 
0.0347 
0.0347 
0.0427 
0.720 
0.756 
1.236 
6143.64 
46178. 46 
§2749,7 
134484, 2 
0.986 
1.510 

1. 06E+07 
6.96E+06 
0.49 
fea 
Dade 
0.79 





Table 7. Characteristics of 3600 rpm synchronous machines 
nuaber of pole pairs 1 2 3 4 ;) 6 
power, hp mcd meezdiies dala Zojia  2alid 
efficiency factor 180000 180000 180000 180000 180000 180000 
voluae factor 1286.1 1286.1 1286.1 1286.1 1286.1 1286.1 
shaft rp@ 3400 3600 3600 3600 3600 3600 
stator current density L.16E+07 1. 20E+07 1.20E+07 1. 19E+07 1.20E+07 1.20E+07 
synchronous frequency 60 120 180 240 300 360 
rotor radius 0.1413 0.1824 0.2191 0.2736 0.2846 0.2877 
gap diaension 0.0139 0.0120 0.0123 0.0178 0.0209 0.0124 
back iron depth 0.0989 0.0638 0.0511 0.0479 0.0398 0.0336 
Stator slot depth 0.0492 0.0638 0.0486 9.0479 0.0398 0.0336 
rotor slot depth 0.0496 0.0541 0.0572 0.0298 0.0595 0.0385 
stator slot factor 0.748 0.750 9.750 0.750 9.750 0.750 
rotor slot factor 0.735 0.409 0.928 0.679 0.691 0.654 
envelope voluae 1.233 0.814 0.870 1.065 1,127 1.102 
envelope weight 9385.14 5377.73 9216.22 4979.26 5177.33 3190.73 
hysteresis loss 14030.95 15587.02 22251.57 31042.24 36284.98 48590. 52 
eddy current loss 9028.6 11172.5 23924.3 44501.0 65021.1 104486.5 
stator capper loss 190131.9 228605.1 165930.4 251863.9 196473.3 194008.4 
full load efficiency Orche wdaga, 0,907) O.903 0.985 9.983 
active length 3.298 = 1.493) 1.372 0.870 (00975 166 
full load current density P.ASE+OG 1.94E+07 1. 31E+07 2.47€+07 1.45E+07 1.73€+07 
no load current density 3.3BE+06 7.41€+06 6.99E+06 1.61E+07 1.126407 1. 28E+07 
xs/turns-squared, p.u. 177 1.89 1.10 0.71 0.43 0.48 
internal volts/turn, p.u. Dell 2.62 1.88 1.94 1.30 fess 
overall length 3.88 Deal 2.30 2.06 AAA, Me] 
overall diaaeter 0.61 0.64 0.66 0.77 Oe? Ome 





—— 


Table 8. Characteristics of 7200 rpm 


number af pole pairs 
power, hp 

efficiency factor 
voluae factor 

shaft rpa 

Stator current density 
synchronous frequency 
rotor radius 

gap dimension 

back iron depth 

stator slot depth 

rotor slat depth 

stator slot factor 
rotor slot factor 
envelope volume 
envelape weight 
hysteresis loss 

eddy current loss 
Stator capper lass 

full load efficiency 
active length 

full load current density 
no load current density 
xS/turns-Squared, p.u. 
internal volts/turn, p.u. 
overall length 


overall diameter 


17305.86 22097. 70 
12404.5 31678.5 
149625.0 159505. 1 


L.65E+07 2.71&+07 
6, 14E+06 1. 16E+07 


Synchronous machines 


4 . ) 
(a ee OE eee 
180000 180000 180000 
1286.1 1286.1 1286.1 
7200 7200 7200 


1.20E+07 1.20E+07 1.20E+07 


480 600 720 
0.2044 0.2219 0.1932 


0.0093 0.0046 0.0055 


0.0358 0.0311 0.0225 


0.0358 0.0311 0.0225 


0.0714 0.0138 0.0292 
0.750 0.790 06.750 
0.633 0.468 0.570 
0.544 0.572 0.564 


5136.52 2950.31 $427.07 
28944.03 30808.73 30465.33 64447,34 

62239.7 88332.6 180863.2 277168.2 
74543.2 143253.5 117812.4 


0.990 0.981 0.977 
1.081 1.082 1.982 


3.92E+07 7. B1E+06 3.708407 1.72E+07 
1.46E+07 5.00E+06 1.99E+07 1.28E+07 


0.74 1.08 0.48 
1.96 1.86 1.34 
1.94 1.99 2.73 
Doe! 0.98 0.49 





Table Q. 


nuaber af pale pairs 
power, hp 

efficiency factor 
voluae factor 

Shatt rpa 

stator current density 
synchronous frequency 
rotor radius 

gap dimension 

back iran depth 
stator slot depth 
rator slot depth 
Stator slat factor 
rotor siot factor 
envelope voluse 
envelape weight 
hysteresis loss 

eddy current lass 
Stator copper loss 
full load efficiency 


active length 


full load current density 


na load current density 


x5/turns-Squared, p.u. 


internal volts/turn, p.u. 


overall length 


averall diaseter 


| 

25775 
180006 
1286.1 
180 

1. 18E+07 
5 

0.2641 
0.0331 
0.1849 
6.0592 
0.0500 
0.741 
0.558 
16.305 
1i9110.5 
10465.3 
187.5 
2.25E+06 
0.895 
14.913 
1. S5E+07 
3. b1E+06 
2.00 
2.12 
16.10 
1.08 


A7 


2 

29775 
180060 
1286.1 
180 
1.12&+07 
6 

0.3738 
0.0978 
0.1308 
0.1308 
0.1867 
0.750 
0.971 
Tver 
36500. 0 
6915.1 
247.8 
2.06E+06 
0.903 
Selda 
1.25€+07 
6. 14E+06 
1.28 
2.04 
5202 
1.47 


3 
25173 
180000 
1286.1 
180 


4 
25775 
180060 
1286.1 
180 


1.20E+07 9.25E+06 


9 
0.4588 
0.0723 
0.1071 
0.1071 
0.1934 

0.750 
0.422 
9.048 
48896.8 
9833.5 
528.6 
1.76E+06 
0.916 
2.590 
1.53E+07 
7.25€+06 
1.04 
1.83 
4.71 
1.49 


12 
0.5869 
0.0517 
0.1027 
0.1027 
0.0545 

0.742 
0.535 
11.378 
36598.0 
18267.4 
1509.4 
1. S4E+06 
0.934 
2.576 

1. 91£+07 
B.93E+06 
1.38 
2.14 
4.85 
1.65 


Characteristics of 180 rpm synchronous machines 


n) 

Panis 
180000 
1286. 1 
180 
1.20E+07 
15 
6.9907 
0.0198 
0.0771 
0.0580 
0.1079 
0.750 
0.287 
10.306 
00899.5 
23999. 1 
2150.3 
1. 186+06 
0.941 
3.708 
1.42&+0/ 
7.27€+06 
1.19 
1.96 
aa? 
1.41 








Table 10. Characteristics of more 180 rpm 


machines 


nuaber of pole pairs 
power, hp 

efficiency factor 
voluae factor 

shaft rpa 

stator current density 
synchronous frequency 
rotor radius 

gap disension 

back iron depth 

stator slot depth 

rotor slot depth 

stator slot factor 
rotor slot factor 
envelope voluse 
envelope weight 
hysteresis loss 

eddy current loss 
stator copper loss 

full load efficiency 
active length 

full load current density 
no load current density 
xs/turns-squared, p.u. 
internal volts/turn, p.u. 
overall length 


overall diaseter 


6 

25779 
180000 
1286.1 
180 

» 20E+07 
18 
0.6176 
0.0647 
0.9720 
0.0720 
0.1652 
0.790 
0.368 
11.693 
49739.8 
22385.0 
2406.8 
1.B6E+06 
0.911 
fel da 
1.74E+07 
1. 29E+07 
0.48 
[So 
4.99 
1.65 


-_ 
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] B 7 

JT Sea Pati 
180000 180000 180000 
1286.1 1286.1 1286.1 
180 180 180 
1.20&+07 1.20&+07 1.20E+07 
21 24 2] 
0.9729 0.7152 9.6943 
0.0335 0.0229 0.0138 
0.0573 0.0626 9.0540 
0.0567 0.0626 0.0540 
0.0855 0.2046 0.0740 
0.730 8.790 0.750 
0.790 0.374 *° 0.283 
10.486 12.870 12.385 
93234.1 92286.9 33189.2 
27947.6 31793.9 42718.3 
3903.7 4597.9 6889.4 
1.246+06 7.69E+05 1.14E+06 
0.938 0.960 0.942 
3.418 2.043 2.949 
1.11€&+07 7.09E+06 1.98E+07 
7.98E+06 4.18E+06 1.07E+07 
0.93 0.89 1.08 
1.39 1.69 1, 86 
0.84 2200 0.38 
1.44 1.73 1.43 


synchronous 


10 

fis is) 
190000 
1286.1 
180 

1. 20E+07 
30 
0.7276 
0.0396 
0.0511 
0.0509 
0.1420 
0.790 
0.619 
14.273 
96440.9 
42086. 1 
7941.6 
1, 11E+06 
0.943 
2.369 
9.69 +06 
7.76E+06 
0.36 
1.29 
3.44 
1.74 





3.4. Verification 


Data on a large turbogenerator, Big Sandy Unit Two, was 
available in reference sixteen; it was used to verify the 
Synchronous machine design program. Big Sandy is rated at 
907,000 kVA, power factor 0.9, at a rated voltage of 26 kV. 
When the machine parameters were input to the synchronous 
design program, it produced a machine very close to Big 
Sandy. It was judged that the design program would yield 


good results. 
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Chapter Four. Permanent Magnet Machines 


Permanent magnet machines are very similar to 
synchronous machines. The main difference lies in the 
method used to produce the field flux wave. Instead of a 
field current causing the wave, permanent magnets provide 
the flux. No exciter is used with these machines and there 


are no brushes. 


4.1. Magnet material 


Many different elements may be used to manufacture per- 
manent magnets. Past designs used ceramics, aluminum- 
nickel-cobalt-iron-titanium (ALNiCo), and samarium-cobalt 
(SmCo). However, ceramic magnets do not produce sufficient 
residual flux (see Appendix C for an explanation of terms), 
and any magnet based on cobalt is high in cost and may be in 
limited supply. Recently, magnets of neodymium-iron-boron 
(NdFeB), have entered the marketplace. None of the con- 
stituents of the NdFeB magnets are strategic materials; it 
is expected that availability and cost will improve. 

NdFeB magnets have a high Maximum Energy Product (MEP) 
that may be used advantageously by the machine designer. 
Data from Sumitomo Special Metals [17] indicates their 
NEOMAX line have MEPs as high as 37 MGOe, which is higher 
than the 30 MGOe of the SmCo magnets marketed as REC-30 by 
TDK Corporation ({18]. High MEP is not the only criteria for 
magnet selection; flux stability, cost-to-performance ratio, 
ease of machine assembly, and other characteristics may en- 
ter the decision process. This study needed the best per- 
formance of its machines, so NdFeB was selected on the basis 
of its high MEP. Thermal stability was assumed to be satis- 
factory if the thermal considerations in Chapter Two were 
met. 

Cost has been mentioned above. Magnet material is sig- 


nificantly more expensive than magnet steel and copper. 
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Magnets in the quantity used by a large production run of 
25,000 hp machines might cost as much as $120 per pound,!4 
compared with the 58¢ per pound of M19 steel.!15 Obviously, 
permanent magnet machines will be more expensive, but the 
cost of magnet materials may be made part of the optimiza- 
tion process. The degree of magnet overhang, discussed in 


Appendix C, also affects cost. 


4.1.1. Magnet cost factor 


A change to the objective function was made to incor- 
porate the cost of magnet material relative to magnet steel. 
The ratio of the above costs was taken and the result called 
the magnet cost factor, km. The objective function was 
modified to 


Effective weight = (weight + ka(magnet weight)) + ke(l - effcy) + kv(voluee + ka(eagnet volume) 


An initial value for km of 170 was used, and several 
machines designed. Then a value of km = 25 was tried. The 
machines with km = 170 indeed had less magnet material in 
them, but at a cost. The change to km = 25 resulted in a 
larger machine (23.5%) with a lower stator current density, 
20% more magnet material, and about a 1.5% increase in 
machine efficiency. That 1.5% translates to a lot of fuel 
aboard a ship, so it was decided to use km = 25. The extra 
magnet material will add about $22,000 per machine. 


4.2. Assumptions 


The largest obstacle to assembling high-power permanent 
magnet machines is their inherent residual magnetism. If 


14. Estimate by Mr. Yokokura, President of Sumitomo Special 
Metals of America. 


15. Book price for 26 gauge M19 steel from Mr. Dagg of USX. 
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the magnets possess all of their properties at assembly, it 
will be extremely difficult to place the rotor (which con- 
tains the magnets) inside the stator. The rotor would be 
strongly attracted to the iron of the stator. Of course, 
the magnets may be magnetized after machine assembly, but it 
may be difficult to achieve MEP without elevated tempera- 
tures inside the machine. The assumption is made here that 
the magnets will be magnetized prior to assembly. The 
detailed design of the machine will have to include con- 
sideration of the jigs and fixtures necessary for assembly. 

The only other assumption worthy of mention is that the 
load line may be modeled as described in Appendix C. 


4.3. Machine description 


The same rpm and pole-pair combinations were used as in 
the synchronous machines. The magnets on the rotor are ar- 
ranged in a cylindrical-wedge configuration, as shown in the 


figure below. 


Figure 7. End view of a permanent magnet machine 


MAGNET PMIATERIAL 


MAGNET STEEL 


BINDER 





The rotor slot factor (lr) used in the other types of 


machines is called here the magnet slot factor and refers to 
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the width of magnet per pole pitch. It may vary between 25% 
and 75% of the pole pitch, the same as for rotor slots in 
other machines. 

The rotor slot depth, ds, does not exist in this 
machine. In this case, the rotor radius, r, is added to ln, 
the magnet radial dimension, to find the actual width of the 
rotating core. (In the synchronous and induction machines, 


dr was included in the rotor radius. ) 


4.3.1. Efficiency 


For the higher rpm machines, efficiency within a par- 
ticular rpm group decreased with an increasing number of 
poles. The most efficient machines, at about 99%, had four 
poles. This is higher than the synchronous or induction 
machines, largely because there are no rotor copper losses. 
The twelve-pole efficiencies were about 98%, which is not 
too large a spread. 

The 180 rpm machines had a fairly flat efficiency curve 
(excepting one anomaly) up to about a 28-pole machine, where 
efficiency started to vary widely. There, the conflict be- 
tween the number of pole-pairs and maximum rotor radius 
started to become significant. The flat efficiency was 
about 96.5%, which is less efficient than the higher rpm 


machines. 


4.3.2. Weight and volume 


The higher rpm machines had a general tendency toward 
lower weight and volume as rpm increased. Within an rpm 
group, the four- and six-pole machines had the lowest weight 
and volume. The smallest machines were larger than the 


synchronous machines. 
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EFFICIENCY 


Figure 8. High rpm permanent magnet machine efficiency 
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EF FICiENCT 


Figure 9. 180 rpm permanent magnet efficiency 
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VOLUME (cubic meters) 


Figure 10. High rpm permanent magnet machine volume 
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VOLUME (M73) 


Figure 11. 180 rpm permanent magnet machine volume 
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Figure 12. High rpm permanent magnet machine weight 


PERMANENT MAGNET MOTOR WEIGHT 


Zo Gore 


(12400 rpm 


WEIGHT (KG) 
(Thousands) 





eae oO Niece Jeo momGet ee es 4 5 6 1 2 S455 6 1 235 45 6 


NUMBER OF POLE—PAIRS 


58 





WEIGHT (KG) 
(Thousands) 
© 


Figure 13. 180 rpm permanent magnet machine weight 
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The weight and volume of the 180 rpm machines increased 
with the number of pole-pairs. Again, the tip speed limita- 
tion is linked to this increase, as with the synchronous 
machines. There was a wide variation around the general in- 
crease, and these machines are large relative to all others. 
They are not competitive as ship propulsion motors because 
of their size. 


4.4. Machine characteristics 


The following tables give the machine characteristics 
of the permanent magnet machines designed for this study. 
The stator slot factor, ls, tended to the maximum of 0.75 
for the higher rpm and 180 rpm machines. This occurred as 
the optimization algorithm tried to minimize envelope volume 
and weight. The rotor slot factor, lr, was a constant 0.378 
for all machines. The rotor slot factor was calculated to 
produce load-line (MEP) operating flux, as derived in Appen- 
dix B. With a different magnet material selection (and a 
consequent change in operating point flux), a different lr 
would have resulted. 

The magnet overhang tended toward the maximum limit, in 
an attempt to achieve the highest flux levels. Permanent 
magnet machines cannot rival the flux level produced by the 
field winding of a synchronous machine, but the optimization 
algorithm did its best. 

The per-unit synchronous reactance-~per-turns-squared 
was limited to 3.0 in these machines. This reactance tended 
toward the limit, but was lower with an increasing number of 
pole-pairs in the 180 rpm machines. It was very difficult 
to achieve valid designs with xsmax = 2.0, as in the 
synchronous machines. If an xsmax greater than 2.0 is unac- 
ceptable, these machines will be less competitive. 

The amount of magnet material varied from 50 kg to a 
few hundred kg in the higher rpm machines to 700-4000 kg in 
the 180 rpm machines. The cost of 700 kg of NdFeB is about 
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Table 11. 


nusber of pole pairs 
power, hp 

efficiency factor 
voluse factor 

sagnet factor 

Shaft rps 

stator current density 
synchronous frequency 
rotor radius 

gap digsension 

back iron depth 

Stator slot depth 
magnet radial disension 
stator slot factor 
rotor slot factor 
envelope voluae 
envelope weight 
hysteresis loss 

eddy current loss 
stator copper loss 
full load efficiency 
active length 

magnet weight 

magnet voluae 
xs/turns-squared, p.u. 
magnet overhang 

air gap flux density 
overall length 


overall diaseter 


] 

29775 
180000 
1284.1 
23 

1800 

7. 49E+06 
30 

0. 2887 
0.0326 
0.1332 
0.0687 
0.0342 
0.750 
0.378 
3.548 
24280.7 
20581.5 
3688. 1 
171150.0 
0.990 
2.328 
405.184 
0.055 
2.994 
0.001 
0.619 
3.750 
1.047 


2 
25775 
180000 
1284.1 
25 
1800 


3 
25775 
180000 
1286. | 
29 


1800 


P 
29775 
180000 
1286. 1 
25 
1800 


") 
25775 
180000 
1286.1 
25 
1800 


1.20€+07 7.32€+06 9. 14E+06 8.00E+04 


40 
0.3396 
0.0382 
0.0870 
0.0742 
0.0401 

0.750 
0.378 
2.908 
12528.8 
20298.5 
7274.8 
265272.5 
0.985 
0.876 
262.576 
0.035 
2.3357 
0.298 
0.687 
2.548 
1.078 


61 


90 
0.3643 
0.0047 
0.0530 
0.0456 
0.0070 

0.750 
0.378 
3.509 
19364.5 
92315.8 
28124.3 
93554.8 
0.991 
2.833 
157.355 
0.021 
2.996 
0.315 
0.642 
4.345 
0.939 


120 
0.5964 
0.0188 
0.0742 
0.0333 
0.0198 

0.750 
0.378 
6.304 
20910.3 
80505.9 
37704.9 
98867.0 
0.988 
0.933 
273.763 
0.037 
1.165 
0.326 
0.722 
3.473 
1.450 


150 
0.5544 
0.0067 
0.0502 
0.0263 
0.0071 

0.790 
0.378 
3.847 
24826. 4 
121495.4 
108857. 1 
69323.3 
0.985 
1.888 
178.613 
0.024 
1.411 
0.331 
0.471 
4.140 
1.275 


Characteristics of 1800 rpm magnet machines 


b 

25775 
180000 
1286.1 
23 

1800 
1.04&+07 
180 
0.6048 
0.0123 
0.0482 
0.0277 
0.0130 
0.646 
0.378 
6.319 
20948. | 
124348. 0 
133717.2 
85761.4 
0.982 
1.257 
246. 460 
0.033 
0.920 
0.329 
0.700 
3.786 
1.390 








Table 12. Characteristics of 2400 rpm magnet machines 


nuaber of pole pairs 
power, hp 

efficiency factor 
voluae factor 

magnet factor 

shatt rpa 

stator current density 
synchronous frequency 
rotor radius 

gap disension 

back iron depth 

Stator slot depth 
@agnet radial dimension 
stator slot factor 
rotor slot factor 
envelope voluae 
envelope weight 
hysteresis loss 

eddy current loss 
Stator copper loss 
full load efficiency 
active length 

magnet weight 

gagnet voluae 
xs/turns-Squared, p.u. 
magnet overhang 

air gap flux density 
overall length 


overall! diaaeter 


l 

29779 
180000 
1286.1 
25 

2400 
12000000 
40 
0.3594 
0.0371 
0.1640 
0.0697 
0.0390 
0.348 
0.378 
4,497 
26948. 1 
54028.8 
8130.4 
234324. 1 
0.986 
1.956 
384.207 
0.052 
2.417 
0,018 
0.617 
3.278 
1.260 


y 
y 


29779 
180000 
1286.1 

29 

2400 
12000000 
80 
0.00035 
0.1883 
0.0420 
0.0249 
0.1978 

0.790 

0.378 

1.342 
8797.0 

18160.1 
8677.9 
217110.2 

0.987 

4.392 
10.204 

0.001 

0.282 

0.330 

0.639 

9.938 

Q.o11 


3 

29779 
180000 
1286. 1 
2) 

2400 
12000000 
120 
0.4018 
0.0136 
0.0625 
0.9341 
0.0143 
9.730 
0.378 
HP) 
12333.3 
46136.0 
33069, 4 
127316.4 
0.989 
1.291 
167.444 
0.023 
1.973 
0.331 
0.676 
3.910 
1.024 


4 
29779 
180000 
1286.1 
25 

2400 
10202760 
160 
0.4408 
0.0088 
0.0504 
6.0325 
0.0092 
0.730 
0.378 
5:19) 
Sones 
68103.9 
65087.9 
90617.2 
0.988 
1.387 
128.106 
0.017 
1.965 
0.319 
0.672 
De2Z 
1.065 


& 
J 


23/78 
180000 
1286.1 

25 

2400 
11959030 
200 
0.3361 
0.0032 
0.0290 
0.0234 
0.0033 

0.749 

0.378 

2.349 

13660. 4 
86724.9 
103604. 3 
112885. 8 

0.984 

3.061 
13.723 

0.016 

2.879 

0.3538 

0.640 

4.431 

0.783 


6 

falda 
180000 
1286.1 
ZS 

2400 
11997423 
240 

0. 4760 
0.0169 
0.0379 
0.0300 
0.0173 
0.678 
0.378 
3.438 
11877.1 
91478.9 
131141.0 
91721.4 
0.984 
1.129 
217.676 
0.029 
0.751 
0.324 
0.692 
3.168 
1.121 








Table 13. 


nuaber of pole pairs 
power, hp 

efficiency factor 
voluae factor 

magnet factor 

Shaft rpa 

stator current density 
synchronous frequency 
rotor radius 

gap dimension 

back iron depth 

stator slot depth 
magnet radial dimension 
stator slot factor 
rotor slot factor 
envelope voluae 
envelope weight 
hysteresis loss 

eddy current loss 
stator copper loss 
full load efficiency 
active length 

magnet weight 

magnet voluae 
xs/turns-squared, p.u. 
magnet overhang 

air gap flux density 
overall length 


overall diaseter 


] 

20775 
180000 
12864. 1 
29 

5000 
10639824 
90 
0.2480 
0.0227 
0.1121 
0.0405 
0.0238 
0.730 
0.378 
(ary 
15091.7 
22692.8 
6777.4 
197512.6 
0. 990 
2.3595 
290.192 
0.034 
2.934 
0.007 
0.618 
3.974 
0.847 


2 

TVR) 
180000 
1286.1 
29 

5000 
9957320. 
100 
0.2394 
0.0075 
0.0510 
0.0369 
0.0079 
0.730 
0.378 
1.834 
12594./ 
56848. 2 
22010.1 
108787.2 
0.991 
5.719 
123.693 
0.017 
2.995 
0.009 
0.619 
4.738 
0.669 


K 

29/79 
180000 
1286.1 
29 

500 
12000000 
130 
0.2239 
0.0035 
0.0317 
0.0235 
0.0037 
0.790 
0.378 
1.771 
12381. 0 
07339,4 
v1374.8 
132086. 7 
0. 988 
32490 
71.981 
0.012 
2.997 
0.526 
0.627 
6.415 


0.560 


4 

29779 
180000 
1286.1 
7 

3000 
11999567 
200 
0.0003 
0.1603 
0.0175 
0.0074 
0.1683 
0.790 
0.378 
2.143 
8204.8 
45284.3 
94098. 3 
170766.7 
0.986 
16.723 
13.689 
0.002 
0.043 
0.325 
0.623 
18.039 
0.371 


s) 

20/79 
180000 
12864. 1 
29 

5000 
12000000 
290 
0.3116 
0.0027 
0.0267 
0.0224 
0.0028 
0,747 
0.378 
1.990 
11507.2 
91118. 4 
136066./ 
97680.0 
0.983 
3.005 
0. 769 
0.008 
feet 
0.311 
0.637 
4.274 
0,727 


Characteristics of 3000 rpm magnet machines 


b 

20773 
180000 
1286.1 
29 

$000 
11400662 
300 
0.4837 
0.0375 
0.0428 
0.0312 
0.0394 
0.681 
0.3/8 
3.643 
7968.3 
99620.6 
171347.8 
74849.0 
0.983 
0.732 
338.030 
0.044 
0.351 
0.336 
0.736 
Za ia 
1.191 





Table 14. 


nuaber of pole pairs 
power, hp 

efficiency factor 
voluae factor 

magnet factor 

Shaft rpa 

stator current density 
synchronous frequency 
rotor radius 

gap dimension 

back iron depth 

Stator slot depth 
magnet radial disension 
Stator slot factor 
rotor slot factor 
envelope voluae 
envelope weight 
hysteresis loss 

eddy current loss 
Stator copper loss 
full toad efficiency 
active length 

Magnet weight 

magnet voluae 
xs/turns-Squared, p.u. 
@agnet overhang 

alr gap flux density 
overall length 


overall diaaeter 


] 

29779 
180060 
1286. 1 
20 

3600 
12000000 
&0 
0.3044 
0.0228 
0.1494 
0.0332 
0.0240 
0.790 
0.378 
2.475 
19448.9 
29048.1 
10410.4 
163951.5 
0.990 
1.594 
Pape Ue 
0.030 
2.960 
0.336 
0.664 
2.799 
1.012 


Z 

23779 
180000 
1286.1 
20 

3600 
12000000 
120 
0.2596 
0.0088 
0.0574 
0.0442 
0.0092 
0.003 
0.3/8 
1.602 
9853.6 
36290.0 
26012.0 
114166.3 
0.991 
2.276 
113.409 
0.019 
2.977 
0.324 
0.641 
3.386 
0.740 


3 

29775 
180000 
1286.1 
29 

3400 
11966619 
180 
0.35006 
0.0121 
0.0462 
0.0411 
0.0127 
0.750 
0.378 
1.427 
7072.7 
37383.1 
40193.3 
104037.5 
0.991 
1.279 
104.478 
0.014 
2.06] 
0.327 
0.665 
2.081 
0.800 
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4 
29779 
180000 
1286.1 
29 
3600 


5 
Wd 


29779 
180060 
1286.1 

20 
3600 


14677983 11998747 


240 
0.2526 
0.0025 
0.0268 
0.0208 
0.0026 

0.735 
0.378 
1.694 
11238.1 
94962.0 
121798.6 
99474.6 
0.985 
4.322 
98.202 
0.008 
2.981 
0.333 
0.630 
deage 
0.605 


300 
0.3700 
0.0063 
0.0332 
0.0223 
0.0067 

0.730 
0.378 
2.043 
9644.1 
92749.9 


Characteristics of 3600 rpm magnet machines 


b 

20779 
180000 
1286.1 
2 

3600 
14975558 
340 
0.4934 
0.0173 
0.0371 
0.0240 
0.0182 
0.749 
0.378 
2.801 
8977.4 
104238.5 


166203.6 224148.8 


73899.8 
0.983 
1.469 

M.42k 
0.012 
1,490 
0.340 
0.461 
3.202 
0.864 


66008. | 
0.980 
0.907 

175.220 
0.024 
0.993 
0.339 
0.708 
2.83 
1.064 





Table 15. Characteristics of 7200 rpm magnet machines 


nuaber of pole pairs 1 2 3 4 y, ) 
power, hp 20775 = 2377525775 2877928775 25775 

efficiency factor 180000 180000 180000 180000 180000 180000 
voluae factor 1286.1 1286.1 1286.1 1286.1 1286.1 1284.1 

magnet factor A 29 25 29 25 25 

shaft rpe 7200 7200 7200 7200 7200 7200 

Stator current density 12000000 12000000 8393050. 12000000 12000000 12000000 
synchronous frequency 120 240 340 480 600 720 
rotor radius 0.2266 0.1938 0.2622 0.2627 0.2467 0.2418 

gap di@ension 0.0146 0.0089 0.0095 0.0041 0.0020 0.0022 

back iron depth 0.1045 06.0423 0.0392 0.0289 0.0209 0.0171 

stator slot depth 0.0285 0.0284 0.0315 0.0244 0.0114 0.0104 
magnet radial dimension 0.0153 0.0062 06.0099 0.0043 0.0021 0.0023 
Stator slot factor 0.749 0.750 0.750 0.7350 0.748 0.750 
rotor slot factor 0.378 0.378 0.378 0.378 0.378 0.378 
envelope volume 1.247 0.811 1413 0.931 ieogo | leo 
envelope weight 8186.9 5319.9 6052.8 4898.9 8797.5 9221.7 
hysteresis loss $0513.8 3$8176.1 65438.7 72593.1 172996.0 218032.2 

eddy current loss 21B871.7 94727.9 140716.0 208133.7 620002.2 937689.5 





stator copper loss 90777.3 71500.7 38515.8 32482.8 49628.2 49815.1 
full load efficiency 0.993 0.992 0.987 0.983 0.958 0.941 
active length 1.952 2.384 1.621 1.5942 3.963 4.508 

magnet weight 110.276 57.403 87.115 36.486 = 42.995 50.295 

aagnet volume 0.015 0.008 0.012 0.005 0.006 0.007 
xs/turns-squared, p.u. 2.946 2.884 1.258 2.246 1.608 1.142 
magnet overhang Vo sso leeon0n507)  U.cze 0.500 0.519 

air gap flux density 0.648 0.635 0.647 9.690 0.631 0.629 
overall length 2.578 3.207 2.747 2.627 4.968 9.493 

overall) diameter 0.748 0.941 0.4685 0.640 0.562 0.545 





Table 16. Characteristics of 180 rpm magnet machines 


nueber of pole pairs 
power, hp 

efficiency factor 
voluae factor 

@agnet factor 

Shaft rpe 

Stator current density 
synchronous frequency 
rotor radius 

gap digension 

back iron depth 

Stator slot depth 
@agnet radial digension 
stator slot factor 
rotor slot factor 
envelope voluae 
envelope weight 
hysteresis loss 

eddy current loss 
Stator copper loss 
full load efficiency 
active length 

sagnet weight 

magnet voluae 
xs/turns-squared, p.u. 
magnet overhang 

alr gap flux density 
overall length 


overall diaaeter 


| 
29775 
180000 
1286.1 
23 

180 

7, BOE+06 
5 
0.4399 
0.0413 
0.1994 
0.0932 
0.0433 
0.463 
0.578 
29.375 


l. 


2 
23/79 
180000 
1286. 1 
29 

180 
1BE+07 
6 
0.1787 
0.0036 
0.0377 
0.0282 
0.0038 
soul 
0.378 
36.694 


211948.9 284972. 1 
20311.6 53970.7 


564.0 


1934.5 


7.90E+05 2. 04E+96 


0.940 
12.100 
4072.4 

0.950 

ribw piss 

0.003 

0.619 
14.198 

1.548 


6&4 


0.881 
171.774 
2079.6 
0.281 
2.990 
0.034 
0.619 
172.918 
0.4% 


5 

20779 
180000 
1286. ! 
23 

180 

1. 20E+07 
9 

0.9745 
0.0201 
0.1491 
0.0327 
0.0212 
0.730 
0.578 
34.709 
106312.7 
47089.7 
2551.0 
6. 97E+09 
0.963 
Salvia 
1941.0 
0.262 
2.981 
0.532 
0.674 
7.238 
2.595 


+ 

29/79 
180000 
1286.1 
29 

180 

1. 09E+07 
12 
0.7389 
0.0166 
0.0779 
0.0476 
0.0174 
0.750 
0.378 
20.521 
107071.7 
41032.7 
2941.1 
7.Q1E+05 
0.963 
4,960 
1048.8 
0.142 
2.983 
0.010 
0.618 
7.651 
1.762 


n) 

2ovi9 
180000 
1284.1 
25 

180 
1.20E+07 
15 
0.8993 
0.0214 
0.0865 
0.0630 
0.0225 
0.740 
0.378 
22.314 
77934.3 
36926.0 
3308.5 
7. 11E+05 
0.962 
1,866 
1087.7 
0.147 
2.870 
0.3539 
0.704 
ages 
2.140 





Table 17. Characteristics of more 180 rpm magnet machines 


nuaber of pole pairs 
power, hp 

efficiency factor 
voluae factor 

@agnet factor 

shaft rpe 

Stator current density 
synchronous frequency 
rotor radius 

Gap digension 

back iron depth 

Stator slot depth 
Magnet radial dieension 
Stator slot factor 
rotor slot factor 
envelope voluee 
envelope weight 
hysteresis loss 

eddy current loss 
Stator copper loss 
full load efficiency 
active length 

@agnet weight 

magnet voluae 
xs/turns-squared, p.u. 
@agnet overhang 

air gap flux density 
overall length 


overall diameter 


b 

29775 
180000 
1284.1 
23 

180 

6. 03E+06 
18 
0.6047 
0.0139 
0.0426 
0.0337 
0.0146 
0.036 
0,378 
445/78 
313139.4 
186214.5 
20021. 3 
3. F1E+05 
0.970 
24,182 
3800. 1 
0.314 
0.046 
0.010 
0.619 
26.715 
1.390 


&7 


? 8 9 10 

olds = 207750 25775 2577S 
180000 180000 180000 180000 
1286.1 1286.1 1284.1 1286.1 
29 29 23 23 

180 180 180 180 

1. 20E+07 1.20E+07 9.42E+06 1.20E+07 
21 24 ra 30 
1.2439 0.9696 0.9792 0.9526 
0.1161 6.0089 0.0031 0.0079 
0.1117 0.0546 0.0469 0.0430 
0.0627 0.0312 0.0255 0.0340 
0.1219 0.0093 0.0033 0.0083 
0.750 0.748 0.575 0.748 
0.578 0.378 0.378 0.378 
We024 = 29,253 43.668 = 26.292 
78601.8 120336.5 234296.3 102286.2 
67918.3 96060.7 214344.0 101073.4 
8469.3 13770.9 34568.5 18111.9 
7. ASE+0S 3.34E+05 3. 91E+09 9. 51E+05 
0.960 0.968 0.968 0.967 
0.651 3.523 7.417 3$.16/ 
3431.9 932.5 704.3 714.1 
0.464 0.126 0.095 0.097 
Gece UE ee. 7 10 ne 2250) 
0.329 0.336 0.340 0.325 
0.859 0.669 0.644 9.471 
6.579 7.474 11.360 7.042 

S Uoveemecetso ge 2.107 9) 12.077 





$150,000, rendering the 180 rpm machines less economic to 
build. A machine such as these may cost between four and 
eight million dollars. Only those machines with magnet 
costs that are reasonable with respect to the other material 
costs should be candidates for design. 

In all the permanent magnet machines, stator current 
density went to the maximum. The overall length and overall 
diameter of the more reasonable machines would allow them to 


fit in machinery spaces aboard a ship. 


4.5. Verification 


No high-power permanent magnet machines were discovered 
during the search to find a benchmark. Because of the high 
material cost and the competition afforded by synchronous 
and induction machines, it seems none have been built. 
several paper studies were found [13, 19, 20, 21, and 22}, 
and the parameters resulting from this computer modeling 
seem to agree with them. The machine size is what was ex- 
pected, given the lower air gap flux density. The ef- 
ficiency was higher than the synchronous and induction 


machines. All-in-all, this modeling gave good machines. 
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Chapter Five. Induction machines 


The stator of an induction machine is the same as those 
of synchronous and permanent magnet machines. The rotor is 
Significantly different. There is no independent mechanism 
to produce a rotor flux wave. The rotor winding is shorted, 
whether it is wound or cast, so that as the stator flux wave 
passes over the rotor, currents are induced in the winding. 
These currents produce only a small reaction flux, but it 
still tends to align with the stator flux wave. When at 
operating speeds, the rotor speed is a bit slower than the 
stator flux wave speed, and the difference in speeds is 
called slip. Typically, slip is a few percent of the stator 
frequency. The rotor currents are at slip frequency. If 
the rotor and stator speeds were the same, slip would be 
zero, there would be no tendency to align and torque would 
be zero. Then, the rotor would lag behind the stator until 
current was induced in the rotor winding by the passing 
stator flux wave and torque was again produced. 

If solid bars are used as the rotor winding, they are 
shorted at the ends of the rotor by end rings, to form what 
is called a “squirrel cage" rotor. If actual turns are 
used, the winding may be shorted through external resis- 
tances to affect the starting and torque-slip characteris- 
tics of the machine. Fitzgerald et al [23] and Alger [24] 


discuss induction machine characteristics in some detail. 


5.1. Assumptions 


A squirrel cage rotor was assumed for these machines. 
Copper was designated as the material for the rotor bars. 
However, these machines will be fed from a frequency 
changer, so only one layer of bars was used and the effects 
of magnetic diffusion ignored in the analysis (see 
Appendix D for a derivation of the components of an induc- 


tion machine equivalent circuit). The number of rotor bars 


69 





was arbitrarily set at 71. This quantity should not cause 
undesirable harmonics, as it will not be an integral mul- 
tiple of the number of poles or stator slots in any machine. 
The number and width of the rotor bars were inextricably 


entwined and could not be separated in the analysis. 


5.2. Machine description 


An attempt was made to design the same rpm and pole- 
pair machines as was performed for the synchronous and per- 
manent magnet models, but problems in limiting rotor current 
density allowed only a few of the machines to be designed. 
For example, no 7200 rpm machines were designed and 180 rpm 
machines could only be designed with up to twelve poles. 

All of the induction machines are listed in the tables 
starting on page . Only medium confidence should be placed 
in the induction machine designs, as there were some conver- 
gence difficulties in slip. (Slip is not listed in the 


tables for that reason. ) 
5.2.1. Efficiency 


The higher rpm machines showed a slight increase in ef- 
ficiency as rpm increased. There was much movement around 
the average value of 97.5%. The movement decreased as rpm 
increased. With only six machines, it is hard to detect a 
trend in 180 rpm machine efficiency. Apparently, efficiency 
did increase with the number of pole-pairs, with all ef- 
ficiencies below 90%. Developmental studies for this thesis 
showed that off-design-point efficiencies for the 180 rpm 
machines were sometimes below 70% for the endurance speed 


condition. 
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Figure 14. High rpm induction machine efficiency 
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Figure 15. 180 rpm induction machine efficiency 
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Figure 16. High rpm induction machine volume 
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Figure 17. 


180 rpm induction machine volume 
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Figure 18. High rpm induction machine weight 
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Figure 19. 180 rpm induction machine weight 
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5.2.2. Weight and volume 


The best of the high rpm machines rivaled the 
synchronous machines in weight and volume. The worst were 
very bad. Weight and volume decreased with an increase in 
the number of poles, but not necessarily with the increases 
in rpm. For the 180 rpm machines, both weight and volume 
decreased dramatically as the number of poles went from two 
to six, with much lower decreases after that. The 180 rpm 


machines were uncompetitive in the synthesis process. 


5.3. Machine characteristics 


The previously mentioned rotor current density dif- 
ficulty showed in the rotor slot factor, which was at the 
limit of 0.75 for almost every motor. The stator slot fac- 
tor gradually grew with the increase in poles, arriving at 
0.75. The length and diameter of both the 180 and higher 
rpm machines is such that they would fit in machinery 


spaces. 


5.4. Verification 


Induction machines were expected to be close to 
synchronous machines in volume, weight, and efficiency. 
They were, and this comparison served as the verification 
for the induction machine model. Because the confidence 
level in the designed machines is only medium, more work 
would be needed to verify that these machines would have the 
advertised properties if built. 
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Table 18. Characteristics of 1800 rpm induction machines 


nuaber of pole pairs 
power, hp 
efficiency factor 
VOluge factor 

shaft rpa 

Oriaary aep-turns 
synchronous frequency 
rotor radius 

gap diaension 

back iron depth 
stator slot depth 
rotor slot depth 
Stator slot factor 
rotor slot factor 
envelope volume 
envelope weight 
hysteresis loss 
eddy current joss 
Stator copper loss 
full load efficiency 
active length 

rotor copper loss 
maxigua torque 
terainal volts/turn 
air gap volts/turn 
R1/turns-squared 
Xi/turns-squared 
Xa/turns-Squared 
X2/turns~squared 
R2/turns-squared 
overall length 


overall diaaeter 


| 

29775 
180000 
1286.1 
1800 

6. L0E+05 
30 
0.9236 
0.0262 
0.1894 
2.1467 
0.205! 
0.285 
0.750 
148.843 
376603.7 
1B0991.2 
32425.7 
1. 03E+06 
0.939 
0.191 
8427.] 
185470 
234.31 
63.29 


2 
29779 
1B0000 
1286.1 
1800 


3 
AVES 
180000 
1286.1 
1800 


7, 74E+05 2, 44E+05 


60 
1.0197 
0.0039 
0.0173 
0.9423 
§.3399 

0.641 
0.790 
46.240 
B0669.2 
Ze 
6653.2 


90 
0.4638 
0.0036 
0.1073 
0.5379 
0.1546 

0.399 
0.750 
7,508 
29447.0 
68838. 7 
37006.8 


1.45E+06 498648. 3 


0.928 
0.090 
8943.0 
185470 
45.07 
32.92 


0.969 
0.799 
19627.2 
185470 
266.49 
126.58 


4 
20770 
180000 
1286.1 
1B00 
1.91£+05 
120 
0.4392 
0.0029 
0.0093 
0. 3490 
0.145! 
0.640 
0.750 
3.880 
21499.7 
39013.4 
27964. 1 


c 
to} 


29779 
180000 
1284.1 

1800 

1, B6E+05 
130 
0.4040 
0.0049 
0.0405 
0.1902 
0.1353 

0.699 

0.750 

3.760 

13333. 1 
42456. 7 
38040. 3 


6 

29779 
180000 
1284. 1 
1800 
1.79E+05 
180 
0.35945 
0.0051 
0.0440 
0.1388 
0.1176 
0.749 
0.748 
3.164 
erost 
48214.8 
91839.3 


165391.9 264414.2 320360.8 


0.987 
0.930 
22347.3 
185470 
233.96 
145.45 


0.981 
1.000 
30471.4 
185470 
193.76 
145.98 


0.977 
1.081 
37148.9 
185470 
189.53 
193.34 


9.23E-07 B.09E-07 2.80E-06 1.90E-06 2.54E-06 3. 32E-06 
3.70E-04 3.89E-05 9.60E-04 9.53E-04 6.B2E-04 6. 17-04 
2,91E-03 4.41E-03 1.23E-02 1. 286-02 6.23E-03 9. 15&-03 
2.03E-05 3.46£-05 3.76E-04 6,.52E-04 7.79E-04 9. 08E-04 
4,50E-08 9.97E-09 3.89£-07 3.37£-07 6.62E-07 B.49E-07 

3.990 4.184 2.62% 2.682 2.644 2.680 

6.972 3.166 1.824 1.593 1.283 1.169 
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Table 19. Characteristics of 2400 rpm induction machines 


nuaber of pole pairs 
power, hp 

efficiency factor 
voluae factor 
shaft rpa 
prigary aap-turns 
synchronous frequency 
rotor radius 

gap dimension 

back iron depth 
stator slot depth 
rotor slot depth 
stator slot factor 
rotor slot factor 
envelope voluae 
envelope weight 
hysteresis loss 

eddy current loss 
stator copper loss 
full load efficiency 
active length 

rotor copper loss 
waxigua torque 
terainal volts/turn 
air gap volts/turn 
Ri/turns-squared 
X1/turns-squared 
Xa/turns-squared 
X2/turns-squared 
R2/turns-squared 
overal] length 


overall diaseter 


2 3 4 x) 6 

29773) 0 239773) 23775 2877520778 
180000 180000 180000 180000 180000 
1286.1 1286.1 1286.1 1286.1 1286.1 
2400 2400 2400 2400 2400 

3. 94E+09 3. 7BE+09 2.456405 1. 7BE+09 1.93E+09 
80 120 160 200 240 
0.7541 9.7998 0.4803 0.4126 0.3841 
0.0069 0.0020 0.0094 06.0087 06.0081 
6.0992 0.0013 0.0834 0.0375 0.0448 
0.6169 0.5475 0.2755 0.2341 0.156! 
0.2514 0.1979 0.1601 0.1375 0.1197 
0.454 0.477 0.417 6.413 0.730 
0.730 0.790 0.639 0.568 6.750 
23.889 21.180 6.201 4.554 2.942 
60526.1 30880.0 19931.3 17243.1 11499.0 
78227.5 26411.4 80426.4 86278.8 %54014.3 
37381.4 18716.2 76864.2 103071.8 77433.6 
1. 10E+06 334484.0 376335.4 306069.2 282967.9 
0.940 0.980 0.972 0.974 6.978 
0.142 0.189 0.933 0.909 0.832 
9819.4 7426.0 17138.3 21989.8 27994.3 
139102.5 139102.9 139102.5 139102.9 139102.5 
103.23 113.69 229.60 346.92 212.939 
glee col A ele eo? 179.70 16.78 
1.046-06 7.82E-07 2.08E-06 3.23E-06 2.92E-06 
1,90E-04 2.32E-04 7.89E-04 1.67E-03 7.40E-04 
3.B8E-03 1.25E-02 3.37E-03 4.31E-03 3.35E-05 
6.22E-05 1.56E-04 4.63£-04 1.07E-03 9. 36E-04 
2.84E-08 4.44€-08 2.97E-07 7.69€-07 &.73E-07 
3.186 863.380) = 2.492) 2.594 =... 420 
2.946 2.693 1.697 1.426 1.186 
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Table 20. Characteristics of 3000 rpm induction machines 


nuaber of pole pairs 2 3 4 a, b 
power, hp ipioeeesiia “a2ayla \2os7a. | 25/75 
efficiency factor 180000 180000 180000 180000 180000 
volume factor 1286.1 1286.1 1286.1 1286.1 1286.1 

shaft rps 3000 3000 3000 3000 3000 

prigary aap-turns S.10E+05 3.62E+05 2.80E+05 1.76E+05 1.37E+05 
synchronous frequency 100 150 200 250 300 
rotor radius 0.9198 9.6366 0.5580 0.3660 0.3325 

gap disension 0.0138 0.0036 0.0070 0.0070 0.0039 

back iron depth 0.1799 0.0134 0.0178 0.0503 0.0388 
stator slot depth 0.8926 0.3785 0.3668 0.2274 0.1288 
rotor slot depth Oviusee Oneree 0.1960 0.1220 0.1093 
stator slot factor 0.340 0.607 0.545 0.476 0.721 
rotor slot factor 6.750 0.750 0.750 0.750 0.750 
envelope volume 22.412 10.122 98.000 3.379 2.065 
envelope weight 72906.0 19541.1 15961.5 13000.9 9016.0 
lysteresis loss 185978.6 25156.1 44235.6 74610.2 32848.0 
eddy current loss 111088.5 22539.3 92845.5 111415.1 94701.2 
Stator copper loss 322718.0 $59218.8 221369.1 269959.2 226785.8 
full load efficiency 0.968 0.979 0.983 0.976 0.980 
active length 0.383 0.189 0.307 0.818 1.006 

rotor copper loss 8574.1 9089.1 9006.9 18292.9 25548. 1 
Maximua torque W11Z82 0 11282 = 41282 = 111282 =—:111282 
terminal volts/turn 421.67 101.44 174.86 361.45 251.07 
air gap volts/turn 119.20 72.02 102.80 179.30 200.92 
Ri/turns-squared 1,12E-06 9. 146-07 9.39€-07 2.91E-06 4.01E-06 
X1/turns-squared 1. 316-03 1.96€-04 5.04€-04 1.78E-03 1.09E-03 
Xa/turns-squared AAAJE@US 6. 95E-05 3,.82E-03 3.33E-05 9.00E-03 
X2/turns-squared 1.82E-04 1.72E-04 3.25€-04 9.97E-04 1.54€-03 
R2/turns-squared 1.69E-07 5.15E-08 1.08E-07 6.66€-07 1.01E-06 
overall length Pilioeeetg0 92.90) 2.510 2,502 
overall diameter 3.210 2.064 1.899 1.301 1.008 
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Table 21. 


nuaber of pole pairs 
power, hp 
efficiency factor 
voluae factor 

Shatt rpa 

prigary aap-turns 
synchronous frequency 
rotor radius 

gap digsension 

back iron depth 
stator slot dept’ 
rotor slot depth 
stator slot factor 
rotor slot factor 
envelope voluae 
envelope weight 
hysteresis loss 
eddy current loss 
stator copper loss 
full load efficiency 
active length 

rotor copper ljass 
@axiaua torque 
terainal volts/turn 
alr gap volts/turn 
R1/turns-squared 
Xi/turns-squared 
Xa/turns-squared 
X2/turns-squared 
R2/turns-squared 
overall length 


overall diaaeter 


i 

29779 
180000 
12864. 1 
3600 
2.72E+09 
&0 
0.4964 
0.0061 
0.3428 
1.3891 
0.1694 
0.290 
0.730 
41,461 
164441.9 
261300.3 
93647.7 
541600. 1 
0.965 


140.25 
 J4E-06 
1.47E-03 
2.S9E-02 
1. 20E-04 
1.98E-07 

2.403 

4,469 


—- 


2 

20779 
180000 
1284. 1 
3400 
1.70E+05 
120 
0.4350 
v.0027 
0.0133 
1.0052 
0.1308 
0.338 
0.790 
17.940 
94127.3 
145818.0 
104519.7 


3 
een 
180000 
1286.1 
3600 
2.99E+09 
180 
0.9100 
0.0033 
0.0238 
0.3606 
0.1717 
0.585 
0.790 
§.05 
14922. 3 
S4381.6 
36966, 2 


107764.3 288857.8 210564.4 


0.981 
0.700 
122023 
92739 
702.41 
218.89 


0.981 
0,247 
10787.1 
92730 
143.82 
91.39 


1,25E-06 1.07E-06 8.08E-07 1.37E-06 
S.93E-03 3. 70E-04 S.O7E-04 1, 13E-05 
4, A9E-02 9,66E-03 2.80E-05 3. 94E-03 
4, 65E-04 2.54E-04 3.05E-04 7, 24E-04 
4, 97E-07 1.03E-07 9,48E-08 2.79£-07 


2.431 
2.913 


81 


2.320 
1.805 


4 =) 

LE MENT 
180000 180000 
1284.1 1206.1 
3460 3400 
2.99E+09 2.05E+05 
240 500 
0.5309 0.4338 
0.0088 0.0079 
0.0238 0.0258 
0.3236 0.2792 
0.1768 0.1446 
0.688 0,993 
0.790 0.750 
6.936 4.331 
14074.1 11549.4 
42814.0 62801.7 
61376.6 112537.8 
173450. 1 

0.983 0.982 
0.249 0,480 
9994.8 11162.4 
727353 72735 
142.32 279.41 
94.86 149.83 
2.406 2.247 
1.773 1.494 


Characteristics of 3600 rpm induction machines 


) 

29779 
180000 
1286. | 
3400 

1. 25E+05 
340 
0.3144 
0.0030 
0.0367 
0.1181 
9.1049 
0.750 
0.749 
Weize 
7993.8 
02473, 1 
112835.4 
192320.6 
0.981 
0,964 
23436. 
92739 
261.98 
218.01 
4, 256-06 
1. 18E-03 
1. 26E-02 
1. 84E-03 
1.07E-06 
2.234 
0.944 





Table 22. 


nuaber of pole pairs 
power, hp 
efficiency factor 
voluae factor 

shaft rpa 

primary aaps-turas 
synchronous frequency 
rotor radius 

gap diaension 

back iron depth 
stator stot depth 
rotor slot depth 
stator slot factor 
rotor slot factor 
envelope volume 
envelope weight 
hysteresis loss 
eddy current joss 
stator copper joss 
full load efficiency 
active length 

rotor copper joss 
@axiaua torque 
tergainal volts/turn 
alr gap volts/turn 
Ri/turns-squared 
Xi/turns-squared 
Xa/turns-squared 
X2/turns-squared 
R2/turns-squared 
overall Length 


averali diaseter 


i 

29/79 
180000 
1286.1 
180 
1034077 
3 
0.9954 
0.0025 
0.0227 
0.6714 
0.3519 
0.590 
0.790 
47.668 
217524. 1 
4551.5 
81.4 

4, 44E+06 
0.808 
0.829 
1.086 +05 
1854700 
51.849 
29.510 
1. 336-06 
2. 08E-05 
1. 21£-02 
1.89E-05 
1.16&-07 
4.817 
3.584 
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fe 
25779 
180000 
1286.1 
180 


4 
29779 
180000 
1284, 1 
180 


264971.9 94697335.6 


6 
0.7234 
0.0060 
0.0240 
0.7116 
0.2412 

0.0635 
0.790 
36.837 


q 
0.6896 
0.0199 
0.1968 

2006 
0.2299 
0.734 
0.790 
20.316 


168344.2 98690.0 


161944 
080.4 


19336.0 
1940.4 


4 
25779 
180000 
1286. 1 

180 


c 
J 


29779 
180000 
1286.1 

180 


710008.9 912934.4 


We 
0.6614 
0.0334 
0.1198 
0.2231 
0.2039 

0.750 
0.683 
19.490 
94954.9 
24244.4 
1741.4 


15 
0.7411 
0.0239 
0.1038 
0.2200 
0.1845 

0.594 
0.696 
21.3597 
95084.9 
351609, 4 
2832.1 


1. 96E+06 S.05E+06 S.15E+06 2.59E+06 


0.901 
2.04] 


0.894 
1.854 


0.850 
2.486 


0.872 
2.159 


Characteristics of 180 rpm induction machines 


) 

25779 
180000 
1284.1 
180 
524264.9 
ig 
0.6962 
0.0054 
0.0812 
0.1464 
0.1901 
0.718 
0.701 
17.600 
B7958.9 
36548.4 
ores <8 
2. 09E+06 
0.891 
5.090 


1. 29E+05 2.0BE+05 2.23E+05 1.94E+05 2.14E+05 


1854700 
104.401 
93.243 


1894700 
08.731 
43.967] 


1854700 
76.930 
09.124 


1854700 
76.097 
07.910 


1854700 
§7.418 
77.34 


2.05E-06 3. 146-06 4.04E-06 3.2BE-06 6. 51E-06 
9,01£-05 9.356-05 1.12E-04 


1, 586-04 6.01E-05 


4, S1£-03 §.17E-04 4.60E-04 5. 08E-04 2.49E-03 
1. 37£-04 1.27E-04 2. 40E-04 


Bit lesvon a F1E-03 


4.44E-07 4. 336-07 7.40E-07 &.20E-07 1.15E-04 


4,969 
2.930 


4,91 
2.257 


2.266 
2.068 


di2l/ 
2.177 


9.897 
1.859 
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Chapter Six. Nominal ship design 


6.1. Technology sensitivity analysis 


Technology sensitivity analyses, such as this thesis, 
must be able to quantitatively compare similar technologies. 
The pertinent differences must be made apparent through ap- 
propriate analysis. Inherent in the analysis must be the 
consideration of the global system complexity. Naval ships 
are extremely complex and the effects of various tech- 
nologies can be lost in the complexity. One methodology for 
technology characterization in naval ships has been proposed 
by Goddard [25]. This method has been followed to show the 


benefits of electric drive. 


6.2. ASSET 


The Advanced Surface Ship Evaluation Tool (ASSET) was 
developed over several years to be the U. S. Navy’s premier 
ship design computer program. It has its roots in HANDE, a 
hydrofoil design program developed in the 1970’s. The ship 
design spiral is traversed in an iterative fashion until 
convergence on a number of parameters is achieved. Boeing 
Computer Services is the contractor for ASSET, under the su- 
pervision of the David Taylor Naval Ship Research and 
Development Center, at Carderock, Maryland.!6& [It is divided 
into large modules by ship type. These modules include 
Monohulls, Hydrofoils, and SWATH ships. 

The geometry of a particular ship is input to ASSET. 
The following were used as characteristic ship traits: full 
load displacement, certain Ship Work Breakdown Structure 
(SWBS) weight groups,17 endurance fuel load at 20 knots, 


16. Greenwood, R.W., and Fuller, A.L., “Development of a 
Common Tool for Ship Design and Technology Evaluation, ” 
Proc. SNAME New England Section Marine Computers 1986. 


17. See Appendix E for more information on SWBS. 
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draft, maximum and sustained speed for 51,550 installed 
horsepower, and transverse GM. From the changes in these 
parameters during the various computer runs, the effects of 
transmission choices were noted. 

The inputs to ASSET describe the ship that is being 
designed. The outputs cover the range of calculations pos- 
sible in structures, volume, space, machinery, propeller 
characteristics, resistance, powering, and weight. There 
are ASSET performance modules on cost, stability, hydros- 
tatics, seakeeping, manning and space but the usual syn- 
thesis output is of more use during a technology sensitivity 
analysis. 

The descriptions of several ships are contained in an 
ASSET data base. For a particular ship, a Current Model is 
maintained that holds all of the parameters to describe that 
particular ship. In ASSET Version 2.0, over 380 parameters 
are used for each ship description. User control over most 
of these parameters is possible, or control may be given to 
the executive program which will then "design" a ship sub- 
ject to whatever constraints the user desires. 

Some intricacies of ship design are not handled well by 
ASSET. For example, the program is not able to handle 
equipment re-arrangements easily, and almost all equipment- 
level volumes are approximated from studies of past ships.18 
For this reason, some equipment-level weights and volumes 
need to be calculated off-line and input to the program 
through its weight adjustment facility, especially if ac- 
curacy in these areas is important to the study being per- 
formed. 

The baseline ship used in these studies has a full load 
displacement of 5485 LT, carries 272 crew members, is 425 


feet long, and has a primary mission of anti-submarine war- 


18. The Enhanced Machinery Module [27], in the process of 
being made available, will improve the situation dramati- 
cally. Some of the relationships from that module were used 
in calculating electric propulsion weights. 
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fare (ASW). It is armed for that purpose and has equipment 
in keeping with its size. The baseline ship is described in 


more detail below. 


6.2.1. Margins 


A naval ship design has margins in weight, vertical 
center of gravity (KG), space, ship service electrical gen- 
eration, propulsion power, accommodations, and structural 
strength which allow for equipment, mission, and system 
growth over its projected thirty year life. Without these 
margins, the ship would be difficult to modify, because 
whatever might be added in these areas would have to be paid 
for by a removal. For example, if a 50 ton radar system 
were added, the original 40 ton radar system and 10 tons of 
fuel might be removed to leave the ship at its original 
weight. With margins, the 10 tons of fuel might not be 
removed. 

Margins are typically split into Acquisition and Serv- 
ice Life allowances. Acquisition margins recognize the fact 
that ship specifications change over the design cycle and 
during construction. For example, the fourth ship built may 
have a different weapons system than the first, with a dif- 
ferent electrical requirement. If the electrical generation 
plant had to be changed during construction to accommodate 
the new weapons system, the total cost might be prohibitive. 
If an Acquisition margin is built into the original design, 
this may not occur. A Service Life margin makes allowance 
for configuration changes over the life of the ship. 

The ASSET program uses margins when synthesizing a 
ship. The margins are under operator control. The margins 
Suggested by Goddard and used in this analysis are listed in 
the table.18 


19. Goddard, C.H., "A Methodology for Technology Charac- 
terization and Evaluation for Naval Ships," S.M. NAME AND 
O.E. thesis, MIT, 1985, p. 31. 
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Table 23. Recommended technology assessment design margins 


for a monohull surface combatant 


Acquisition Service Life 
Weight@ 12.5% of SWBS 1-7 10% 
KG 12.5% of KG 1-7 20" #432 
Space 0 (no excess volume) 0 
Electrical]> 20% 20% 


Propulsion powerc 10% total EHP prior to prelim body plan 
8% prior to self-propelled model tests 

Accommodations Accom = 1.1 x ship manning at delivery 

Strength 2.24 KSI of marginal stress at delivery 
(Max primary stress for hull material) 


Notes: 

a. The service life weight margin applies only to naval 
architectural limits of the ship (reserve buoyancy, 
stability, structures), not to the final design weight. 

b. In sizing the electric plant, the calculated maximum 
electric load plus these design margins shall be met with 
one generator out of service. The remaining generators 
shall not be loaded in excess of 90%. Note that the service 
life margin is not applied to SWBS group 200 which would be 
expected to remain stable over the life of the ship. 

c. Performance requirements (sustained speed, endurance 
range) are met at delivery full load displacement. 


6.3. Philosophy of effort 


The nature of this technology characterization required 
that certain limits be imposed on the total effort. (If the 
Naval Sea Systems Command were to do this study, many people 
would simultaneously be employed to investigate every 
detail.) Some items were fixed, some were allowed to float 
with the design. 

The hardest item to handle is volume. There are very 
few ways of adjusting volume as easily as weight is ad- 
justed. One way is through the use of Marginal Volume Fac- 
tors, which equate a weight penalty with every increase in 
volume. (See Howell [26]. ) 

The differences among transmission systems appear 


primarily in the machinery spaces and fuel tanks of a ship. 
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ASSET can handle the tankage changes without aid. However, 
it does no equipment arrangement or space analysis inside 
the machinery spaces, leaving that for arrangement experts 
to do off-line. It was decided to keep the machinery space 
volume the same in all ships of this study, no matter how 
that volume was divided into spaces. Changes in equipment 
volume will be noted in the analysis and left for the advan- 
tage of others in machinery space rearrangement. 

Moderate to high technical risk has been accepted in 
specifying equipment cooling. Current densities at the 
limits of cooling technology are used, assuming that liquid 
cooling of both the stator and rotor can be performed. 
Lower current densities would be required if such cooling 
were not possible, resulting in slightly larger, heavier, 
less efficient machines. This last statement was proven 
during the course of the thesis research, as the first 
(chronologically) current densities used were two-thirds 
larger than those listed herein. 

Other risk areas include the use of advanced vacuum 
switchgear and the assumed efficiencies of reduction gears 
and power converters. These are low risk items; the 


technology is well understood and commercially available. 


Table 24. Ship design items held constant during analysis 


Endurance speed 20 knots 

Endurance rangea 5500 nautical miles 
Machinery box volume 109,670 ft3 
Installed horsepower 51,550 hp 

Payload weight, volume, and electrical requirements 

Length 425 ft. 

Beam 55 ft. 

Ship electrical load 2030 kW (24 hr. avg) 
Ship molded lines 

Manning 


Deckhouse and superstructure geometry 
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Table 25. Ship design items allowed to float during analysis 


Maximum and sustained speeds 

Endurance fuel load 

Ship full load displacement 

Ship draft 

Ship resistance and powering 

Ship arrangement outside machinery spaces 


Note: 

a. Endurance range was allowed to float as a comparison 
between two electrical transmission ships and the mechanical 
transmission baseline ship. It was held constant for the 


rearranged ship. 


88 





Most of the ship synthesis has been left to ASSET 
though parameters from Goddard were used where possible. 
ASSET designs a reasonable, generic ship with good seakeep- 
ing characteristics. As shown above, some ship characteris- 
tics were frozen to ensure transmission comparisons were not 
performed with different ships. Leaving the ship synthesis 
to ASSET allowed concentration on the specifics of the 


propulsion plant. 


6.4. Baseline ship 


The baseline ship has a mechanical transmission, i.e., 
two power trains consisting of a gas turbine, clutch and 
coupling, reduction gear, shafting, and propeller. There 
are two machinery rooms, each containing one gas turbine. 
The gas turbine used as the model is the General Electric 
LM-2500, rated at 25,775 brake horsepower. This is a very 
common marine gas turbine. It and its predecessors are 
powering the latest classes of naval combatant, such as the 
DD963, FFG7, and DDG51. 

The locked-train double-reduction gears are reversible, 
allowing the use of fixed-pitch propellers. There is no 
mechanical cross-connect allowed between the shafts. Except 
for the power level, this is the gear system being employed 
on the DDG51 class. Gears of this sort are about one per- 
cent inefficient29 per reduction stage. Since these are 
double-reduction gears, an efficiency of 98% was used. 

An endurance speed of 20 knots has been specified. 

This is in keeping with standard fleet practices. An en- 
durance range of 5500 nautical miles permits ocean crossings 
without refueling. The lack of a cross-connect capability 
between the two shafts means at least two gas turbines will 


20. Inefficiency = 1 - efficiency. Information on stage in- 
efficiency is from a conversation with Mr. Samuel Shank, the 
author of the ASSET Enhanced Machinery Module [27]. 
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be on-line during endurance cruising. This is inherently 
inefficient; some operators alleviate the inefficiency by 
declutching one shaft and free-wheeling that propeller, 
reaching desired speeds by loading up the operating gas tur- 
bine. 

A measure of initial static stability is the ratio 
GM/Beam. GM is the vertical distance between the center of 
gravity and the metacenter of the ship. Typical values for 
this ratio are 8-10%. A lower value (6.5%) has been ac- 
cepted for purposes of comparison with the variants. A 
large ship redesign effort would have been necessary to 
bring GM/B into a better range. 

The electrical generation plant consists of three gas 
turbines, each driving a 2000 kW generator. The data used 
for the gas turbines was taken from the Detroit Allison 501K 
turbine-generator set used aboard the DD963 ships. 

Both the deckhouse and main hull are constructed of 
High Tensile Strength (yield stress = 50,000 psi) steel. 
Active stabilizing fins and a sonar dome are included in the 
design. The payload is listed in Table 26. A coarse layout 


of the machinery spaces is shown in Figure 20. 


Table 26. Payload for baseline and variant ships 


FFG7 Command and control suite 
Satellite, UHF, and HF communications 
SLQ32V3 electronic countermeasures 

NIXIE acoustic countermeasure 

SPS-49 two-dimensional air search and tracking radar 
SPS-55 surface search and tracking radar 
SQR-19 towed array surveillance system 
MK92 missile and gun fire control system 
Harpoon fire control system 

LAMPS III helicopters and support system 
JP-5 aviation fuel 

MK32 over-the-side torpedo system 

MK13 guided missile launching system 
MK75 76mm gun 

Close-In-Weapon-System 

Small arms 

Appropriate ammunition and reloads 
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Baseline ship subdivision and machinery arrange- 


Figure 20. 
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6.5. Backfit ship 


jts¢ ane mechanical transmission of the baseline ship 
were replaced with an electrical transmission without chang- 
ing the subdivision of the ship, it would be as if an older 
ship had been updated, or "“backfitted,'" with new technology. 
This is the idea behind the first variant ship. 

Propulsion motors and generators were added into the 
propulsion plant, in the original spaces. Shafting still 
runs from the machinery spaces to the propellers. A reduc- 
tion gear is still necessary for higher rpm propulsion 
motors, but electrical cross-connect may improve the en- 
durance fuel efficiency. 

Some rearrangement within the machinery space is neces- 
Sitated by the backfit. The machinery spaces in the 
baseline ship are not long enough to contain the stack-up 
length of gas turbine, generator, motor, and reduction gear, 
Without "folding" the power train. This may be accomplished 
by changing the design of the reduction gear or by placing 
the gas turbine and generator (which require a mechanical 
connection) side-by-side (or transversely) with the motor 
and reduction gear (another mechanical connection), using 
transmission line to electrically connect them. Since there 
are a variety of ways to rearrange the machinery box, and 
since the chosen method has no effect on the analysis, the 
rearrangement was not specified. 

The propulsion generators and motors may operate at 
different rotational speeds, and therefore different elec- 
trical frequencies, so power converters must be used between 
them. Power converters change the frequency of the power 
being transferred between the generators and motors, through 
the use of cycloconverters or thyristors. They add another 
inefficiency to the transmission. A reasonable estimate of 


the efficiency of an 18 MW power converter is 97%.21 


21. Professor John Kassakian, MIT, private communication. 
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6.6. Rearranged ship 


The second variant ship takes advantage of the benefits 
of machinery space rearrangement. The same machinery space 
total volume is preserved but split into five spaces. The 
propulsion motors are placed very near the propellers, 
resulting in much shorter shafting runs and an increased 
GM/B ratio. The decrease in shaft length means a decrease 
in shaft weight and more space available outside the 
machinery box. The shafts previously ran through shaft al- 
leys that may be returned to other uses. This rearrangement 
is almost certainly not the optimum one and can be improved 
in the sense of space efficiency. It is an arrangement, 
however, that can demonstrate the benefits to be expected of 
a ship designed for electric drive. Power converters again 
connect the generators and motors. Transmission line forms 
the connections, at a much lower weight than shafting. 

Many other choices in large components are possible for 
this rearranged ship. For example, three propulsion gener- 
ators and gas turbines driving two or four propulsion motors 
might have been chosen. The number of prime movers was 
retained from the mechanical baseline ship, however, to make 
the comparison of transmissions realistic. Too many changes 
might have obscured the fundamental differences in 
efficiency, weight and volume. 

Both geared and direct-drive propulsion motors were 
used in this variant. When geared motors are used, the 
reduction gears are also placed near the propellers. A 


coarse layout of the rearranged ship is shown in Figure 21. 


6.7. Weight and volume algorithms 


Few components in this thesis are exact commercial 
models. The weight, volume, and other characteristics are 
taken from those for which data was available. The equa- 


tions for shafting and reduction gears were taken from the 
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ASSET theory manuals [28], while the switchgear and braking 
resistor equations came from the ASSET Enhanced Machinery 
Module [27], which is not yet generally available. These 
equations represent much study by ship and equipment desig- 
ners, incorporating equipment which is commercially avail- 
able. Where possible, the ASSET equations were verified 
against other studies and actual equipment [29, 30, 31]. 
For example, the machinery in the FFG7 was used as a model 
and verification for reduction gears and shafting. 

Appendix E gives more explanation, as well as presenting a 


computer program used to generate weight and volume figures. 
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Layout of machinery spaces on rearranged ship 


Figure 21. 
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Chapter Seven. Machine design and system synthesis 


7.1. Machine matrix 


During this study, synchronous, permanent magnet, and 
induction machines were designed at shaft speeds of 180, 
1800, 2400, 3000, 3600, and 7200 rpm. The number of pole- 
pairs for the 180 rpm machines ran from one to twenty-five. 
For the higher rpm machines, from one to six pole-pairs were 
used. If every machine could be a generator or a motor, a 
165x165 machine matrix results. From these 27,225 possible 
combinations, two were chosen and input to the two variant 
ships for synthesis in ASSET. 

These particular rpms were chosen partially because of 
the choice of the gas turbine. The LM-2500 operates ata 
full-load speed of 3600 rpm, making multiples and “nice” 
fractions of that speed desirable. <A 3600 rpm, two-pole 
machine has a synchronous frequency of 60 hertz, the stand- 
ard in the United States. A four-pole machine at 1800 rpm 
is also a 60 hertz machine. If "nice" frequencies result 
from rpm choices, results may not be obscured. The low rpm 
machine is tied to propeller rpm. For the baseline ship, 
maximum speed propeller rpm is 170 rpm. A ten rpm allowance 
for "battle override’ gives a requirement for 180 rpm. 

The reduction in the number of machine combinations is 
a bit more difficult to explain. First, it was observed 
that reduction gears add greatly to the weight and volume of 
the transmission and detract from its efficiency. Second, 
induction generators are notoriously difficult to control. 
It was then decided that generators would not be induction 
machines and any generators used would operate at the same 
shaft speed as the gas turbine, eliminating a possible 
reduction gear. The matrix then measured 12x165 and had 
1980 combinations. 

From this point on, the decision theory espoused in 


Schweppe and Merrill [13] was used, specifically using "knee 
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sets" to eliminate uncompetitive designs. The considera- 
tion, or Figure-of-Merit, was always minimum weight, volume, 
and inefficiency. These external characteristics are those 
“seen by a ship design (since propulsion current and volt- 
age were specified in terms of turns -- fertile ground for 
another tradeoff study). Initially, individual motors of 
all of the higher rpms were considered together to select 
the best geared propulsion motors. Since the weight and 
volume of a reduction gear varies with shaft rpm, the weight 
and volume of the reduction gear was included with that of 
the machine to select the best machine-gear system. Direct- 
drive machines were also selected. The 3600 rpm synchronous 
and permanent magnet machines were also considered 
separately as propulsion generators. 

Once the initial selection of machines was made, the 
matrix measured 5x11, or 55 combinations. All of the com- 
binations were plotted in knee curves that showed the 
volume, weight, and inefficiency of each transmission. The 
volume and weight of shafting, braking resistors, cooling 
systems, power converters, and the inefficiency of the power 
converters were common to all combinations and were not in- 
cluded at this level. The inefficiency of any reduction 
gears was included where appropriate. Three of the gener- 
ators were synchronous machines and two were permanent mag- 
net machines. Of the motors, two were 1800 rpm synchronous 
machines, two were 1800 rpm permanent magnet machines, three 
were 1800 rpm induction machines, and two were 180 rpm 
synchronous machines. The 1800 rpm motors clearly dominated 
the higher rpm machines, largely because of the differences 
in reduction gear weight and volume. 

Since there was no single dominating combination, ten 
of the 55 combinations were selected. These ten were among 
the best at least twice on the knee curves. These ten com- 
binations were composed of only synchronous machines. 

Programs to calculate off-design-point efficiency were 


written. Each motor and generator was evaluated at the 
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power level and rpm appropriate for the sustained and en- 
durance speed conditions of the ship, using the delivered 
horsepower (DHP) and propeller rpm taken from the ASSET syn- 
thesis run on the baseline mechanical ship. Some of the 
combinations had very low endurance efficiencies. There was 
a sharp division evident between the geared (lower weight, 
much higher volume) and the direct-drive systems. 

The ten combinations, with their maximum, sustained, 
and endurance speed transmission inefficiencies, were again 
made the subject of knee curves. A simple scoring scheme 
was devised to rank the combinations according to their 
grouping on these last knee curves. If a combination was in 
the best group on a particular knee curve, it was given two 
points for that curve. If it was in the second best group, 
it received one point. If it was in neither the best or 
second best group, it received no points. When the scores 
were totaled, two combinations stood out. One was a geared 
drive system and one a direct-drive system. These two com- 
binations were used in both of the variant ships and are the 


subject of the next chapter. 


7.2. Knee curves 


Figures 22 through 24 show the knee curves for the 
propulsion generators. The first letter of the generator 
ID indicates whether it is a synchronous machine or a per- 
manent magnet machine. All of the generators are 3600 rpm 
machines. The generators selected were SB (four poles), SC 
(six poles), SD (eight poles), PB (four poles), and PC (six 
poles). Generator SA was not selected because of its poor 
showing on the volume-weight curve, even though it was com- 
petitive on the volume-efficiency curve. 

The 180 rpm, direct-drive propulsion motor curves are 
in Figures 25 through 27. They were not combined with 
geared motors because one of the points being investigated 


was whether or not geared motors were "better" than direct- 
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drive motors. The clumping of the machines necessitated 
other graphs on different scales to distinguish between the 
machines. Machines 1-25 are synchronous, 26-50 are per- 
manent magnet, and 51-56 are induction machines. The number 
within a group indicates the number of poles in the machine, 
e.g., machine 32 has (32-25)x2=14 poles. Valid designs with 
over twelve poles were not achieved for the induction 
machines. Machines 3, 8, 11, and 16 were selected for fur- 
ther work. These are all synchronous motors. 

The higher rpm motors were combined with their reduc- 
tion gears to form system knee curves. In all of these knee 
curves, machines 1-6 are 1800 rpm, 7-12 are 2400 rpm, 13-18 
are 3000 rpm, 19-24 are 3600 rpm, and 25-30 are 7200 rpm. 
Figure 28 is the volume-efficiency curve for synchronous 
machines, showing the distinct grouping of the machine-gear 
systems due to the high volume of the reduction gears. Note 
the high values of the permanent magnet machines in the 
volume-weight curve Figure 29. The clumping of induction 
motors around the low inefficiencies is shown in the weight- 
efficiency curve of Figure 30. From these curves, the 
motors on page 97 were selected. 

The initial motor and generator combinations were made 
and plotted on more knee curves (Figures 31 to 33). On 
those curves, the high-volume or high-weight nature of the 
combinations can be seen. Since the multiple-attribute 
decision theory embodied in knee curves does not say how to 
select between high-volume or high-weight, the best of each 
were selected. Combinations 1, 2, 8, 9, 12, 13, 19, 20, 30, 
and 31 were chosen. The off-design-point efficiencies were 
calculated and all of the information was plotted. 

Figure 34 is a bar-graph of the maximum, sustained, and en- 
durance speed transmission efficiencies of the various com- 
binations, including the reduction gears, if any, and power 
converters. The final combination knee curves are sum- 
marized in the scoring scheme of Table 27, which was ex- 


plained on page 98 . Combinations 12 and 20 were chosen to 
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use in the ships of the study. 

This is a good method to choose among the possible 
machines. During the course of this thesis, the above path 
was followed through several complete iterations and a few 
partial ones. As stated in Schweppe and Merrill, knee 
curves serve very well to eliminate uncompetitive options, 


allowing concentration on the better ones. 


Table 27. Final combination knee curve scores 


Combo ID Firsts Seconds Total 


1 3 1 7 
2 3 0 6 
8 1 0 Z 
9 0 1 1 
EA 6 i 13 
13 4 0 8 
19 if 1 3 
20 3 4 10 
Sil 0 2 2 


Conclusion: test combinations 12 and 20 
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Figure 22. 
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Figure 23. Curve of volume-weight for 3600 rpm generators 
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Figure 24. Curve of weight-efficiency for 3600 rpm gener- 
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Figure 25. Curve of volume-efficiency for 180 rpm motors 
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Figure 26. Curve of volume-weight for 180 rpm motors 
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Figure 27. Curve of weight-efficiency for 180 rpm motors 
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Figure 28. Curve of volume-efficiency for geared motors 
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Figure 29. Curve of volume-weight for geared 
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Figure 30. Curve of weight-efficiency for geared motors 


SYSTEM KNEE CURVE 


INDUCTION MACHINES, 25775 HP 
0.08 


0.07 
0.06 
0.05 


0.04 


(1 — EFFICIENCY) 


eS 


0.02 


0.01 





400 


0 100 200 300 
(Thousands) 
WEIGHT (KG) 


109 





(1 — TOTAL EFFICIENCY) 


Figure 31. Volume-efficiency curve for initial PM and PG 
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Figure 32. 
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Figure 33. Weight-efficiency curve for initial PM and PG 


combinations 
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Figure 34. Final combination transmission efficiencies 
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Chapter Eight. Analysis 


This chapter presents an analysis of the ships with 
electric transmissions. Standard naval architectural 
methods have been used to observe and comment on the variant 
ships, comparing them to the mechanical baseline ship. Con- 
clusions and recommendations follow the analysis. 

The names used to describe the various ships imply 
their internal arrangement, equipment, and ASSET Design Mode 
Indicator (DMI). The two DMIs used were ENDURANCE, when en- 
durance range was held constant at 5500 NM, and FUEL WT, 
when the usable fuel weight was held constant at 996.6 LT.22 


The ship names are as follows: 


MECH 23 BASELINE: Baseline, mechanical transmission 

ELEC 23 BASELINE 12: Backfit ship, geared motors 

ELEC 23 BASELINE 20: Backfit ship, direct-drive motors 

NEW MR ELEC 12: Rearranged ship, geared motors 

NEW MR ELEC 20: Rearranged ship, direct-drive motors 

CONSTANT FUEL ELEC 23 BASE 12: Backfit ship, geared 
motors 

CONSTANT FUEL ELEC 23 BASE 20: Backfit ship, direct- 
drive motors 


8.1. Direct effects 


The direct effects of an electric transmission are the 
changes in weight and volume of the propulsion system, as 
well as the transmission efficiency. Included are the 
weight and volume of the propulsion motors and generators, 
transmission lines, cooling systems, switchgear, power con- 
verter, exciter, braking resistor, any reduction gears and 
their associated gear oil, and shafting. These items are 
listed in Table 31. A positive difference from the baseline 


ship means a heavier and/or larger ship. 


21. Professor John Kassakian, MIT, private communication. 


22. Not all fuel in a ship is usable. There are nooks and 
corners of fuel tanks that are inaccessible to the fuel sys- 
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Note that the only variant that has a lower direct 
weight effect than the baseline is NEW MR ELEC 12. The ac- 
cumulation of weight increases in the others makes them 
heavier, while NEW MR ELEC 12 has lower motor and shafting 
weight than the rest. Geared drive is always lighter than 
direct-drive, largely due to the high weight of the direct- 
drive propulsion motors. With respect to volume, direct- 
drive is always smaller than geared drive, because of the 
reduction gears. All electric transmissions are larger than 
the mechanical baseline, but the smallest variants, within 
motor type groups, are the rearranged ships. 

Table 28 contains the maximum, sustained, and endurance 
speed transmission efficiencies of the two generator-motor 
combinations. Note that the off-design-point efficiency of 
the direct-drive combination is significantly lower than the 
geared combination, even though it does not have the added 
inefficiency of reduction gears. This is in large part due 
to the poor efficiency of the slowly rotating direct-drive 


motor at the endurance speed. 


Table 28. Transmission efficiencies 


Combo Maximum Sustained Endurance 


— =e ame cee cme os ame cee cee cee com ose ome cme ome ome eee eee ee eee eee ee ee eee 


12 0.9307 0.9266 0.8817 geared combination 
20 0.9209 0.9093 0.7754 direct-drive combination 
Endurance efficiency with one generator driving two motors. 


Table 29. Propulsion generator efficiencies 


— oe oe oe es ee eee ee ee ee eee eee eee 


SC 0.9891 0.9870 0.9737 geared combination 
SC 0.9891 0.9872 0.9768 direct-drive combination 


Table 30. Propulsion motor efficiencies 
PM ID Maximum Sustained Endurance 


S4 0.9898 0.9876 0.9526 geared combination 
SL8 0.9598 0.9496 0.8184 direct-drive combination 
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The above tables show the efficiencies of the motor and 
generator used in two particular combinations. The 
efficiency of the motor, reduction gears, and power con- 
verter have a direct effect on the efficiency of the gener- 
ator, as they change the loading point of the generator. 
Generally, motors and generators are more efficient when 
they are loaded more closely to their design point. (The 
same is true of gas turbines.) The inefficiencies of the 
power converters and reduction gears, if any, are included 


in the transmission efficiencies. 
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Table 31. Direct volume and weight effects 
CONSTANT CONSTANT 
FUEL FUEL 
MECH ELEC ELEC NEW NEW ELEC ELEC 
23 23 23 HR HR 23 25 
BASELINE BASELINE BASELINE ELEC ELEC BASE BASE 
12 20 12 20 12 20 
Electric Propulsion weights: in LT 
PMs 0 18.08 102.94 18.08 102.94 18.08 102.94 
PGs 0 10.27 10.27 10.27 10.27 10.27 10.27 
Trans lines 0 0.18 0.18 0.77 0.77 0.18 0.18 
Cooling sys 0 3.98 o.90 5.98 398 3.98 D960 
Switchgear 0 1.56 1.56 1.56 1.56 1.56 1.96 
Power converter 0 7.16 7.16 7.16 7.16 7.16 7.16 
Exciters 0 3.90 3.90 BRAY 3.50 3.90 sca 
Braking resistor 0 10.00 10.00 10.00 10.00 10.00 10.00 
Reduction gears 78.90 41.19 0.00 41.19 0.00 41.19 0.00 
Shafting 69.00 66.84 66.84 41.72 41.72 66.84 66.84 
W298 (op fluid) 14,7 14.7 0 14.7 0 14.7 0 
Direct effect 162.60 179.46 208.43 154.93 183.90 179.46 208.43 
Diff fa Baseline 0 16.86 45.83 -7.6] 21.30 16.86 45.83 
Electric Propulsion voluses;: in cubic feet 
PMs 0 108.86 909.13 108.86 909.13 108.86 909.13 
PGs 0 61.49 61.49 61.49 61.49 61.49 61.49 
Trans lines 0 20.9 20.9 6.39 6.39 0.33 0.33 
Cooling sys 0 200 200 200 200 200 200 
Switchgear 0 70.2 NOs 10.2 10.2 70.2 70.2 
Power converter O 1089.73 1089.73 1089.73 1089.73 1089.73 1089.73 
Exciters 0 245.02 245.02 245.02 245.02 245.02 245.02 
Braking resistor O 1422.88 1422.88 1422.88 1422.88 1422.88 1422.88 
Reduction gears 2731.04 1425.75 0.00 1425.75 0.00 1425.75 0.00 
Shafting 917.71 JO1.a 901.5 306.06 306.06 901.5 501.5 
Direct effect 3248.75 9146.33 4520.85 4936.38 4310.90 9125.76 4500.28 
Diff fa Baseline 0.00 1897.58 1272.10 1687.63 1062.15 1877.01 1251.55 


8.2. Indirect effects 
Indirect effects are again composed of weights and 
volumes, but these are the ripple effects of the propulsion 


For example, if a transmission is 
it should be expected 


that less onboard fuel would be needed to achieve the same 


system through the ship. 


more efficient at endurance speed, 
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endurance range as a less efficient transmission. This is 
indeed the case. Another important indirect effect is the 
change in full load displacement, which is tied to the 
powering characteristic of the ship.23 The following tables 
list the characteristics and indirect effects of the various 
ship configurations. 

Every electric transmission had lower maximum and sus- 
tained speeds than the baseline ship (by about 0.43 knots), 
but also lower EHP requirements at those speeds. The lower 
EHPs are a reflection of lower drafts (less resistance). 

The lower speeds show that the transmission efficiencies of 
the variant ships are lower than the mechanical baseline 
ship. There are more components in the electrical power 
trains, hence the lower efficiencies. The speed difference 
of 0.43 knots may be regarded by some as significant; it is 
about the speed reduction to be expected by a fouled bottom. 

The endurance range of all ships except those with con- 
stant fuel load is 5500 NM. The fuel load in the others 
varies greatly, showing the benefit of electrical cross- 
connection. In the two constant fuel ships, there was an 
increase in the endurance range of 1350 and 1400 NM, respec- 
tively, for the geared and direct-drive transmissions. This 
is an indication of fuel cost savings from the electric 
transmission. If a ship refuels every three steaming days 
(receiving a third of its tank capacity), steams 100 days 
each year, and fuel is priced at about $18 per barrel, this 
represents about a $600,000 savings per ship per year. 

The initial static transverse stability of the variants 
was degraded by the change in propulsion equipment. As pre- 
viously stated, the 6.5% GM/Beam ratio of the mechanical 
baseline ship is not as large as desired for an actual ship, 


but provided a benchmark to measure relative changes. 


tem. Typically, 95% of the onboard fuel is usable. 


23. For the same molded lines, ships with higher displace- 
ments will have greater wetted surface areas and higher 
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Table 32. General ship characteristics 


aes @ ae @ ae 8 Se 2 2 2S SS SS SS SS SF SF SF SSF F888 SF OK FSF SET STF SFT STF EFSF EFSFTF ST SEETSTSSISSFISESEISSTF@ SOT SOS SETHWSETSBFI SBEF EB FO FEOF FFE SBE ER GBs 


LBP 
Ship disp FL 


Diff fa Baseline 


Nav beaa 

Nav draft 
Depth sta 10 
Accomaodations 
GNt/B 

M1 

Full load K6 
LCB/LBP 
LCF/LBP 

Wetted surface 
Cp 

Waterplane area 
Hull voluae 
Dkhs voluae 
Total voluse 
Endur fuel vol 
Machy box vol 


MECH 
23 


BASELINE 


420 
0489 

0 

ah) 
16.44 
38 

272 
0.065 
950.01 
23.23 
0.516 
0.969 
23347.7 
0.4 
17641 
618504 
121928 
740032 
44345 
199671 


EEGE EEG 

23 23 
BASELINE BASELINE 

12 20 
429 425 
0297.9 47.1 
-187.5 =67.7 
Joe 04 ud. 02 
16.07 16.31 
38 38 
272 272 
0.059 0.098 
973.41 “Rie 
23.87 23.79 
0.914 0.919 
0.960 0.969 
23010 = 23230.5 
0.596 0.598 
17973.9 17624.8 
618504 618504 
121928 =: 121528 
740032 = 740032 
36092 39647 
109671 109671 


17939.8 
618504 
121528 
740032 

36367 
109624 


NEW 
MR 
ELEC 
20 


0.064 
967.39 
23.93 
0.914 
0.069 
2313.7 
0.997 
17602 
618904 
121528 
740032 
39396 
109624 


CONSTANT 
FUEL 
ELEC 

23 
BASE 
12 


0.063 
990.07 
23.34 
0.916 
0.960 
23346.6 
0.4 
17640.9 
618304 
121928 
740032 
44364 
109471 


CONSTANT 
FUEL 
ELEC 

23 
BASE 
20 


0.916 
0.969 
23425.8 
0.401 
17693.2 
618504 
121528 
740032 
44364 
109671 





Table 33. Powering 


CONSTANT CONSTANT 


FUEL FUEL 
MECH ELEC ELEC NEW NEW ELEC ELEC 
23 an 23 MR HR 25 23 
BASELINE BASELINE BASELINE ELEC ELEC BASE BASE 
12 20 12 20 12 20 


om oe oy a op ae oP op op oe at eet oe Oe ge oe et et om em eet oe oe oe ee ee ee eS eS ES Se ee SE ESE SE SSS STB es se Se Oe eT ee ewe es ee SF eS eee 8 ee ee ees eee es |S 


Powering: 
Vmax oon 1S 29.02 28.83 Ale] 28.9 28.82 28.7 
EHP @ Vaax 34821 33131 3274} 33194 32763 33064 32700 
Vsus 27.94 27.82 27.4 27.88 27.46 27.65 27.48 
EHP @ Vsus 28145 26684 26149 26/700 26181 26436 26136 
Endurance 2000 2000 2500 3300 Jd00) 6550 6084 
EHP @ Vend 6851 6452 6772 6597 6498 685i 6905 
HP i 21590 31590 v1590 v1390 21990 wl oo0 31990 
kWi 6000 6000 4000 §000 4000 4000 6000 
Avg 24 hr load 2030 2030 2030 20350 2030 2030 2030 





In the constant fuel ships, the high weight of the 
propulsion motors and propulsion motors and propulsion gen- 
erators, combined with a smaller amount of vertically-lower 
fuel reduced GM/B by over one-half percent, a not incon- 
siderable amount. For the rearranged variants, GM/B 
decreased less than the two previous ships. They also have 
less low fuel, but the propulsion generators are lower than 
in the backfit ships and the propulsion motors are very much 
lower. Only avery small decrease in GM/B was seen in the 
constant fuel ships because the constant fuel load compen- 
sated for the increased high weight of the electric trans- 
mission. The longitudinal metacentric height, GMl, in- 
creased for all variants, though it seems it should have 
decreased with the decrease in waterplane area and draft. 
The machinery space volume was the same for all ships. 

There were no big surprises in the area of weight. The 
structural weight (W100) encloses the same volume in every 
ship, so it was about constant. The propulsion plant weight 
(W200) varied with the type of transmission. Weight groups 
W300, W400, W500, W600, and W700 were virtually identical in 
every ship, and the variable loads were dominated by the 
change in fuel weight. The Design and Builders Margin is a 
function of the light ship weight (summation of W100 through 
W700), so the margin weight moved with the light ship 
weight. The miles-~per-gallon figure of NM/LT of fuel shows 
the endurance efficiency of electrical cross-connection. 

Only a few comments need be made regarding Table 35. 
The structural weight fraction shows the changes in full 
load displacement, remembering that the W100 weights were 
all about the same. The same may be said for the weight 
fraction of the W300 through W700 groups and payload weight. 
Higher propulsion plant weights in the variant ships drove 
up the W200 fraction, except for NEW12. The fuel weight 
fraction shows the same behavior as the miles-per-gallon 


figure. 
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Table 34. Ship weights 


eae e222 eB SS 8S 88 8S OS SS SS SS SS SS SS SS SF S SF S SS SS SS SF 8S SO GSC S&C SF SS SS SS SS OC TC SS SO SF SF SS SC C8 SSF SS SS SS ee | Ce es = 


Weight: 

DFM onboard 
Usable fuel wt 
Diff fm Baseline 
NM per LT fuel 
Payload weight 
#100 

W200 

W300 

#400 

W500 

Wb00 

W700 

Loads 

D&B margin 
Disp FL 


MECH 
23 


BASELINE 


1049.1 
996.6 
0 
Jedd 
Type 
1684.8 
343.6 
2aoud 
$02.2 
615.4 
426.6 
95.7 
1317.1 
443.1 
2485.2 


ELEC ELEC 

23 23 
BASELINE BASELINE 

12 20 
865.3 Poved 
822 890.6 
see ea (Tees 
6.49 6.18 
aie w71.2 
1486 1692 
sep 387.1 
236.9 256.9 
$02.2 302.2 
613.9 614.3 
426.6 426.6 
rca 70.9 
1120590 119322 
444.1 469.3 
S24 ied aural 


840 


-189.1 
6.73 
af1.2 
1660.4 
326.1 
236.3 
$02.2 
613.9 
426.6 


1115.6 
437.7 
0254.7 


951.6 


“H7.8 
6.21 
my eg 
1644.6 
361.3 
236.5 
302.2 
614.2 
426.6 
foe? 
1187.2 
442.9 
o3u5. 4 


FUEL 
ELEC 
25 
BASE 
12 


1049.1 
996.6 
0 
6.97 
a/1.2 
1686.2 
AHL AE, 
236.5 
$02.2 
615.4 
426.6 
eg 
1305.1 
464.3 
0484.1 


CONSTANT CONSTANT 


FUEL 
ELEC 
23 
BASE 


1049.1 
996.6 
0 

6.11 
vil.2 
1692.1 
387.1 
236.9 
302.2 
615.4 
426.6 
95.9 
1304.9 
469.9 
9090.2 





Table 35. Naval architectural analysis indices 


CONSTANT CONSTANT 


FUEL FUEL 
MECH ELEC ELEC NEW NEW ELEC ELEC 
23 18 us MR AR 23 23 
BASELINE BASELINE BASELINE ELEC ELEC BASE BASE 
12 20 12 20 12 20 


= ee OS SS ew SE Se SH SS SSS SS SS SS SS SS SSF SOS SS ses 2 SOS SCs S28 SC S28 SP Ss SF SS SS SS SO SP S28 S28 28S Sees e222 SS Se ose 


B/T 3.345 5.425 3.373 3.431 3.401 3.349 Saoc7 

GA1/LBP Zolod 2.290 2,237 2.505 Falls feed 2.220 
Wi00/Df1 0.3507 0.318 0.312 0,017 0.311 0.307 0.306 
8200/D1 0.063 0.066 0.071 0.062 0.067 0.064 0.070 
WS00/Dt1 0.043 0.045 0.044 0.045 0.044 0.043 0.043 
W400/DF1 0.005 0.097 0,034 0,038 0.094 0.039 0.099 
W500/Df1 0,112 0.116 0.115 0.117 0.115 0.112 0.411 
W600/DF1 0.078 0.081 0.079 0.081 0.080 0.078 0.077 
W700/Df1 0.017 0.018 0.018 0.018 0.018 0.017 0.017 
Wfuel /Dfl 0.182 0.199 0.164 0.156 0.165 0,182 0. 180 
Rpayload/Df1 0.104 0.108 0.103 0.109 0.107 0.104 0.103 
Wid/Df1 0,240 0.212 0.220 0.213 0.222 0.238 0.236 
Vab/Vtot 0.148 0.148 0.148 0.148 0.148 0.148 0.148 

HPi /Df1 9.598 Flot 7.016 7.848 9.629 7.400 Weaez 

HPi tVaax /DF) ei3e/6 86 282.57 8=«s 274.35 0©=— 286.47 0278.29 270.91 = 267.35 
kai /DF1 1.094 1155 1.108 1,146 heiZ) 1.094 1.085 
Nt/Df1 0.050 0.031 0.050 0,032 0,051 0.030 0.049 
Vtot/Dt1 134.914 139.695 136.610 141.370 138.236 134.941 133.816 
Wi00/Vtat 9.100 9.105 aided 0.026 7.045 0.104 ve122 
W200/HPi 14.930 13.287 16.821 14.170 15.700 15.291 16,821 
Vab/HPi 2.127 ele Z.127 Z.127 2.127 Pole) ZeV2T 
W300/kWi 0.039 0.039 0.039 0.039 0,039 0.039 0.039 
Vab/ (HPi+k#i ) 1.906 1.906 1.906 1.905 1,909 1.906 1.706 
W300/Vtot 1.863 1,837 1.859 1,837 1,899 1.863 1.863 
W600/Vtot 1.291 1.291 ara) 1.291 1,291 eee! 1,291 
Dti/Vtot 16.603 16.035 16.397 15.845 16.204 16.600 16.739 
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For this technology characterization, everything 
devolves to total values. What is the total effect on the 
ship, once the individual pieces are put together? Table 36 
gives the bottom line. The electric propulsion plants are 
larger and heavier than their mechanical drive cousin; 
however, the extra weight and volume are more than compen- 
sated by the savings in fuel weight and volume. If a ship 
is designed from the beginning to be an "optimized" electric 
drive ship, over 6300 cubic feet of volume and 250 LT may be 
saved. The savings might be used for other systems, to 
reduce the overall size and cost of the ship (maybe allowing 
a larger buy, since 30 ships times 250 LT is a 7500 LT 
ship), or to extend the naval architectural limits of the 
ship design. 

If a ship is backfitted with this technology, it is un- 
likely that tank volume can be recovered. However, the dual 
benefits of increased time-on-station and better fuel 
economy are realized. In this case, the choice between 
geared or direct-drive systems can be made by selecting the 
system with the most leverage, i.e., if the ship is volume- 
limited, use the lower volume direct-drive system (since the 
shafts are already in place). 

To put the volume and weight savings in perspective, 
note that 6300 cubic feet and 250 LT translates to twenty 
Tomahawk missile cells. The ship would be volume limited, 
with about 200 LT of weight savings still unused. This is a 
significant addition to the firepower of any ship, and the 


unused weight allows for ship growth. 
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Table 36. Total differences 


See e282 @2 SS SS S2 S28 2S SS BSF FE S88 PF V8 F888 SF SF SF8E8 S&S Sl FF S28 SS SSF SSF SSS FTF SE SB FSFE SFE SSE S88 S282 SF SF SF SB SCS FS SSF SSF S88 S88 S882 OS os 


Total voluses: 
Fuel voluse 
Fuel diff 
Prpln volume 
Proln diff 
Total 

Total diff 


Total weights: 
Df] difference 
prpin diff 
fuel diff 
other 


MECH 
23 


BASELINE 


44365 

0 
3248.79 
0.00 
47613.75 
0.00 


0.00 
0.00 
0.00 
0.00 


ELEC 
23 
BASEL INE 
12 


36992 
aIITS 
0146.33 
1897.98 
41738.33 
-9879.42 


-187.90 
16.86 

- 183.80 
-20.96 


ELEC 
23 
BASELINE 
20 


39647 
~4718 
4920.85 
1272.10 
$4167.89 
3445.90 


Bored 0 
45.83 
~111.60 
ENS. 


NEW 
AR 
ELEC 
12 


36367 
-7998 
4936.38 
1687.43 
41303.38 
-$310.37 


-250.30 
7.6] 
Sted e10 
“35.95 


NEW 
HR 
ELEC 
20 


39396 
-4969 
4310.90 
1062.15 
43706. 90 
-3906.85 


-131.60 
21.30 
-117.50 
~39. 40 


CONSTANT 
FUEL 
ELEC 

1B 
BASE 
12 


44364 

~{ 
0125.76 
1877.01 
49489.76 
1876.01 


-0.90 
16.86 
0.00 
~17.76 


CONSTANT 
FUEL 
ELEC 

23 
BASE 
20 


44364 

=I 
4500. 28 
1201.95 
48864.28 
1250.33 


49.20 
45.83 

0.00 
-0.43 





8.3. Closure 


8.3.1. Conclusions 


This thesis has demonstrated the usefulness of electric 
drive transmissions in reducing ship weight and volume. 
Electric drive transmissions are better than mechanical 
drive transmissions on a ship basis. They provide, besides 
the weight and volume advantages, substantial arrangement 
flexibility and the opportunity to use new technologies in 
the ship design arena. The technical risk associated with 
these different technologies is minimal, as there is much 
industrial experience with electric machines, advanced 
Switchgear, and the like. If the weight and volume reduc- 
tions are reinvested in the ship design through more optimum 
arrangements and subdivision, a substantially more efficient 
ship may be realized. Such a ship could successfully com- 
pete with the best of current ships. 

Small, light, high-power motors can be designed to a 
fair degree of detail with a computer optimization scheme if 
a meaningful objective function can be devised. For a ship 
system, the objective function should contain measures for 
volume, weight, efficiency, and relative cost (if a par- 
ticular material is significantly more expensive than other 
used). A steepest-descent scheme can be combined with a 
Monte Carlo scheme to quickly converge on the objective 
function. 

The use of electric drive, and its consequent electri- 
cal cross-connect, can reduce the endurance fuel load by as 
much as 17.5%. When used in combination with an improved 
machinery arrangement and subdivision, that percentage can 
rise to 18%. If the fuel load stays constant, the endurance 
range may increase as much as 25%. 

On both an equipment weight basis and a ship weight 
basis, systems composed of a direct-drive propulsion gener- 


ator (with the same shaft rpm as the prime mover) and a 
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geared propulsion motor are better than those systems using 
no gears. Regarding volume, a non-geared system has lower 
equipment volume but higher ship volume due to the lower en- 
durance fuel efficiency. 

Geared propulsion motors have better off-design-point 
efficiencies than those in direct-drive systems, primarily 
due to their higher rpm. A reduction in output power (in a 
motor) of 75% means only a few percent reduction in 
efficiency for a geared motor, while the same power reduc- 
tion means a 20% or more efficiency reduction for a direct- 
drive motor. 

Permanent magnet machines do not appear attractive for 
ship propulsion systems. They are both heavier and larger 
than candidate systems using synchronous, and, to a smaller 
extent, induction machines. Their low air gap flux density 
is the main detractor. Current permanent magnet materials 
cannot develop the energy product to compete with other al- 
ternatives, even though the NdFeB magnets are now in the 
marketplace. Induction machines may be useful as propulsion 
motors, but in this thesis they did not appear so. There- 
fore, ship propulsion generators and motors should only be 


synchronous machines. 


8.3.2. Recommendations 


The same modeling approach used in this thesis should 
be taken with variable reluctance machines (VRM). Although 
no VRMs have been built at ship propulsion power levels, it 
is not inconceivable that they could serve in such a 
Capacity. 

The induction machine model used here needs refinement, 
especially in the area of limiting maximum rotor current 
density. All of the machines need an analysis of their 
transient and dynamic characteristics. 

The recent advent of liquid hydrogen temperature super- 


conducting materials may signal an era where conventional 
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machines are overshadowed by the smaller, higher-flux 
machines possible with superconducting technology. However, 
if these new materials fail to provide the required current 
density, design of satisfactory machines may not be pos- 
sible. 

Integrated electric ship service propulsion plants may 
be beneficial additions to electric drive technology. Their 
influence on the systems suggested here may be an area of 


interest for future ship designs. 
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Appendix A. Definitions of machine variables and constants. 


bd 


BETA 


BR 
BR1 


BREM 
BRNGS 


BSAT 
CP 
CRHO 


CW 


dcore 
DCU 
DMAG 
doa 
dr 

ds 
eaf 


effcy 


ew 


freq 


magnet operating point flux density, Tesla. Used 
in permanent magnet machines. 
hysteresis loss factor of M19 magnetic steel. The 
figure 2.5 was used. 
air gap flux density, Tesla. Used only in per- 
manent magnet machine model. 
air gap flux density, 1.05 rms Tesla. 
residual induction, 7.89 kilogauss, of M19 mag- 
netic steel,at BSAT and Ti. 
remanence flux of NdFeB, 1.21 Tesla. 
weight percentage of rotor shaft bearings, 1.03, 
or three percent of rotor weight. 
max flux density anywhere, 1.5 rms Tesla. 
Stator coil pitch. The figure 0.8 was used. 
copper electrical resistivity, 1.724E-8 ohm- 
meters. 
copper weight, kg. 
density, 7.65e3 kge/m3, of M19 magnetic steel. 
back iron depth, meters. 
copper density, 8968.0 kg/meter?. 
density of NdFeB, 7.4e3 kg/meter?. 
overall machine diameter, meters. 
depth of rotor slots, meters. 
depth of stator slots, meters. 
p.u. internal voltage, used in syn only. 
efficiency, defined as 
output power 
Oe ee = i Sin Sa 
output power + ph + pe + i2r + i2rr 

effective weight of machine, a combination of 
weight, volume, and efficiency. Used as the ob- 
jective criteria for the optimization scheme. 
machine synchronous frequency, Hz. 


air gap dimension, meters. 
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gams 


gamr 


GMIN 
HC1 


vet 

i2 
aor 
i2rr 
jr 
jrnl 
js 
JSMAX 
ke 


ks[([] 
kv 


le 


ler 


lm 


loa 


lr 


lrat 


ls 


stator geometric factor, non-dimensional. Used to 
find convergence on active length. 
rotor geometric factor, non-dimensional. Used to 


find convergence on active length. 

minimum machine air gap, 0.002 meters. 

coercive force, 0.48 oersteds, of M19 magnetic 
steel, at BSAT and T1. 

induction machine primary current, amps. 
induction machine load current, amps. 

copper loss, watts, on stator. 

copper loss, watts, on rotor. 

full load rotor current density, amps/r. 

no load rotor current density, used in syn only. 
full load stator current density, amps/rm. 
maximum stator current density, 12e6 amps/meter 
efficiency weighting factor. 


magnet material cost weighting factor, defined as 


ee ee oe EsenoL material _ 
$ per pound of magnetic iron 

harmonic winding factors. 
volume weighting factor. 
machine active length. 
rotor winding space factor, used in induction 
machine. 
combined length of rotor end windings. 
radial dimension of permanent magnet. 
overall machine length, meters. 
ratio of rotor slot width to slot pitch. 
manent magnet machines, defined as ratio of magnet 
width to rotor “slot” pitch. 
length ratio, used in permanent magnet program to 


In per- 


calculate the effect of magnet overhang. 
ratio of stator slot width to slot pitch. 
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MAX TIP_SPEED maximum allowable rotor tangential velocity, 


minpwr 
MU 
MUR 


nr 
NU 


overhang 


Pp 
PC 


pe 
PF 


Pi 
ph 
PSI 


ripl 
r2pl 
rcv 
relpl 
RHO 


riv 


rpm 


RSF 


SCV 


Siv 


slip 


200 meters per second. 


minimum required power of the machine, horsepower. 
MU of air, 4ne-7 Henries/meter. 
NdFeB relative reversible permeability. The 


figure 1.05 was used. 

number of rotor bars in induction machine. 
anomalous loss factor of M19 magnetic steel. The 
figure of 2 was used. 

percentage of magnet overhang. 
number of machine pole-pairs. 
permanence coefficient of NdFeB. The figure 1.1 
was used. 

eddy current loss, watts. 
power factor of all machines. The figure 0.8 was 
used. 


3.141592654. 


hysteresis loss, watts. 

Trated/Tpullout for induction motor. The figure 
0.55 was used. 

rotor radius, meters. 

per length primary resistance, ohms/meter. 

per length secondary resistance, ohms/meter. 


rotor copper volume, meters?. 


per length Thevénin equivalent resistance. 
electrical resistivity, 52 micro-ohm-cm, of M19 
magnetic steel. 
rotor iron volume, meters?. 
maximum rotor radius, 2.0 meters. 
machine shaft revolutions per minute, referred to 
as rotor speed. 


rotor slot space factor. The figure 0.35 was 


used, 
stator copper volume, meters?. 
stator iron volume, meters?. 


guessed slip of the induction machine. 


pal 











patos ath if 
a a at La “) ney | " | 


re 








slipl 
slip2 
smax 
SSF 


T1 
tmaxpl 
viapl 
va 
vagpl 
vol 
VOLALL 


volmag 
vtpl 

Ww 

wr 


wm 


wt 


WTALL 


wt_iron 
wtmag 
x1lpl 
x2pl 
xbeltpl 
xelpl 
xmpl 
xrdpl 
xrspl 
Xs 
xsepl 
xslot 


xsdpl 


larger derived slip. 
smaller derived slip. 
maximum machine slip. 
stator slot space factor. The figure 0.35 was 
used. 

0.014 inch, of M19 magnetic steel. 


per length maximum torque. 


thickness, 


per length Thevénin equivalent voltage. 

'P+jQi, VA rating, volt-amps. 

per length air gap voltage. 

machine envelope volume, meters3, with a margin. 
volume allowance for frame and foundation. A ten 
percent allowance was used. 

volume of permanent magnet material, meters?. 

per length terminal voltage. 

electrical frequency, radians per second. 
width of rotor slots. 

mechanical angular velocity of rotor, radians per 
second. 

weight of copper and iron in a machine, plus a 
margin. 

weight margin for frame and foundation. A ten 
percent margin was used. 

iron weight, kg. 

weight of permanent magnet material, kg. 

per length primary impedance. 

per length secondary impedance. 

per length belt impedance. 

per length Thevévin equivalent impedance. 

per length air gap magnetizing impedance. 

per length rotor differential leakage impedance. 
per length rotor slot leakage impedance. 

p.u. synchronous impedance. 

per length stator end turn impedance. 
slot impedance. 


per length stator differential leakage impedance. 
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xsspl per length stator slot leakage impedance. 


xz2pl per length zig-zag impedance. 


133 





Table 37. Listing of various functions used in computer 


programs 


/$%% a pseudo random nusber generator $3%/ 
tdefine MULTIPLIER 25173 

Hdefine MODULUS 32768 

tdefine INCREMENT 13849 

Hdefine MODFLT 32768.0 

double randoa() 


{ 
extern long int seed: 
seed=(MULTIPLIERSseed+INCREMENT) 2 MODULUS; 
return(seed/MODFLT); 

} 

double swf (n) /$ stator winding factor %/ 
int n; /% haraonic order $/ 
{ 


double kp, kb; 

extern double cos(), sin(); 

kp=cos(0. 31428n); /$ pitch factor, assuaes 0.8 coil pitch %/ 

kb=(sin(0.92368n))/(0.52368n); /t breadth factor, froa 
Kirtley’s "Basic Formulas ... " and assumes an 
electrical winding angle of 40° $/ 

return(kp$kb); 

XY 


F 


double abs (q) 
double q; 
{ 
if (q ¢ 0.0) 
g= (-1.0)8q; 
return(g); 
} 


double sinh(u) 
double u; 

{ 

double exp{), ans; 


ans = 0,98fexp(u) - exp(-u))s 


return(ans) ; 
} 


double cosh(u) 
double u; 
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{ 
double exp{), ans; 


ans = Q.ob(exp(u) + exp(-u)); 
return(ans); 


} 

double rwf(n) /t rotor winding factor, same comments as swfi) $/ 
int n; 
{ 


double kp, kb; 

extern double sin(); 

kp=1; 
kb=(5in(0.92368n))/(0.52368n); 
return(kp$kb)s 


4 
4 


#define tol 1.0e-4 
Fdefine pi 3.141592654 


double besi (p,x) /‘t bessel Iforder, arg) $/ 
double x; 
int p; 
{ 
double bi; 
int bial); 
double exp(),sqrt(),bis(); 
if (p<0) abort (" besi: negative index"); 
if {x<0) abort (" besi: negative arguaent"); 
if (x>60) bDi=exp(x)/sqrt{2%pisx): 
1f((xCd2)iibialp,x,&bid) bi=bis(p,x); 
return (bi); 


ul 
4 


double bis(p,x) 
double x; 
int p; 
{ 
double fabs({); 
double xx; 
int 1,¢k,k; 
double b1i=05 
double t=1; 
xx=x/2.5 
for (i=lyi<ptigitt) t=t8xx/1; 
if (t¢1.0e-34) bi=05 
else 
{ 
bi=t; 
XX=XXEXX5 
for (k=l (kCLOO1) &&( (fabs (t)-fabs (bi Stol)) 0); k++) 
{ 
fk=k$ (otk): 





tatdxx/tks 
bi=bitts 
} 


\ 
F) 


return (bi); 
} 


int bialp,x,pbi) 
double x,%pbi; 
int p; 
{ 
int ist, fn,k, fk; 
double xx,bi,t; 
double fabs(),sqrt() ,exp(); 
fn=4iptps 
t=1.0; 
bi=1.0; 
xx=0,125/x3 
for (k=Ls (k<30)&&( (fabs (t)-fabs(bistol)) 0) sk++) 
{ 
#k=(20k-1) 8 (2tk-1)3 
tstdxxd(fk-fn) /k: 
bi=bitts 
} 
if (k==31) ist=1; 
else 
{ 
ist=0; 
bi=bisexp(x)/sqrt(2,08pix): 
} 
ipbi=bi; 
return (ist); 
} 


double besk(p,x) /$ aodified bessel function Kp(x) $/ 
double x; 
int p} 
{ 
double bk; 
double exp(),sqrt(),kO0) ,ki(): 
if (p¢Q) abort (" besk: negative index °); 
if (x¢0) abort (" besk: negative arguaent “); 
if (x>60) bk=exp(-x)/sart(2.08x/pi);: 
else if (p==0) bk=kO(x); 
else if (p==1) bk=kilx); 
else { 
double g0,q1,9); 
int J; 
gO=kO(x)s 
glkilx); 
for (j=23j)<ptisjtt) 
{ 
Qj=28(j-1) 8gi/x+q0; 
g0=g1; 





gl=qi; 


bk=9); 
} 
return (bk); 
} 


double kO(x) 
double x; 
{ 
double bk; 
double log(),sqrt(),exp(); 


if (x<1) 

{ 

double a,b,2,¢,90,x2j.f,hj,rj3 

int 3; 

b=, 98x; 

a=.97721366+t1 0g (b) ; 

c=b8b; 

g0=-a; 

x2)=13 

f=15 

hj=0; 

for (j=13j«73 j++) 
{ 
rj=1.0/); 
x2)=x2)%c3 
f=f$rjir): 
hjahjtrj3 
g0=g0+x2jst8(hj-als 
) 

bk=90; 

} 

else 

{ 

double t{12),a,b,c,pa; 

int 1; 

a=exp(-x)s 

b=1.0/xs 

c=sqrt(b); 

tCOJ=b; 


for (L=1slCi2she+) t£1J=tl1-1)8b; 


pa = 1.2555141-.15666428t0); 

pa t= .0BB111288t(1)-.091390958t12); 

pa t= .15445968t(3)-.22998508t(4); 

pa t= .$792418t(5)-.52472778t04); 

Pa t= .99/9S6BSt(7)~-.426263358t18); 

pa t= ,21B45188t(9)-.066809778tL10]+, OO91B9SBSst LiL); 


bk=asctpa: 
} 
return (bk): 


\ 
4 
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double ki (x) 


double x; 


{ 


double bk; 
double log{),exp(),sqrt(); 


if (x<{) 


else 


{ 

double a,b,c,gi,x2j,f,hJj,rj; 
int J; 

b=x/2: 
a=.9/7721966+1 09 (b) ; 
c=b3b; 

x2J=bs 

f={3 

hy=t; 
Ql=t.0/x+x2j3(.Sta-hj)s 
for (j=23j<95 j++) 


{ 
x2)=x2)8c3 
rj=t.0/)s 
f=f3rjtr);: 
hj=hjtr); 
gl=qitx2j3f3(.5+la-hj)3j); 
} 

bk=q1; 

} 

{ 

double a,b,c,tl!2),pa; 

int 1; 

a=exp(-x)s 

b=1.0/x; 

c=sqrt(b); 

t{OJ=b; 


for (l=ts1Ci2s]4+) tlPJ=tll-£18b; 


pa = 1.2533141+.4499978t(0]-. 1468583st LL]: 
pa += .12804278¢(2]-.17344328t(3)+. 284761984]; 
pa += -.45943428t(5]+. 6283381 stL]-. 66322958 [7]; 
pa += .50502398t(8]-. 25813048t(91: 

pa t= .O7BRQ00LStL10] -.010824188tL 11]: 


bk=asc3pa; 
\ 


return (bk); 


} 


double besip (p,arg) /$ deriv of bes I(order, arg) %/ 


int ps 


double arg; 


{ 


double x,y,2; 
x = besi (p-1,arg)s 


138 





y =p t besi (p,arg)/arg; 
c= xk Vy 

return (2); 

} 


double beskp(p, arg) /t deriv of bes K(order, arg) %/ 
int p; 

double arg; 

{ 

double x,y,z3 

x = -besk (p-l,arg); 

y=- p t besk (p,arg)/arg; 

eek * YS 

return (z): 
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Appendix B. Synchronous Machines and General Relations 


An equivalent circuit for a synchronous machine is 


Kaf 





Figure 35. Synchronous machine equivalent circuit 


where Ra is the stator resistance and Xs is the synchronous 
reactance. The internal voltage of the machine, Eaf, is 
developed between the stator and rotor across the air gap. 


The current direction is shown as if the machine is a motor. 


Vta = Rala + jXsIa + Eaf 


Eaf represents a mutual coupling between the stator and 


rotor, and 


wM Ie 4 Yo ki ke 1 r Ns Nf 


If is field current, Nf is the number of series field turns 
and ke is the fundamental rotor winding factor. Since it is 
never desired to specify the number of turns on either the 
stator or the rotor, a scheme has been devised so that 
derivations are conducted in volts-per-turn, ampere-turns, 


and ohms-per-turn-squared, which results in power in watts. 
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Figure 36. Phase belt conductor area 


For the stator, the transverse area can be divided into 
phase belts. The area labelled ‘A’ is the area occupied by 


e 


phase ‘a’, in one direction. Then, Ns Ia = Aa Ja, where Ja 


& ’ 


is the current density in phase ‘a’ conductors. The area 
has some effective conductor area, subject to the need for 
conductor cooling passages and insulation area. Therefore, 


for stator and rotor currents, analysis yields 


fie and it e= ==> == 


It is also desired never to specify the number of slots on 
either the stator and rotor. Accordingly, slot space fac- 


tors are defined as 
slot 
| width | 


| slot pitch | 


Figure 37. Slot space factors 
stator slot width rotor slot width 
~ stator slot pitch rotor slot pitch 





For a typical turbogenerator conductor bar, the copper area 
is about thirty-five percent of the conductor envelope area. 
Variables titled SSF (stator slot space factor) and RSF 
(rotor slot space factor) embody this thirty-five percent. 


The conductor area of a single stator phase in one direction 


(r+g) em ls ds SSF 
Rens cg Se ee and, for the rotor, 


Af = ---------------- There is only one phase on the 
2 rotor. 


where rho is the electrical resistivity of copper and lt is 


the turn length. If a stator turn can be modeled as 






(cincuMFEREN cE) cP 


P 





| ACTIVE LENGTH | 


Figure 38. Stator turn 


and the circumfrence is en(rtg), then 


where J(4/3) is (1/sin(60°)) and CP is the coil pitch. Then 


Ra 2 rho { 1 + (2an(rt+g)CP/p)V¥ (4/3) } 


Ns2 wm (r+g) ls ds SSF 
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The stator and rotor copper losses are 
Rala* = rho lt Aa Ja? RfI*¢* = rho lite At Je? 
The rotor is full pitched. 


The vector diagram for the synchronous machine equiv- 


alent circuit is (assuming Ra ==> QO) 





Figure 39. Synchronous machine vector diagram 


By use of the law of cosines, 


EFaf? = Vt2 + (XsIa)? - 2VtXslIa cos(r + QO? ) 


where t is the power factor angle. If this equation is put 
into per unit form, with eaf = Eaf/Vt and xs = XsIa/Vt, then 


eaf* = 1+ xs* + 2 xs sin(T) 


xS can be calculated, allowing the calculation of eaf. 

There is simple linear relationship between eaf, the no-load 
rotor current density (jrnl), and the full-load rotor cur- 
rent density (jr). It is jr = (eaf)(jrnl), because Eaf is 
directly proportional to It and It is directly proportional 
to Jé. Ampere’s Law states 


p Hr « dl = | J « n dA 


If one chooses an integration path around half of the rotor, 
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If the no-load condition is chosen, and it is recognized 
that Br is essentially constant with increases in Vt 


(because when Vt increases, so must EFaf), then 


2Ho r lr dr RSF ke jrnl 


and jrnl = --------------------- 
2o r lr dr RSF ke 


The synchronous reactance is a measure of flux linkage 
within the machine. It is Xs = w Ls = w (Lal + 1.5 Laa), 
where Laa is the armature single phase inductance and Lal is 
the slot leakage inductance. The factor of 1.5 is derived 
from the 120° separation between the three phases of these 


machines, as 


XsIa = w (IaLal t+tlaLaa + IpnLab + IcbLac) 


Due to the symmetry of the machine, 


Lab = Lba = Lac = Lea = Lbec = Leb = -0.5 Laa 


Then, XsIla = w la (Lal + 1.5 Laa). 


Self-inductance is 


4 Ho ks? Ns* lr 


from Ampere’s Law and L = (flux/current). If a single con- 
ductor per slot is postulated and the effects of slot teeth 


are ignored, then the leakage inductance is 


Jb LI Sa a (Pself + Pmutual) Ns? 
(pole) (phase) 
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where P is reluctance. For this conductor configuration, 


Ho ] Vo 
Pself = a ds ws and Pmutual = aa ds ws, 


with ws equal to slot width. Since ws is not known, we use 
the stator slot space factor, ls, multiplied by the number 
of slots to yield 


(slots) ws = ls tw (rtg) Therefore, 


5 Ns? H/o 1 ds is w (rtg) 


The real power developed by the machine is 


Pwr = 3 Vt Ia pf, where pf is the power factor. Through the 
use of Lenz’ Law and Ampere’s Law, terminal voltage may be 


expressed as 


Using our previous relation for Ia, the expression for Pwr 


is 
on r lw Br ks Ja SSF (rtg) 1s ds pf 


Finally, winding factors need to be derived. The wind- 
ing breadth factor, kb, is 


sin(m nP/2) 
m sin(n [/2) 


where f is the electrical angle between adjacent slots, n is 
the harmonic order, and m is the number of slots per pole 


per phase. The winding pitch factor, kp, is 


145 





kpn = Sin (n «/2) where « is the electrical angle 


between sides of the coils (pitch angle). For a three phase 
winding, x = 2m p. If a0O.8 coil pitch is assumed (which 
will rid the machine of certain harmonics during balanced 
operation) then « = 180°(1 - 0.8) = 36° and 

kpn = sin(0.3142 n). ASSUmptions are needed to calculate 
kon without specifying the number of stator slots or turns. 
If the winding angle is specified as did Kirtley [32], then 


sin(n @w/2) 


A reasonable electrical winding angle is 60°, since most of 
the stator periphery will contain turns. The breadth factor 
devolves to kbn = sin(0.5236 n)/(0.5236 n), for which the 
fundamental harmonic factor equals 0.955. The winding fac- 
tor is the product of the breadth and pitch factors. For 


the rotor, a pitch factor of one is assumed. 
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Table 38. Listing of synchronous machine design program 


finclude "stdio.h* 
tinclude “def.h° 


/$ prograa nage: syn.c for synchronous, round rotor aachines %/ 


lang int seed; /% start point for random nuaber generator %/ 
double b£263£26), hCi1](26), ks{8), kr{8); 
/tb is "best" array, his “hold” array, ks/kr are winding factors %/ 
agin() 
{ 
double design point(), rnd_walk(), swt(), rwf(), ke, kv, ainper, 
stepsize, randoa(), abs(), freq, rpa; 
int p, iteration, 1, }, best, print_out(), loops; 
FILE %fopen(), &fp; 


printf (*\nReading input data froa SYN.DAT .. ."); 


fp=fopen("syn.dat","r"); /% input seed for rando@ nuabers %/ 
fscanf(fp, "4d", &seed)s 


fscanfifp, "td", &p); /t input nuaber of pole pairs $%/ 

fscant(fp, “Z1¢*, &minpwr); /¥ input machine power, derived fa ASSET %/ 
pinpwr $=746.0; /t convert to watts %/ 

fscantitp, “Z1f*, &kel; /t CERs for Effective Weight %/ 

fscant(tp, "Zt", kkv); 

fscantlfp, "214", &rpal; /t machine max shaft rpe %/ 

fclose(fp); 


printf ("\nHow many loops do you want? “); 
scant("Zd", &loops); 


printf (*\n\nDoing program calculations . . .\n"); 


for (i=ls i ¢ Bs 1+=2) /t harmonic winding factors %/ 
( 
ksCi J=swt (i); 
krCilj=rmf (ids 
} 


freg=rpatp/40.0: /t wax electrical frequency %/ 
/t MAIN BODY OF THE PROGRAM $/ 


for (i=l; 1 (= loops; ++i) 
( 
stepsize=0.1; 
iteration=0; 
design point(ainpwr, p, ke, kv, freq); 
/t put stuff in the hold array %/ 
while (iteration <= 10) 
{ 
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rnd walk(minpewr, p, stepsize, ke, kv, freg); 

/t stagger around %/ 
best=0: /t index to best EW of the lot $%/ 
for (j=l j<=10s +3) 

if (h0j)018) ¢ h{bestj(18)) 
best=j: /$ find the best sachine %/ 

if (abs((h(O3018) - hlbestI£183)/h(O3018)) < 0.005) 

/t smal] improvement in EW &/ 


f 
5 


stepsize/=2.0: 
+titeration: 
} 
else /t transfers best to 0 position %/ 
{ 
for (j=l: j <= 253 ++#)) 
hCQICj) = h€bestj(jd; 


} 
} 
for (j=l; J <= 203 ++) 
b(1)0j)J=hCbestilj]; /t keep the best sachine %/ 


1 
d 


best=1; 
for (1=1; 1 <= loops; +1) 
if (b0i]018) ¢ b€best3(i8)) 
best=i: /t find and keep the best of the best &/ 
ainpwr/=746,0; /% turn back into hp %/ 


print_out(best, p, ainpwr, ke, kv, rpm); /8 output to disk file %/ 


fp=fopen(“syn.dat", “w"); /t output seed %/ 
forintf(fp, "ad", seed): 

fprintf(lfp, “\ndd", p)s 

fprintflfp, "\ndlf", ainper); 

fprintflfp,"\ndlf", ke); 

fprintflfp, “\ndlf", ky): 

fprintflfp, "\ndlf", rpa)s 

fcloselfp); 

} 

/t END OF MAIN PROGRAM; ALL THAT FOLLOW ARE FUNCTIONS %/ 


double design point(minpwr, p, ke, kv, freq) 
/t determines a random design point %/ 
double sinpwr, ke, ky, freq; 
int p; 
{ 
double r, jrnl, jr, js, 15, Ir, dcore, ds, dr, g, w, 1, x5, eaf, i2rr, 
va, ph, pe, i2r, vol, wt, effcy, ew, xsi, xsd, x57, xsal, 
cW, SiV¥, Scv, riv, rcv, loa, doa, wt_iron, find siv(); 
extern double sgrt(), randoal); 
int c=0, d=0; 


while (d ‘= 1) { 
while (c '= 4) 
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{ 
r=randoe() $kKMAX; 
if ((28P srt treq/p) < MAX_TIP_SPEED) /% check rotor tip speed %/ 


breaks 
} 
w=29P]8freq; 
lr=randoa()80.5 + 0.29; /t rotor slot factor %/ 
dr=randoa()%r/5.0; /% slot no deeper than 20% of rotor radius %/ 
dcore=(BR&r)/(RSAT&p) /t back iron depth %/ 
ds=randoa() $0, 9%dcore; /t slot depth ¢ 902% of body depth %/ 
while (c '= 1) /% gap dimension %/ 
{ 
gsrandom()3(0.1%r - GMIN) + GRIN; 
if (g > 0) 
break: 
} 
ls=randoa(}80.5 + 0.29; /t stator slot factor %/ 
js=randoa() tJSMAX: /t full load stator current density %/ 


jrnl=(BRegsp)/(BSMUSRSFSrsdrslrskr(1)); /& noload rotor current density %/ 


xsl=ks(1)8ks(1); 
xso=(ks{9)8ks{5)/25); 
x57=(kS(7)tks(7)/49): 
xSal=(58lstPl$ds/id); 
XS=(MUSSESSFR(r+g)Blr+giedstPlsls & (128(xsi + xs5 + x57) /(PI8gtp) + xsal)) 
/(12er tBRtks(1));5 
/% p.u. synch iapedance %/ 

if (xs > 2.0) 

continue; /t don’t want xs too big %/ 
else 

this 


eaf=sqrt(1 + xstxs + 28x580.6)5 /% 0.4 is sin(t), pwr factor angle, 
eaf iS p.u. internal voltage at full load %/ 


jr=eattjrnl; /% jr full load, linear with eat %/ 


]=(minpwrtp)/(28BRewsksll sr sSSFedssPls(rtg)8lstjstPF); /% active length %/ 


va=Z2bPlordl sw8BReks(1 8 jstSSFs (r+q) sdstls/p: /% va rating %/ 
Siv = find siv(l, r, q, ds, dcore, 1s); /% stator iron voluae %/ 
riv = 18Pléré(r - 28drélr);: /$ rotor iron voluae %/ 


SCV = 28P[8(rtg)tdst]ss(] + 2. 30948P 18 (r+q) SCP/p): 
/t stator copper volume %/ 


rey = 28Plérsdrdlrd(] + 2, 50948P]8r/p); /t rotor copper voluae %/ 
cw = (rcv + scv) SDCU; /% total copper weight %/ 
loa = 1] + 48(rtq); /t length-over-all %/ 

doa = 28(r+gtds+dcore); /t over-all-diaaeter %/ 
vol = YOLALLS(loasPltdoatdoa/4); /t machine envelope voluae %/ 

wt = WIALLS(cw + DS(BRNGSSriv + siv)): /% aachine weight in kg %/ 
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wt iron = De(riv + siv); /% iron weight only &/ 


ph = 31. 86225BETAsfreqsBRIGHCléwt iron/D; 
/t hysteresis loss in watts, uses iron weight of sachine %/ 
pe = (106236. 9dNUSBSATEBSATIT IST freqtfregtwt_ iron) /(RHOMD); 


/t eddy current loss in watts, uses iron weight of aachine %/ 


i2r = 2,08SSFtdssCRHOsPLs (1 + 2.30948P Lb irtg)&CP/p) tjstjst(r+g) tls: 
/% stator copper loss in watts §/ 
/t revised 1-12-87 %/ 


i2rr = 2,00RSF&drértlretjrtjrtCRHOSPls (1 + 2, 30948PL&r/p): 
/t rotor excitation losses, 1-12-87 %/ 
i2rt=12rr; /t total copper losses %/ 


effcy=(sinpwr)/(minpwr + ph + pe + 12Fr)3 
ew=wt + keb(l-effcy) + kvavol; /&t Effective weight ¢/ 


h{OICiI=js; h{O}(2)=freq; hlOIL3I=0; h(OIC4)=r3 hCOILSIJ=93 

hCOM{6]=dcore; hlOIC7]=ds; h{OIC8I=dr; HCOILPI=ls; HCOILLOI=Ir; 

h{OMC11J=vol;  hlOILI2I=wts hCO}013)=ph; hCOIL14J=pe; h{OIE1S)=12r; 

h{OIC16)=va; HCOIC7I=effcy; hlOIC18)=en; hfOICI9I=1; h{OIL20)=jrs 

hCOME21=jrnl;  hlOIC2ZI=xs5 h{O}C23]=eat; hCO1L24)=loa: h{OQI(25)=doa; 
/t this section just changed al] the variables in the “hold® array %/ 


return; 
}} 


double rnd_walk(sinpwr, p, stepsize, ke, kv, freq) 
/t walks about design point 10 times %/ 
double stepsize, ainpwr, ke, kv, freq: 
int p: 
{ 
double r, jrnl, jr, js, ls, Ir, dcore, ds, dr, gq, w, 1, xs, eaf, i2rr, 
va, ph, pe, i2r, vol, wt, effcy, ew, xsi, xso, x57, xsal, 
Cw, Si¥, SCY, riv, rcv, loa, doa, wt_iron, find_siv(); 
extern double sqrt(), randoa(); 
int i=]; 


while (i <= 10) 
{ 
/t read in the walk around the design point %/ 


js=h(O]C1]8(1 + stepsizes(randoa() - 0.9)); 
if (js > JSMAX) 
js = JSMAX: /t reset to liait ¢/ 
w=20PT8freq; 
r=nCO](4]0(1 + stepsizet(randoa() ~ 0.5)); 
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if ((wir/p) > MAX_TIP_SPEED) 

continues /% go to next try if violated %/ 
g=hlO}LS]8(1 + stepsizeslrandoal() - 0.5)); 

if (g < GHIN) 


g=GNIN; /t reset to the liait %/ 
dcore=(BRtr) /(BSATSp); /t aost efficient use of iron %/ 
ds=hlOIC7I8(1 + stepsizesl(randoa() - 0.5)); 
if (ds > dcore) /t can’t have too-deep slots %/ 
ds=dcore; /t reset to the liait $%/ 


dr=hCO}(GI8(1 + stepsizet(random() - 0.5)); 
Is=hCOJL9IS(1 + stepsizes(random() - 0.5)); 
ifetls)> 0.70) 
1s=0.753 /t reset to the liait %/ 
if (ls ¢ 0,29) 
1s=0. 2953 
Ir=hCOJLLO]8(1 + stepsizes(randoa() - 0.5)); 
if (Ir > 0.79) 
Ir=0.793 /‘¥ reset to the liait %/ 
rf lin 0220) 
1r=0.253 


/% computation section of the walk $/ 
jrnl=(BR&qgtp) /(88MUSRSF Sr SdrSlr$kr£01)); /$ noload rotor current density $/ 


xsl=ks(1]$ks{1); 
xso=(ks{S)8ks5{5]/25) s 
xSJ=(ks{7)$ks{7)/49): 
xsal=(S8lssPi$ds/18); 
XS=(MUSJSSSSFS(rtg) Slrtg)$ds$PI$ls $ (128lxsl + x55 + xs7)/(PI$g%p) + xsal)) 
/(i2srsBRsks{i)); 
/% p.u. synch impedance $/ 
if (xs > 2.0) 
continue; /$ can’t have xs too big $/ 


eaf=sqrt(1 + xs$xs + 28x580.6); /8 0.4 is sin(t), pwr factor angle, 
eaf 15 p.u. internal voltage at full load ¢/ 


jr=eaftjrnl: /% jr full load, linear with eaf $/ 
l=(minpar$p)/(28BRewiksCi]srsSSFsdstPis(r+g) &lstjssPF); /% active length $/ 
va=Z8PT8r bl &wSBRoksl1 IS jstSSF 3 (r+q) Sdstls/p: /¥ va rating $/ 


siv = find siv(l, r, gq, ds, dcore, 1s); /% stator iron voluae $/ 


riv = 18Pl8rt(r - 28dr8lr); /t rotor iron voluae $/ 
SCV = 28P]8(rtg) tdstlss(] + 2. 30948PI8 (r+g)8CP/p): 

/% stator copper voluae %/ 
rev = 29Pl$rddrhlrs(l + 2.50948PL8r/p); /$ rotor copper voluae $/ 
cw = (rcv + scv) S8DCU; /$ total copper weight $/ 
loa = ] + 48(r+q); /$ length-over-all 8&/ 
doa = 28 (rt+gtdstdcore) ; /t over-all-diaaeter $/ 
vol = VOLALL& (loatPI$doasdoa/4); /$ aachine envelope voluae %/ 
mt = WIALLE(cw + DS(BRNGStriv + siv)); /t machine weight in kg $/ 





wt iron = D8(riv + siv); /t iron weight only %/ 


ph = SL. Bd2253BETASfreqtBRIsHCi swt iron/D; 
/t hysteresis loss in watts, uses iron weight of aachine %/ 


(106236. 9INUSESATSBSATIT 1ST is fregtfreqdwt_iron)/(RHO8D); 
/t eddy current loss in watts, uses Iron weight of @achine %/ 


pe 


i2gr = 2, 08SSF&ds#CRHOSPls (1 + 2.3094SP 18 (r+g)8CP/p) sjstjstirt+g) tls: 
/¢ stator copper loss in watts %/ 
/t revised 1-12-87 %/ 


l2rr = 2.08RSFadrtrtlrsjrtjrtCRHOsPrs (1 + 2.30943PI8r/p)s 
/t rotor excitation losses, 1-12-87 %/ 
12rt=i2rrs /t total copper losses %/ 


effcy=(ginpwr)/(ginpwr + ph + pe + i2r); 
ew=wt + ked(l-effcy) + kvévol; /t Effective weight %/ 


h{i}01)=)s; hlill{2)=freq; hlill3)=e0; h{ijl4)=r; hlil{S)=9; 

h{il($)=dcore; Alil{7)=ds; hi }(8)=dr; hCi3(9)=1s; h{i}(10J=Ir; 

hCrWCiiJ=vols = ALi NC12)=nts h{iJ013)=phs h{i JC14)=pe; hi M15 )=12rs 

h{iJ(16)=va; hCidli7J=effcy; hli)(18)=en; ACiICiA9I]=1;5 hCr)(20)=)jr; 

hCid(2t)=jrnls = hlid£{22)=xs; hCil(23)=eats = hlid{24)=loa; hlill{25)=doa; 
/t this section just changed all] the variables in the “hold® array ¢/ 


ths ‘t go to the next hilt) ¢/ 
} 


return; 
} 


print_out(best, p, @inpwr, ke, kv, rpa) 
int best, p; 
double ainpwr, ke, kv, rpa; 

f 

char outfilet14); 

FILE %fpo, tfopen(); 

mt 1} 


printf ("\n#hat is the nage of the file where you want the output? "); 
scanf("ds", outfile); 


fpo=fopen(outfile, “w")s 


fprintf(tpo,"Zd", p); 
forintf(fpo,“\ntlt", @inpwr)s 
fprintf(fpo,"\ndlf", ke)s 
fprintf(fpo, "\ntlt", kv)s 
forinti(fpo,"\nzlf", rpa)s 





for (1=f3 1 ¢= 20; ++) 
fprintf(fpo, "\ndlf*,blbestjlid); 
fprintf(fpo, “\n"); 
fcloselfpo); 
} 


double find sivll, r, g, ds, dcore, 1s) 
double 1, r, g, ds, dcore, Is: 

{ 

double one, two, three, four; 


one = (rtgtdstdcore)tirtgtdstdcore) - (r+q)%irtq)s; 
two = 28Pitirtg)tdstis; 
three = (r+qtdstdcore)tirtgtdstdcore) - (rtgtds)t(rtqtds): 


four = 13(Pitone - two) + PIs4tir+g) three; 
return(four)s 
} 





Table 39. Listing of synchronous efficiency program 


#include “*stdio.h" 
firiclude “def.h° 


/$ program name: seft.c to find efficiency of synchronous aachines %/ 
/t works with only a single aachine %/ 


aain() 

{ 

FILE Sfopen(), %fp; 

Goubher, irnl, jr, JS, 15, Ir, dcore, d5, dr, g, w, |, x5, eat, i2rr, 
va, ph, pe, izr, vol, wt, effcy, ew, xsl, xs0, x57, xsal, ks(8), 
Siv, riv, wt_iron, find siv(), paett, pajs, parpa, dhp, rpa, @inper, 
ke, kv, freq; 

extern double swt(), sqrt); 

int e=0, #, p, 13 

char intile(t4); 


POpeiy—13 1 « 8; it=2) /t haraonic winding factors %/ 
{ 
ksCi J=swt (i); 
} 


xslaks(h)8ksCl); 

K§9=(k509)8k5(95)/25); 
x57=(ke07)8k5(71/49); 
aSal=(StlssPl$ds/18); 


printfi*\nCalculates efficiency of a single sotor.\n"); 
while (e '= 1) 

{ 

f = 0; 

orintf(°What 15 the nage of the input file? °); 
scanf("ds", intile); 

fp = fopen(infile, “r°)s 


fscant(fp, "Ad", &p)s /‘t input nuaber of pole pairs %/ 
fscantifp, “Alf", &einpwr)s 
ainpwr t= 746.0; /$ now in watts %/ 


fscanfitp, “Zlt", &ke); 
fscantitp, "Alt", &kv); 
fscanf(fp, "Alf", &rpa); 
fscantitp, “Zlf", &js)s 
fscantltp, "Alf", &freq); 
fscantitp, “Alf, &w); 
fscantifp, "Alf", &r); 
fscanfifp, "Alf", &q); 
fscantltp, "Alf", &dcore); 
fscantitp, “Zlf*, &ds); 
fscantitp, “4lf*, &dr): 
fscantifp, "Zl", &ls): 
fscantifp, "Alf", &lr); 
fscantifp, “Al#*, kvol); 
fscanfltp, "Alf", &wt); 





fscantifp, “Z1f", &ph); 
fscanflfp, “Zlf", &ped; 
fscanf(fp, “ZLlf", &izr)s 
fscantitp, "Alf", &va)s 
fscanfifp, “Zlf", keffcy); 
fscanf(fp, “Zlf*, kew): 
fscant(fp, "£1", &l); 
fscanfifp, “£1#", &jrds 
fscant(fp, “Z1f", &yrnl); 


fclose(fp); 

siv = find siv(l, r, g, ds, dcore, 1s); /¥ stator iron volume %/ 
riv = L8Plard(r - 2kdrtlr)s /¥ rotor iron volume %/ 
wt iron = Dt{riv + siv); /¥ iron weight only %/ 

while (f '= 1) 

{ 


printf("\ndhat is the sustained speed machine horsepower? *); 
scant(*Zlf", &dhp)s 

dhp %= 744.0; 

printf(*What is the sustained speed aachine rpa? "); 
scant("Zlf", &parpa)s 


freq = parpatp/40.0; /¥ max electrical frequency %/ 

pajs = jstdhp/ainpwr; /t FA stator current %/ 

KS=(MUNSSSSFR(rtg)dlrtg) SdssPlals § (128( xsl + x55 + xs7)/(PI8gtp) + xsal)) 
/(12arsBRaks{1)); /t p.u. synch impedance %/ 


eaf=sqrt(l + xshxs + 28xst0.6); /$ 0.6 15 sin(t), pwr factor angle, 

eaf iS p.u. internal voltage at full load %/ 
jr=eaftirnl; /t jr full load, linear with eaf %/ 
ph = 31. 862259BETAMfreqtBRisHCi§wt iron/D; 

/t hysteresis loss in watts, uses iron weight of machine %/ 
pe = (106236. 9SNUSBSATSBSATS T1811 ttreqtfreqtwt_iron) /(RHO8D); 

/¥ eddy current loss in watts, uses iron weight of sachine %/ 
12r = Z.Q855FedstCRHOSPl$ (1 + 2.50949PI8(r+g)8CP/p) Sjstistirtg) sls: 

/¥ stator copper loss in watts %/ 
i2rr = 2.08RSFadrtr sledge sir sCRHOMPls (1 + 2. 30949P14r/p): 
/t rotor excitation losses, 1-12-87 %/ 

paeff = dhp/(dhp + ph + pe + i2r + i2rr); 
printf("\n Sustained speed efficiency is Z1f", paeff): 


printf("\nWhat 1s the endurance speed machine horsepower? *); 
scanf("Zlf", &dhp); 

dhp t= 744.05 

printf("What is the endurance speed machine rpa? "): 
scanf("Zlf", &parpa): 


pajs = jstdhp/ainpwr;: /¥ PM stator current %/ 

freq = parpatp/40.0: /4 max electrical frequency %/ 

KS=(MUR}SESSES (rtg) a (rtg)adssPI8ls & (128(xs) + x55 + x57) /(PI8gtp) + xsal)) 
/(129r SBRSksC1));5 /t p.u. synch iapedance %/ 


eaf=sqrt(l + xsixs + 28x580.6)s /8 0.6 is sin{t), pwr factor angle, 
eaf is p.u. internal voltage at full load %/ 
jr=eaftjrnl; /t jr full load, linear with eaf %/ 
ph = 31, 852258BETASfreqtBR1sHCiswt_iron/D: 
pe = (106236. 98NUSBSATSBSATSTIsTisfreqtfreqdwt_iron) /(RHO8D); 


ara! 





i2r = 2, 08SSFtdstCRHOsPls (1 + 2.50948PIs(r+g)sCP/p) tystistirtg) tls; 
b2rr = 2.08RSF&drtr slr tyrtyrtCRHOsPlt (1 + 2.50948PL tr /p); 

paeff = dhp/(dhp + ph + pe + i2r + i2rr); 

printf("\n  €ndurance efficiency is Z1f*, paeff); 


printf ("\nSame aachine? °); 
scanfi("id", &f); 


if (f == 9) 
continue; 
else if (f == 2) 
{ 
e=i; 
break; 
} 
} /t end of f-loop %/ 
) /t end of e-loop %t/ 
} /¥ end of @ain prograe %/ 


double find siv(], r, g, ds, dcore, 1s) 
double 1, r, g, dS, dcore, 15; 

{ 

double one, two, three, four; 


one = (rtgtdstdcore)&(rtgtdstdcore) - (r+q)8(r+q); 
two = 28P18(r+g) tds8ls; 
three = (rtgtdstdcore)&(r+gtdstdcore) - (r+gtds)t(r+gtds); 


four = 1%(PItone - two) + PIt4s(r+g)tthree; 
return(four); 


} 
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Appendix C. Permanent magnet machine 


The equivalent circuit for a permanent magnet machine 
is almost identical to that of a synchronous machine. The 
only difference is the source of the internal voltage, which 
develops the field flux wave that interacts with the arma- 
ture flux wave. The field flux wave is a result of per- 
manent magnets built into the rotor to develop magnetic 
poles. 


A typical magnetic circuit, combined with Ampere’s Law, 


shows 


Hmlm + Hee = Ni 





Figure 40. Typical magnetic circuit 


We use a constitutive law of B = !H and assume that any 
steel has »P =o, If flux leakage is about zero, BmAm = BeAg 


Since flux is solenoidal. Then, 


Minimum magnet volume occurs when the magnet’s maximum 


energy product (MEP), HmBm, is a maximum. If current is 
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zero and only magnitudes are used, 


The load line of the magnet is developed as 


lm g Bm 

i a tases a and 
Am Ho Hm Ag 

Bm lm Ag 
aos Se NO SHS See ee =e 

Hm g Am 


The permanance coefficient, Pc, is the slope of the load 


line. On a magnet diagram 


iE es 


load line 


He 


Figure 41. Magnet operating point diagram 


A good algorithm for machine design is to adjust the 
geometric dimensions to operate at the MEP, on the load 
line. If operation at MEP is assumed, the needed slope is 
determined, the dimensions are randomly generated, PoHg is 
calculated, and the design is maximized for VPoHge and mini- 
mized for magnet volume, then a search technique has been 
delineated. 

Magnets may “overhang” the active length at either end 
to account for manufacturing imprecision and to permit a 
smaller armature diameter. This overhang affects the 


developed flux. 
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Figure 42. Permanent magnet machine diagram 


The flux per pole, @, with overhang and a given Ja and Jn, 
is the same as would exist in a configuration in which 

Ja’ = Jm’ = Ja + ND and there were no overhang. N iS a non- 
linear parameter promulgated in graph form by Ireland [33]. 


A good power fit for N is 


5 Jm 
N = 0.38558 Se a 


5 Jm 
where 5h" = overhang and 0O £ overhang £ 0.34. Then, 
Ja § Jm 
D 
ee and Owith = Ow/o Y 

Ja 
an © N 
D 


This flux-with-overhang is applied to the problem as would 


be the usual flux. What is the usual flux? A permanent 
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magnet hysteresis diagram shows the residual flux density to 


be defined at a single point. B 


aM Why 


Figure 43. Permanent magnet hysteresis diagram 


The magnet does not operate at that point, but rather 
on the load line. Also, there is not magnet material at 
every point along the circumference of the rotor. A Fourier 
series is a good way to find the flux-without-overhang. 


Using a developed rotor, 





Figure 44. Magnet material on a developed rotor 


MyM = Bd 
wr p wr 
lr so ee ee K|—e TF Re ee 
wr/p T xr 
= nw x 2nr 
B(x) = £ An cos(------- ae | Da eau 
n=1 4p p 
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2 T ni x 
An = --- Sec ecos.| .——_—— . Ge tOren = 1, 2, 3 
ii 0 I 


Bex) = Bd HOR Sox & (tr ir/p) 
0 for (wr lr/p) i x & (tr/p) 
Bd fOorm(tr,p) Sx S&S (7xr/p) (itilr) 
0 for (wr/p)(1+lr) £ x £ (2mr/p) 


After integration and manipulation, 


The equation represents only the fundamental term of the 


Fourier series. The flux-without-overhang is 
2Bd vw lr vw Ilr xp 
Br = ----- [ sin(--7--) = SESS) iecos(---—) 
TT 


As is usual in steady-state analysis, the magnitude is used. 
The next quantity to find is Bd. If the magnet 


operates on the load line, the operating point flux is 


| BC i tte a ego where ¥r is the relative 
1 + Hr /Pc 
reversible permeability of the magnet and Brem and Pc are as 
previously defined. With this relative permeability, the 
magnet length is lm = g Pr. The magnetic machine can now be 
specified. An end view with dimensions is given in Chapter 


Four. 
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Table 40. Listing of permanent magnet machine design program 


finclude “stdio.h’ 
tinclude "def.h" 


/$ prograa naee: paa.c for permanent aagnet aachines %/ 


long int seed; /t start point for random nuaber generator %/ 
double BL101)£27), h{i13£27), ks€8); 
/¥ bis "best" array, his "hold* array, ks are winding factors %/ 
aain() 
{ 
double design point(), rnd_walk(), swf(), ke, kv, ka, a@inpwr, 
stepsize, randorm(), abs(), freq, rpa; 
int p, iteration, 1, J, best, print_out(), loops; 
FILE tfopen(), fp; 


printf ("\nReading input data from PAM.DAT.. ."); 


fp=fopen("paa.dat’,"r"); /% input seed for random nuabers %/ 
fscanf (fp, "Ad", &seed)s 


fscanflfp, "4d", &p); /‘t input nuaber of pole pairs %/ 

fscanf(fp, "Alf", Sa@inpwr)s /t input aachine power, derived fa ASSET %/ 
ainpwrt=746.0; /t convert to watts %/ 

fscanflfp, “Alf*, &ke)s /t CERs for Effective Weight &/ 


fscanflfp, "Alf", kkv); 

fscanf(fp, "Alf", kk); 

fscanf(fp, "Alf", &rpa); /t machine aax shaft rpa %/ 
fclose(fp)s 


printf ("\nHow aany loops do you want? "); 
scanf("dd", &loops); 
printf ("\n\nDoing prograa calculations . . .\n"); 


for (i=l3 i ¢ Bs i+=2) /t haraonic winding factors %/ 
{ 
ksfij=swé (i); 
} 

freq=rpatp/40.05 /t wax electrical frequency %/ 


/t MAIN BODY OF THE PROGRAM %/ 


for (i=ls 1 <= loopss ++i) 
{ 
printf ("\ni=2d", i); 
stepsize=0. 13 
iteration=0; 
design _point(minpwr, p, ke, kv, ka, freq); 
/t put stuff in the hold array %/ 

while (iteration <= 10) 

{ 

rnd_walk(ainpwr, p, stepsize, ke, kv, ka, freq); 


fen 





/% stagger around %/ 
best=0: /t index to best EW of the lot #/ 
for (j=Ls J<=105 ++4)) 
if (hOj]018) < h€best]018)) 
best=); /t find the best aachine %/ 
if (abs((hlO]018) - h€bestICI8))/(hCOIC18)+.001)) ¢ 0.005) 
/t s@all improveaent in EW t/ 
/% 0.001 takes care of div by zero %/ 
{ 
stepsize/=2.0: 
+titeration; 
} 
else /% transfers best to 0 position %/ 
{ 
for (j=l; ) <= 263 ++)) 
hCO1C3] = h{best)Cj]; 


} 
} 
for (j=l; } <= 265 ++)) 
bli]{j)J=h€best1Cj); /% keep the best sachine %/ 


} 
best=1; 
for (i=l; 1 <= loops; ++i) 
if (bCi]C18)] ¢ bfbest]C18]) 


best=i; /% find and keep the best of the best %/ 
ainpwr /=746.0: /t turn back into hp %/ 
print_out(best, p, minpwr, ke, kv, ka, rpa); /% output to disk file %/ 
fp=fopen(*pas.dat","w"); /¥ output seed %/ 


fprintf(fp,°4d", seed); 

fprintt (ip, “\ndd", p); 

forintflfp,"\ndlf", @inpwr); 

fprinté{tp, “\ndlf*, kel; 

mpeintt (tp, "\nZlt", kv); 

fprintt(fp,"\ndlf", kab; 

fprintf(fp,"\ndlf", rpa; 

fcloseltp): 

} 

/$ END OF MAIN PROGRAM; ALL THAT FOLLOW ARE FUNCTIONS $/ 


double design _point{ainpwr, p, ke, kv, ka, freq) 
/% deteraines a random design point %/ 
double ainpwr, ke, kv, ka, freq; 
int p; 
{ 
double r, js, 1s, Ir, dcore, ds, g, w, 1, xs, va, ph, wtaag, 
pe, i2r, vol, wt, effcy, ew, xsi, xsd, x57, xsal, la, volaag, 
bd, find Ir(), cw, siv, scv, riv, loa, doa, vind siv(), wt_iron; 
extern double sqrt(), randoa(); 
int c=0, d=05 


while (d ‘= 1) ¢ 
while (c '= 4) 


163% 





{ 
r=random() $RMAX: 
if ((2tPItrtfreq/p) < MAX_TIP_SPEED) /& check rotor tip speed %/ 


break; 
} 
w=28F l8treq; 
bd=BREM/ (i + MUR/PC); /% air gap flux, operating point, fro@ magnet 
characteristics %/ 
lr=tind_Ir(bd); /t find magnet slot factor %/ 
while (c '= 1) /% gap dimension %/ 
{ 
g=randoa()$(0.1%r - 6MIN) + BAIN; 
1f (g > 0) 
break; 
} 
1 a=g SHUR; /t radial length of aagnet, see PH-8, 12/15/86 %/ 
ls=randoa()%0.5 + 0.253 /t stator slot factor %/ 
js=randoa() $JSMAX; /t full load stator current density %/ 
dcore=(bdt(r+la)) /(BSAT&p); /% back iron depth %/ 
ds=randoa()%0.9%dcore; /t slot depth ¢ 90% of body depth %/ 


xsl=ksCileks{l); 
x59=(kSC5)8kS(5)/25)s 
xs7=(ks{718ks{71/49): 
xSal=(58]58P]%ds/18); 
xS=(MUS]SESSFR(r+gtla) S(r+gtla) tdssPl$ls 
8 (128(xst + xsd + x57)/(P18q8p) + xsal))/(128(r+la) Sbdsks{1)); 
/% p.u. synch iapedance %/ 

if (xs > 3.0) 

continue; /% don’t want xs too big $/ 
else 

+t: /$ ay escape hatch -- aission coaplete %/ 


l=(minpwrtp) /(2thdswtksCile(rtla) SSF 8dssP)] &(rt+gtla)slssjssPF): 
/t active length %/ 


va=29P 18 (r+la)elewtbdskslilsjstSSFelr+gtla) tdstls/p: /$ VA rating $/ 
siv = find siv(l, r, g, ds, dcore, 1s, la): /t stator iron voluae %/ 
riv = 18Pl&r%r; /% rotor iron voluae %/ 


SCV = 28P1]8(r+gtla)Sdsslss(] + 2.30948P 18 (r+g+]e) SCP/p); 
/t stator copper voluae ¢/ 


cw = scvsDCu; /t total copper weight %/ 
loa = ] + 48(rtgtla); /$ length-over-all &%/ 

doa = 2t(rtgtlatdstdcore)s; /¥ over-all-diameter %/ 
vol = VOLALL#(]oatP]$doatdoa/4) ; /t aachine envelope voluae %/ 

wt = WIALL#(cw + D&(BRNGS#riv + siv)); /t machine weight in kg %/ 
wt iron = D&(riv + siv); /% iron weight only %/ 


ph = 31. B42258BETASfreqtBRisHci Swt_iron/D; 
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/t hysteresis loss in watts, uses iron weight of machine %/ 


pe = (106236. S8NUSBSATSBSATST Ls T i tfregt freqtwt iron) /(RHO8D): 
/t eddy current loss in watts, uses iron weight of machine $/ 


vOlmag=2kPlt&rtislatir; /$ eagnet volume %/ 
wt@ag=volaagsDMAG; /t magnet weight %/ 


i2r = 2. 08SSF&dssCRHOSPLS (1 + 2. 50948P 1S (rtotladsCP/p) Sjstjstlr+qtlad tls; 
/% stator copper loss in watts %/ 
/t revised 1-12-87 %/ 


effcy=(ainpwr)/(@inpwr + ph + pe + i2r); 


ew=wt + kabwtmag + kef(l-effcy) + kvi(vol + katvolaag); 
/$ Effective weight ¢/ 


h{OILiJ=35; hlOJ(2)=freq; hLOILSJ=0; hCOJL4)=r; hlOI(S)=93 
hlOMbJ=dcore; hlOJL[7J=ds; HCOILSI=1e; hCOILII=ls; hlOJL10)=Ir; 
hlOIL11)=vol; hCOJLI2I)=mt: hlOJLiISJ=ph; h{OI(14I=pe; hlOJL15)=12r; 
hlOIL1L6]=va; hCOML1W7IJ=effcy; hlOIL18I=ew: hCOIL19I=15 
HCOWL20)=wtwag; hlOIL21]=volaag; h{OJ(22)=xs3 
h{O}(23I=0.05 /‘t this one makes overhang=0 %/ hl0J(24)=bd; h{OI(25)=loa; 
h{03L26)=doa; 

/% this section just changed all the variables in the “hold” array %/ 


return; 
}} 


double rnd_walk(minpwr, p, stepsize, ke, kv, ka, freq) 
/t walks about design point 10 tises 8/ 

double stepsize, minpwr, ke, kv, ke, freq; 
int p3 

{ 

double r, js, ls, Ir, dcore, ds, g, w, 1, x5, va, ph, pe, i2r, 
vol, wt, effcy, ew, xsi, x50, x57, xSal, le, volmag, wteag, bd, 
find_IrQ), find_brQ), br, overhang, lrat, cw, Siv, scv, riv, loa, 
doa, find sivQ), wt_iron; 

extern double sqrt(), randoa(); 

int 1=1, 3, check=0, ccj; 


while (i <= 10) 
{ 
/t read in the walk around the design point $%/ 


Js=hlOULAIS(L + stepsizes(random() - 0.9)); 
if (js > JSMAX) 
js = JSMAX: /t reset to limit ¢/ 
w=26P I ttreq: 
r=HCOJL4]8(1 + stepsizes(randoa() - 0.5));5 
if ((wir/p) > MAX_TIP_SPEED) 
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continue; /t go to next try if violated $/ 
g=hlOJ(S)IU + stepsizes(random() - 0.9)); 
if (g < GMIN) 
g=OMIN; /% reset to the lieit %/ 


overhang=r andoa()%0. 34; /t magic nuaber "0.34" froa the book 
on per@anent gagnets by Jases Ireland. 
The effect of overhang is to concentrate air 
gap flux, or reduce leakage. %/ 

J=h(01(19); 


]a=gthUR; /t new la, based on the new g ¢/ 
bd=BREM/(1 + MUR/PC): /t agnet operating point flux density %/ 
Ir=find_Ir(bd); /t find new Ir &/ 


Ls=h(O}C9I8(1 + stepsize’(randoa() - 0.5)); 
if (Is > 0,75) 


15=0.75: /% reset to the liait %/ 
ife(1s Ss0.23) 
1570.25; 
for (j=l; } <= 103 ++j) /t aagnet char convergence loop %/ 


{ 
Irat=]/(28(r+qtla))s 
br=bdtfind_br (overhang, rat); 

/t tind effect of overhang on air gap flux ¢/ 
dcore=(br$(r+la))/(BSAT&p): /t ost efficient use of iron t/ 
ds=hCO)(7)8dcore/hlOJ(4); /% keep them in the sae ratio 

as upon exit fa design pt %/ 


x61=ksC1)tks(1)s 

x$9=(k5(9)8ks(5)/25)s 

¥57=(ks07)8k5071/49) : 

xsal=(98]stPlds/18); 

xS=(MUS}SSSSFS(rtgtlae(r+gtladsdssPltls & (128 (x51 + x55 + x57) /(PI&gtp) 
+ xSal))/(128(r+la) Sbreks(1)); 

/t p.u. synch i@pedance t/ 

l=(ainpwrtp)/(28brdatksl i]t (rtl a) eSSFsdstP le (rtgtladtlstjstPr): 
/¥ active length %/ 

} /t end of magnet char convergence loop %/ 


if (xs > 3.0) 


{ 
++check; 
if (check > 25) /t can’t close on decent xs %/ 
{ 
for (ccj=1; cc) <= 103 ++tcc)) 
h€ccj 1£18)=10000000.0; 
/t aake this randoa walk undesirable %/ 
printf(" burp"); 
breaks /% go to next design point %/ 
} 
continues /t can’t have xs too big %/ 
} 
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Va=2PTRlrtla) tl twtbrtksl i eystSSF et (rtgtia) tdstis/p; /t VA rating %/ 


siv = find siv(l, r, 9, ds, dcore, ls, la); /t stator iron volume $/ 
riy = {8Pltrtr; /t rotor iron voluae $/ 
ccv = Z2EPIeirtgtim)sdstist(l + 2.30948PT8(rt+gtin) $CP/p): 

/t stator copper voluge $/ 
cw = scv$DCu; /t total copper weight %/ 
loa = 1 + 48(rtqgtia); /t length-over-all %/ 
doa = 2t(rtgtiatdstdcore) ; /$ over-ali-diameter %/ 
vol = VOLALL8(loatPl$doatdoa/4); /t machine envelope volume %/ 
wt = WIALLS(cw + DS(BRNGS$riv + siv)): /t machine weight in kg %/ 
wt iron = D&{riv + siv)s /t iron weight only %/ 


ph = Si. 8o2258BETAsfreqsBRitHCl wt _iron/D: 
/t hysteresis loss in watts, uses iron weight of aachine $/ 


pe = (106236. 98NUSRSATEBSATSTisTistfregttreqtwt iron) / (RHOSD); 
/t addy current loss in watts, uses iron weight of machine ¢/ 


volaag=2ePTtrtit(l. O+(overhangs (28 (r+gtia))))tlagtir: /t aagnet voluae %/ 
wtaag=voleagsDNAb; /¢ magnet weight %/ 


1dr = 2,08SSFedstCRHOSPis (i + 2,50948PI 8 (r+gtia) SCP/p) Sjstistirt+gtla) tis; 
/t stator copper joss in watts %/ 
/t revised {1-12-87 %/ 


effcy=(ainpwr)/(ginpwr + ph + pe + i2r): 


ew=et + kabwtaag + keS(i-effcy) + kvé(vol + ketvolaag); 
/t Effective weight %/ 


hCillil=jss hCil(2]=freq; hlill3J=n; hCil{4J=r; hi (5)=93 
hlil€4J=dcore; hfilf7J=ds; hCil(8J=les hCil{9]=iss hCil(16]=1r;3 
HCillifd=vol; hlil{i2J=wt; hCil£i3)=ph; h(i J£i4]=pe: hi W{iSJ=12rs 
hCid(t6J=va: HCilli7J=etffcy; hlid€18]=ew: HEide¢9)=1; 


hCil£20)}=wtaag; Ali I(21)=voleaag; hCil(22]=xs;3 

hli }€23)=overhang; h(i J€24)=br; hCil{23]=loa; hli 1£26J=doa; 
/t this section just changed ali the variables in the “hold® array &/ 
fyi: /t go to the next hfil£{) &/ 
check=0; 
} 

return; 

} 


print_out(best, p, ainpwr, ke, kv, ka, rpa) 
int best, p}3 
double ainpwr, ke, kv, ka, rpa; 

{ 

char outfilefi4); 


io 





FILE &fpo, tfopen(); 
iat 1: 


printf ("\nWhat 15 the name of the file where you want the output? *); 
scanf("%s", outfilel; 


fpo=fopen(outfile, "w"); 


fprintt(tpo, "4d", p); 
forintf(fpo, “\nilt", ainpwr); 
fprintft(fpo, "\ndlf", keds 
fprintf (fpo, "\nklf*, kv); 
fprintt(fpo, "\ndlf", kad; 
forintt(fpa,"\nilf*, cpa); 
for (i=l; 1 4= 26; ++1) 
fprinttlfpo, "\ndlf*, blbestilid); 
forintfifpo, “\n"); 
fclose(fpo); 
} 


double find Ir (br) /¥ bracketing %/ 
double br; 

{ 

extern double cos(), sin(); 

int c=0; 

double Ir=0.5, top=1.0, tbr=0.0, abs(), bta=0.05 


while (c!=1) 
{ 
tbr=(2SBREM$ (sin(O.S8Pisir) - cos(O.GsPislr) + 1.0))/PIs 
1f (abs ((br-tbr) /br) <= 0.001) /% check for convergence %/ 
break; 
if (br > tbr) 
{ 
bta=lr; 
Ir=bta + (top- bta)/2.0; 
} 
else ‘tbr < thr &/ 
( 
top=lr; 
lr=bta + (top-bta)/2.0; 
} 
if (abs((top-Ir)/top) <= 0.0005) /% no convergence %/ 
{ 
Arince lo \Nnlt. y) Lf): 
abort("no solution for Ir"); 


i 


} 
return(ir); 
} 
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double find br(overhang, Irat) 
double overhang, lrat; 
{ 
extern double pow(); 
double expon=0. 706133501, factor=0.385576838, enn, yyy; 
/% enn is Ireland’s N, yyy is Ireland’s Y &/ 


enn = factor $ powloverhang, expon);: 
yyy = (Irat + overhang)/(lrat + enn); 
returnlyyy); 

} 


double find siv(l, cr, 9, ds, dcore, Is, Ja) 
double 1, r, g, ds, dcore, ds, Is, 1a; 

{ 

double one, two, three, four; 


one = (rtgtietdstdcore) t(rt+gtiatdstdcore) - (r+gtle)sirtgtla); 
two = 23P18irtgtia)tdstis: 
three = (r+gtlatdstdcore) S(rt+gtlatdstdcore) - (rtqgtlatds) tir+g+tlatds) ; 


four = 18(Pl&one - two) + PL&4s(r+qtla) &three: 
return(four): 
} 


eles 





Table 


41. Listing of permanent magnet efficiency program 


tinclude “stdia.h" 
tinclude "def.h" 


/‘t program name: 


it 


aain()} 
{ 


petf.c to find efficiency of permanent magnet machines %/ 


works with a single machine %/ 


FILE tfopent), fp; 

double r, js, 1s, Ir, dcore, ds, dr, g, w, 1, xs, la, ke, sgeff=1.0, 
va, ph, pe, izr, vol, wt, effcy, ew, xsi, xso, xs/, xsal, ks(8], 
S1V¥, riv, wt iron, find siv(), pmetf, poeff, tett, pgs, pais, 
parpa, gefft=0.98, pcetf=0.99, gr, dhp, rpa, @inpwr, ke, ky, freq, 
wimag, volm@ag, overhang, bd; 

extern double swf(), sqrt(); 


int e=0, f, 


py 4; 


char infilet14); 


siv = find siv(l, r, g, ds, dcore, 1s, Ja); 
riv = ]aPltrtr; 
D¥(riv + siv); 


wt iron = 


/t stator iron voluae ¢/ 
/t rotor iron voluee %/ 
/¥ iron weight only %/ 


printf (“\nCalculates efficiency for a single @achine. \n"); 


while {e '= 1) 
{ 
$ = 03 


printf ("What is the nase of the input file? "); 
scanf("ds", intile); 
tp = fopentintile, “r°); 


fscanfifp, 
fScant itp, 


ainpwr k= 746.0; 


fscanf (fp, 
fscanf (fp, 
fscanf (fp, 
fscanf (fp, 
fscant (tp, 
fscanfifp, 
fscant (tp, 
fscanf (fp, 
fScant (tp, 
fscant(tp, 
fScant (tp, 
fScant (fp, 
fscant tp, 
fscanf (fp, 
fscanfifp, 
fscanf (fp, 
fscant (tp, 
fscant (tp, 
fscanf (fp, 
fscantifp, 


"hd", &p); 
"Alf", &minper)s 
/¥ now in watts %/ 


/¥ input nuaber of pole pairs %/ 


"Alt", &ke)s 
"elt", &ky)s 
oh) Peeks 
iit, kena) ® 
"Alf", &j5); 
"Alf", &freg)s 
nit" 5 Gels 
Shit, ah): 
ehat 46g) 
"Alf", &dcore); 
"Alt", &ds)s 
"Alt, &la)s 
"Alf", &ls)s 
"Alt", kr); 
"Alf", kval) 
Shit vy eML) 
sAlt*, ephi)s 
"hlf®, &pe)s 
Sh tagewueta: 
"Alt", kva)s 
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fscanf(tp, "Alf", keffcy); 
fscantlfp, “41t"*, kew); 
fscant(fp, ““1f", &l); 
fscant(fp, “41f", &ataag); 
fscanfifp, “£1f", &volaag); 
fscanflfp, "Alf", &xs); 
fscanfifp, ““lf", &overhang); 
fscanflfp, "Alf", &bd); 
fcloselfp); 


while (f '= 1) 
{ 
printf("\nWhat is the sustained speed aachine horsepower? °); 
scanf(*Z1f°, &dhp); 
dhp t= 744.0; 
printf ("What is the sustained speed a@achine rpa? °); 
scanf("2lf", &parpa); 
freq = parpmtp/40.0; /t aax electrical frequency ¢/ 
pajs = jstdhp/ainpar; /t PA stator current %/ 
ph = 31.Bb2253BETASfreqtBRitnCitwt iron/D; 
/‘t hysteresis loss in watts, uses iron weight of aachine %/ 
pe = (106234. 93NUSBSATEBSATRTitTitfreqtfreqtwt_iron) /(RHOSD)s 
/t eddy current loss in watts, uses iron weight of aachine %/ 
i2r = 2.08SSFRdssCRHOSPIE (1 + 2.3094SPIR(r+gtl a) SCP /p) &jskisk(rtgtla) sis; 
/t stator copper loss in watts %/ 
paeff = dhp/(dhp + ph + pe + i2r); 
printf("\n Sustained efficiency is 41f", paeff); 


printf("\n\nWhat is the endurance speed machine horsepower? °); 
scanf("Zlf*, &dhp); 

dhp t= 746.0; 

printf("What is the endurance speed aachine rpa? °); 

scanf("Zlf", kpearpa); 

pajs = jstdhp/einper; /¥ PM stator current %/ 

freq = parpatp/40.0; /t wax electrical frequency %/ 
ph = 31. B62258BETASfreqtBRItHCitmt iron/D; 

pe = (106236. 93NUSBSATSBSATSIL eT itfreqtfreqtwt_ iron) /(RHORD)s 

12r = 2.0855F&dstCRHOSPLS (1 + 2. 3094RPL 8 (rtgtla)CP/p) eyskjstir+gtla) tls; 
preff = dhp/(dhp + ph + pe + i2r); 

printf(*\n Endurance efficiency is Z1f*, paeff); 


printf (*\nSa@e aachine? °); 
scanf(*kd", &f)3 


if (f == 0) 
continue; 
else if (f == 2) 
{ 
e= |; 
break; 
} /t end of f-loop t/ 
/t end of e-loap %/ 
} /t end of sain progras %/ 
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double find siv(l, r, g, ds, dcore, 1s, la) 
double 1, r, g, ds, dcore, 1s, la; 

{ 

double one, two, three, four; 


one = (r+gtlatdstdcore) tirtgtlatdstdcore) - (rtgtla)t(rtgtla); 
two = 2kP18(r+gtla)sdstls: 
three = (rt+gtlatdstdcore)t(rt+gtlatdstdcore) - (r+gtlmtds)&(rt+qtlatds) ; 


four = 18(Pl¥one - two) + PIS4s(r+gtla) sthree; 
return(four); 
} 
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Appendix D. Induction machines 


The equivalent circuit for an induction machine is 


shown in the figure. 





Figure 45. Induction machine equivalent circuit 


The stator leakage reactance, Xl, is the sum of the 
stator differential leakage reactance (Xsd), the stator slot 
leakage reactance (Xss), and the stator end turn leakage 
reactance (Xse). In turn, Xsd is the sum of the belt and 
Zig-zag reactances. 

The rotor leakage reactance, X2, is the sum of the 
rotor differential leakage reactance (Xrd), the rotor slot 
leakage reactance (Xrs), and the skew leakage (Xskew). Xm 
is the magnetizing reactance and Rec is the core resistance. 
Ril and R2 are the stator and rotor resistances, respec- 
tively. 

A derivation of the properties of electric machines, 
using Ampere’s Law and the constitutive reiation earlier 
postulated, yields the result that 


iio 





From this important stator inductances are taken. 

Belt leakage reactance (Xbelt) is the sum of the reac- 
tances due to phase-belt harmonics of an "infinite" slot 
winding. In most machines, the most important harmonics 


present are the fifth and seventh, as the third is canceled 


in balanced operation. Then the belt reactance is 
6 Po w Ns? ler Ks 5 ks 7 
aN eee Ces == et =S-=— )* J 
Tg Pp 5 7 


This is the harmonic form of the fundamental mutual reac- 
tance of Appendix B. The winding factor is again the 
product of the pitch and breadth winding factors. 

ZLig-zag reactance is leakage due to all the air gap 
harmonics that would be produced if the winding had one slot 
per pole per phase. For a phase belt of one slot, with each 
slot carrying the same current and equally separated in time 
and space phase, the zig-zag reactance alone would be 
present. Belt leakage occurs because phase belts are ac- 
tually several slots wide. Zig-zag reactance has harmonic 
orders higher than seven, with the same form as Xbelt. No 
even or triplen harmonics will be present. 

The fundamental harmonic of the flux yields the mag- 
netizing reactance, Xs, which can be viewed as that required 


to "energize" the air gap. 


6 Po w Ns 1 gr ks?2 


Figure 46 shows a typical stator slot. The stator slot 
leakage reactance, summing the self and mutual reactances, 


is 


18 wPRo 1 Ns? d2 ds ns 
WO B SSsssssSSS55555 (SS SS aS og as + pNp ) 
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where ns is the number of stator slots and Np is the coil 


throw, in slots. Np = (CP ns)/(6 p), a result available by 


manipulation. we 
d2 
wie _] 
———_—_ 
WS 
Figure 46. Stator slot diagram 
Then, 
3 who 1 Ns2 d2 ds 
Sa aa Ge e == =2 )(1 + CP) 
ns wa 2 WSs 
and 
(1 - ls) 2m (rtg) 
NS = ----c rrr ne ee ne ee ee 


where we have traded the need for the knowledge of the num- 
ber of stator slots for the need to know the width of an in- 
dividual slot. <A reasonable relationship between slot 


dimensions is d2 = ds and w2 = ws. Then, 
3 vo wl Ne? ds (1 + CP) 
ls 20 (rt+g) 


The stator end turn leakage reactance may be estimated 
by treating the two end regions as a single helically shaped 
winding. If the active region of the machine is ignored and 
the helix given air core properties, the inductance can be 


found from standard sources. 


-16 wt Po r= sin* (8w/2) Ns# : pir pir 
P 
where le is the combined length of both end windings. 
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Presume the winding radius to be r, with helix pitch 
Ow = 71/3 and 

le = 2mwr tan(m/3) Np/ns 
Ip’ and Kp’ are the first derivatives with respect to their 
arguments of the hyperbolic Bessel functions Ip and Kp. 
When the three phases are summed, a multiplier of 1.5 will 


be realized. Finally, using the previous results for Np, 


81 wHo Ns? r prr prr 
p 


Stator resistance is Rl = (rho lt Ns)/(As/Ns), as in 
Appendix B, leading to 


6 CRHO Ns¢ (1 + 4(4/3) 20(rt+g) CP/p) 
SSF ds wt (rtg) ls 


When the actual calculation is performed, 1 is not known. A 
guess-and-iterate scheme is used. Iteration continues until 
convergence on 1 is achieved. 

Rotor resistance uses a similar scheme, but the 
presence of rotor bars and end rings instead of turns 


changes it somewhat. <A model of a rotor bar is below. 


we 
op 
ee Az 
dr Oo 
——————— 
wr 


Figure 47. Typical rotor bar configuration 


Induction motor transformer models provide a way to 


find rotor resistances and inductances. The flux density 
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produced by the stator and rotor is 


Ho 
Boe ---— (ke + Fre ) 
= 


The stator mmf is Fs = -j ----------- Ti, where [1 is stator 


current and ws is slip frequency. Rotor mmf is 


where Zs is the rotor surface impedance. The air gap volt- 
age, or voltage across Xm, is defined in terms of the flux 


density and rotor mmf as 


If rotor mmf is now identified with rotor current referred 


to the stator winding, 


3 Ns ks 
Fr = -j cccccc--- Iz, where I2 is rotor current. Then, 
wT p 
Zs W 6 1 Ns2 ks2 
a i  - = 12 
Ws T xr 


Separating Zs into its real and reactive parts and using a 
rotor surface model to describe the relation between rotor 


electric field amplitude and rotor surface current yields 


IP tev 





If magnetic diffusion is ignored, end ring resistance 
can be calculated by comparing losses in the rings and slot. 
The ratio of current densities is found by the ratio of the 
areas. This is then squared and multiplied by the ratio of 


volumes. When summed, 


12 1] Ns? ks? nr rwr 
eee sn ee a ee eon gt } 
nr Tr 1 ler + 
Yslot = CRHO/(dr wr), nr is the number of rotor bars 


and ler is the end ring length, approximated as 
ler = 2m (r - wr/4 - ds/2). The rotor bar width, wr, is 
found by specifying nr and observing that (nr wr) = (27 lr). 

Rotor skew leakage arises when the rotor slots are 
skewed angularly along the axial length to prevent rotor 
cogging. Then, flux does not fully link the bars. When the 
effect is integrated over the rotor, it is seen that 

Xskew = Xm [1 - (2 sin(skew/2)/skew)? ] with the amount 
of skew measured in radians. When typical values of skew 
are input, we see that (Xskew/Xm) * 0.5%. This is a negli- 
gible effect and will be ignored. 

From the previous rotor bar model, it is seen that 

Xslot = wHo (d2/w2 + dr/2wr). Assume that d2 = wr/4 
and w2 = wr/4. Then xslot = wHo (1 + dr/2wr). 

Fitzgerald et al [23] state that only small errors 
result if Re is omitted. Therefore, the core branch may be 
omitted. 


Once the components of the equivalent circuit have been 
calculated, the designer must turn to power and torque con- 
siderations. The internal mechanical power of the machine 
is 


tes sii 
peepee ee: Sees: RZ 


The air gap voltage has been previously defined. The ter- 


minal voltage, Vt, may be found from Vag by means of a volt- 
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age divider and a Thevenin equivalent circuit developed. 





Figure 48. Induction machine Thevenin equivalent circuit 


Ri + j(X1 + Xm) 
Rel + jXel = (Ri + jX1) in parallel with jXm. Then operat- 
ing point torque is 
1 3 Via? (R2/slip) 
Ww (Rel + R2/slip)? + (Xel + X2)4 


Torque is a maximum when the power delivered to (R2/slip) is 
a maximum. By matching load and Thevenin impedances, the 


power is a maximum and a slip-at-maximum-torque is found 


(Rel? + (Xel + X2)2 )9-5 


and the corresponding torque is 


Ww Rel + (Rel? + (Xel + X2)2 )0-.5 


Typical induction motors have the ratio between maximum and 


operating point torque as 


tis 





which is used to find the machine active length to use in 
the circuit component calculations. The equation for rated 
torque can be manipulated to yield a quadratic expression 
for operating point slip, or a Newton’s method convergence 
can be used to find operating point slip. Convergence on 
slip and active length through Newton’s method is used to 
generate the machines of this thesis. 

Finally, rotor copper losses are (1 - slip)P; stator 


copper losses are found using [1 and Rl. 
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Table 42. Listing of induction machine design program 


finclude "stdio.h’ 
tinclude "def.h" 


/t program name: ind.c for induction machines, 4/19/87 %/ 


long int seed; /t start point for randoa nuaber generator %/ 
double 6€26)(33], h€11]£33), xs5(42); 
/t bis “best® array, h is “hold” array, xs are winding factors %/ 
oain() 
{ 
double design _point(), rnd walk(), swf(), ke, kv, minper, 
stepsize, randoa(), abs(), freq, rpa, ksl43); 
int p, iteration, 1, J, best, print_out(), loops, flag; 
FILE fopen(), tp; 


printf("\nReading input data from IND.DAT.. .")} 


fp=fopen("ind.dat","r"); /¥ input seed for randow nuabers %/ 
fscanf (fp, "Ad, &seed) 5 

fscantltp, "4d", &p); /t input nuaber of pole pairs ¢/ 

fscanf(fp, “Zlf", &minpar); /t input machine power, derived fa ASSET %/ 
ainpwr t=746.0; /t convert to watts ¢/ 

fscant(fp, "Alf", &ke); /t CERs for Effective Weight %/ 

fscant(tp, "Alf", &kv)s 

fscant(tp, "Alf", &rpa); /t machine @ax shaft rpa %/ 
fclase(fp); 


printf("\nHow gany loops do you want? "); 
scant("4d", &loops); 
print#(*\n\nDoing program calculations . . .\n"); 


for (i=ks i ¢ 425 1+=2) /t haraonic winding factors %/ 
{ 
kSCi J=swf (i); 
xSCi] = (ksCiltksli))/(itids 
} 
freq=rpatp/40.0; /$ max electrical frequency %/ 


/t MAIN BODY OF THE PROGRAM %/ 


for (i=ls 1 <= loops; ++i) 
{ 
printf("\nid", 1); 
stepsize=0. 1; 
iteration=0; 
flag = 03 
design point(minpwr, p, ke, kv, freq, &flag); 
/t put stuff in the hold array ¢/ 
if (flag == 1) 
{ 
hCOIL1B] = 10000000. 0; 
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best = 0; 
} 
while ((iteration (= 10) && (flag != 1)) 


f 
t 


rnd_walk(ainpwr, p, stepsize, ke, kv, freq); 
/t stagger around %/ 
best=0; /% index to best EW of the lot $/ 
for (j=l; j<=103 ++)) 
{ 
if (hCjI018) < 1.0) 
continue; 
if (h€jJC18) < h€bestI(18)) 
best=); /t find the best machine ¢/ 
} 
if (abs((h{0]018) - h€best]£18))/h£01£18)) < 0.005) 
/‘t small improvement in EW %/ 
{ 
stepsize/=2.0; 
++i teration;: 
} 
else /% transfers best to 0 position %/ 
{ 
for (j=ls j <= 323 ++)) 
h(O]Cj] = h{bestjCj); 


} 
} 
for (j=l; j<= 323 +43) 
bfil€jI=hlbestilj); /$ keep the best sachine $/ 
} 
best=1; 
for (i121; 1 <= loopss ++i) 
( 
if (bfi)£18) ¢ 1.0) 
continue; 
if (bCiI£187 < blbesti£18)) 
best=i; /t find and keep the best of the best %/ 
} 
ainpwr/=746,05 /$ turn back into hp %/ 
print_out(best, p, minpwr, ke, kv, rpa); /t output to disk file &/ 
fp=fopen(*ind. dat", "w"); /$ output seed %/ 


fprintf(fp, "4d", seed); 

fprintf(fp,"\nad", ps 

fprintf(fp, "\ndlf", minpwr); 

fprintf (fp, *\ndlf*, ke); 

fprintf(fp, "\ndlf", kv); 

fprintf(fp,"\ndlf", rpm); 

fclase(fp); 

) 

/$ END OF MAIN PROGRAM; ALL THAT FOLLOW ARE FUNCTIONS 3/ 


double design point(minpwr, p, ke, ky, freq, fla) 


IRS ee 





/¥ deteraines a randoe design point %/ 
double ainpwr, ke, kv, freq; 
mrp, ‘tia; 
{ 
double r, 1s, Ir, dcore, ds, dr, wr, g, w, we, |, i2rr, 
igs XUCIt, kez, X55, ¥SE, XI, xd, xrs, 
xa, vag, vt, via, rel, xel, tax, radical, sb, ifa, i{b, 
x2, Va, ph, pe, i2r, vol, wt, effcy, ew, le, teinper, 
Sa, SC, CW, SIV, SCV, riv, rcv, loa, doa, r2seax, ttaax; 
double ler, fl, £2, £3, slip, sa@ax, jr, js, 
11, i2, xslot, wt_iron; 
extern double sqrt(), randoal), cos(), sin(), cosh(), sinh(), vrat(), 
floor(), besip(), beskp(), tan(), find_sivl), swf(); 
int cslip=0, nr=71, ccr=1, ccg, ccs=1, k; 


w=23F 1$freq: /¥ synchronous frequency in rad per sec ¥/ 
wa=W/D; /$ aech angular velocity in rad per sec %/ 
P= 10.0; /$ initial quess on length $/ 

taax = ainpwr/(watPS]); /¥ pull-out torque ¢/ 


while (cslip != 1) 
{ 
r=randoa()SRMAX; 
if ((29PU%rkfreq/p) ¢ MAX_TIP_SPEED) /% check rotor tip speed $/ 


break; 
} 
lr=random()%0.5 + 0.25; ‘¥ rotor slot factor %/ 
dr=randoa() $r/3.0; /¥ slot no deeper than 332 of rotor radius %/ 
dcore=(BR%r)/(BSATIp); /$ back iron depth %/ 
ds=r andoal)%0.9%dcore; /¥ slot depth ¢ 902 of body depth to start %/ 


g=randoa()3(0.1%r - GMIN) + GMIN; 
if (q ¢ GMIN) 


q = GMIN; 
ls=randow()%0.9 + 0.253 /$ stator slot factor %/ 
fl=0.; 
f2=0.3 
f3=0.; 
for (k=f; k ¢ 1003 ++k) 
{ 
fi = (t./(p + k8nr)) & (L./(p + k&nr))s 
f2 = ({./(p - k$nr)) & (L./(p - kSnr))s 
forta t) + £23 
} 
while (ccr <= 3) /¥ start convergence loop on r $/ 
{ 
if (ccs > 3) /¥ we’ve had trouble with jr %/ 
{ 


r $= jr/JSMAXs 
if (r > (MAX_TIP_SPEEDSp/w)) 
r = MAX_TIP_SPEEDSp/w; 
if (r > RMAX) 
r = RMAX; 
Ir t= jr/JSMAXs 
et ee. 7 a) 
lr = 0.793 
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) 


wre = (B88P 108i 8(r - dr/2.))/(4tnr + Pieler); /t rotor slot width %/ 

ler = 28Pi8(r - wr/4.0 - dr/2.0); /t end ring length 8/ 

ccs = |; 

while (ccs <= 3) /t start convergence loop on ds and dr $/ 

{ 

ccg = 1; 

while (ccg <= 10) /t start convergence loop on length %/ 
{ 
ri = (O8CRHOS(] + 2. 3094SP1 8 (r+) 8CP/p)) /(SSF&dssPit(r+g) $s); 


r2 = (6exsliSCRHORL et) + (2trtrtIr) /(lertp8o)))/(Plersirsdr); 

xbelt = (d8wSlsnUsrs(xsl3) + xsl7J)/ (PL sgtptp)); 

R22 = (ObwSl8hUsr 8 (xsCll) + xsC13) + xsCi7) + xsC19) + x5(23) + 
xS029) + xs{31) + xs(€37) + xsl41]))/ (P18 gtptp); 

KSS = (SUMUSwELedsB(L+CP)) /(1s82eP is ir+g))s 

Te = (Pl&rsCPstan(Pl/3))/(38p): 

xse = (278wS18hUsrtrt(-beskp(p, psPltr/leltbesip(p, pePitr/lel))/ 

(Pltletptp); 

xi = xbelt + xzz + x55 + xse: 

xrd = (OSmMUSwhltxsliJtref3) /(PI8g); 

ers = (OexsCl]ewSTSMUs (wrtO. 8dr) )/ (ir sPlt(r - (dr /2.0) - (wr/8.0))); 

no= ond + xrss 

xa = (OUNUSweL er sxslil)/(Pltgtpsp); 

vag = 2r8BRew8l tswt (i) /p; /t air gap voltage %/ 

rel = (xetxatri)/(ritri + (xltxa)e(xitxa)): 

/% thev equiv resistance %/ 
xel = (xasrdtrd + xidxd + xidxe))/(ritrd + (xdtxa) i (xitxa)); 

/t they equiv inductance ¢/ 
r2smax = sqrt(relérel + (xeltx2)%(xeltx2)); ‘$ r2 at sax %/ 
vt = vag/vrat(r2smax, x2, x1, xa, rl); 

/‘t terminal voltage at saax %/ 
¥la = (vt8xe)/sqrt((xi + xade(xd + xa) + (ritri)); 

/t thevenin equivalent voltage %/ 
ttmax = (L.o%viatvial/(wat(reltsogrt (reltrel+t(xel+x2) t(xel+x2)))); 

/t test gaxiaua torque %/ 
if (abs((teax-tteax)/taax) <= 0.005) 


break; /t we have convergence %/ 
] $= taax/tteax; /‘t reset I, reiterate %/ 
++CCQ;3 
lf (ccg > 10) 
{ 
tfla = 1; 
orintf("\n flag set on length"); 
return; 
) 
} /t end convergence loop on length %/ 


smax = r2/sqrtireitrel + (xeltx2)8(xeltx2)); 


Slip = s@ax/3.0; /t starting point for slip converge %/ 
cslip = 1; 
while (cslip <= 20) 

{ /t start convergence loop on slip %/ 
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tainper = (38viatvlatr2/slip)/((reltr2/slip)t(rel+r2/slip) + 
(xeltx2) 8(xel+x2)); 
if (abs((minpwr-tminpwr)/ainpwr) <= 0,005) 


break: /% we have convergence %/ 
Slip §= tainpwr/sinpur; 
++c511p3 
} /% end convergence loop on slip $/ 


12 = vagtsqrt((r2/slip)t(r2/slip) + x28x2)/((r2/slip)&(r2/slip) + x2%x2); 
jr = ($812) / (PIR  rtdr tr tRSF) ; /t rotor current density %/ 
if (jr > JSMAX) 


f 
\ 


dr t= jr/JSMAX: 
if (dr > (r/3)) 
dr =r/3; /% reset to liait %/ 

) 
bla = (r2/slip) @(xm+x2) - x28r2/slip; 
11b = x28(xatx2) + (r2/slip)8(r2/slip); 
1] = vagtsqrt(ilatila + ribsitb)/(xmd(x2ex2 + (r2/slip)t(r2/slip))); 
js = (3.0811) /(PI8(r+g) tdstlstSoF); 
if (js > JSMAX) 

ds t= js/JSMAX; 
if ((js <= JSMAX) && (jr <= JSMAX)) 

break; 
+4053 
} /% end convergence loop on ds and dr %/ 
if (ccs <= 9) 

break; 
t+CCr; 
Peeccr > v) 

£ 


ttla = 1; 
printt(*\n flag set on rr"); 
returns; 
} 
} /¥ end convergence loop on r %/ 


/t calculations %/ 


saax = r2/sqrtireltrel + (xel+x2)%(xel+x2)); 
Siv = find siv(l, r, g, ds, dcore, 1s); /t stator iron voluae %/ 
riv = ]&P)&ré(r - 2tdrtlr); /t rotor iron voluae %/ 
SCV = 2EP](rtg)tdstlst(] + 2.30948P1%(r+q) 8CP/p); 

/t stator copper voluae %/ 


rcv = wrtdrt(nrtl + 2tler); /t rotor copper voluae %/ 

cw = (rcv + scv) $DCU; /t total copper weight %/ 

loa = 1] + 48(rtq); /t length-over-all %/ 

doa = 28(r+gtds+dcore); /t over-all-diaaeter %/ 

vol = YOLALLS(LoatPI$doatdoa/4); /t aachine envelope voluae %/ 
wt = WALLS (cw + DS(BRNGS#riv + siv)); /% aachine weight in kg %/ 
wt_iron = Di{riv + siv); /t iron weight only %/ 


ph = SL. Bd225sBETAStreqsBRIsHCl swt iron/D; /t hysteresis loss in watts, uses iron weight of 
machine %/ 
pe = (106236. 9tNUSBSATSBSATST1¢Tlttreqt freqtwt iron) /(RHOSD) ; 
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/t eddy current loss in watts, uses iron weight of aachine &/ 
12r=3. 081 181 18r; /$ stator copper loss in watts %/ 
l2rr = slipkeinpwr; /$ rotor copper loss in watts %/ 
vt = sqrt((vagtiltrd)t(vagtiigrd) + (idsehyeidssid)s 9 /& rated voltage ¢/ 
va = Stythil; /t VA rating at terainals %/ 
effcy=(minpwr)/(ainpwr + ph + pe + 12r + i2rr); ew=wt + kes(l-effcy) + kvtvol; (GEE: 
fective weight $/ 


hfOlCij=ii; hCOM2)=freqs; hlO]{3]=w; hCO(4]=r; hCOL5)=q3 
hCOWlS}=dcore; hlOIL7)=ds; h(O](8)=drs hCONLPI=153 hCOICIOI=Ir; 
HCOWL11J=vols = ACOICA2Z)=wt; = ALO 3)=ph; h(01(14)=pe; h(0)(15)=i2r; 
h(01016)=va; hCOI17)=effcys HCOICL8)=en; h(OIC19I=1; h(0)(20)=12rr; 
h(O21]=saax; ALOI(22]=tmax; hlO](23]=vt; HCOW(24]=slips hlOI(25]=vaqs 
hCOV(26]=r1; HCOWE27J=x13;  nlOJ(28)=xa; hCO]C29)=x23 h(01(30)=r2; 
HCOWUSII=loa; = h(09£52)=doa; 

/t this section just changed all the variables in the “hold® array %/ 


return: 
} 


double rnd_walk(einpwr, p, stepsize, ke, kv, freq) 
/t walks about design point 10 tiees %/ 

double stepsize, ainpwr, ke, kv, freg; 
int p; 

{ 

double r, 1s, Ir, dcore, ds, dr, wr, g, #, wa, 1, i2rr, ila, i1b, 
ri, r2, xbelt, xzz, xss, xse, x1, xrd, xrs, r2s@ax, ttaax, 
x@, vag, vt, via, rel, xel, teax, radical, sb, tainper, 
x2, va, ph, pe, 12r, vol, wt, effcy, ew, le, cw, Siv, SCV, 
Sa, SC, rlv, rcv, loa, doa, jr, js, 
ler, f1, #2, #3, slip, smax, il, 12, xslot, wt_iron; 

extern double sqrt(), randoe(), cos(), sin(), cosh(), sinhQ), vratQ), 
besip(), beskp(), tan(), find siv(), swf(); 

int d, 1=1, ccg, cer, ccs, nr=71, k, clips 


#1=0.3 

#2=0.5 

#3=0,3 

for (k=15 k ¢ 100; ++k) 
{ 
fl = (1./(p + ktnr)) & (1./(p + ktnr))s 
{2 = (1./(p - kenr)) & (1.7 (p - k&nr))s 
Porta state f2: 
} 

while (i <= 10) /t ten steps around the design point ¢/ 
{ 


/$ read in the walk around the design point %/ 
w=23P 18 freq: 
r=hCO](4]8(1 + stepsizet(randoa() - 0.5)); 
lf ((wir/p) > MAX_TIP_SPEED) 
continue; /$ go to next i-loop if violated %/ 
g=h(O}(9]8(1 + stepsizek(randoa() - 0.5)); 
if (g < GMIN) 
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g=GMIN; /t reset to the limit %/ 
dcore=(BR8r)/(BSATSp); /‘¥ aost efficient use of iron $/ 
ds=H(OIL7]¢(1 + stepsizes(randoa() - 9.5)); 
dr=hCO]CB]8(1 + ctepsizetirandoa() - 0.5)); 
elope feo.) 
dr = 7/303 /‘t reset to the liait &/ 
Ls=hCOIC9I8(L + stepsizes(randoa() - 0.5))s 
if (1s > 0.79) 
1s=0.75; /t reset to the limit %/ 
if (ls ¢ 0.25) 
15=0, 29; 
br=h(O](10]8(1 + stepsizet(randoa() - 0.9)); 
Taha Oe o} 


lr=0.753 /‘t reset to the limit %/ 
Vile 40. 25) 
lr=0.29; 
wa-w/D; /t aech angular velocity in rad per sec %/ 
1 = h0)f19); /t length starting point %/ 
tax = h{0][22); /t pull-out torque %/ 
acs = };3 
ccr = 1; /‘t first time thru, no adj) in %/ 
while (cer <= 9) /% start convergence loop on r %/ 
{ 
mt (ees > 5) 
{ 


r = jr/JSMAX; 
if (r > (MAX_TIP_SPEEDSp/»)) 
r = MAX TIP SPEEDSp/w; 
if (r > REMAX) 
r = RMAX; 
lr t= jr/JSMAXs 
if (lr > 0.79) 


lr = 0.73; 

} 
wre = (BSPL8iek(r - dr/2.))/(48nr + Pltlr); /t rotor slot width %/ 
ler = 24PU(r - wr/4.0 - dr/2.0); /t end ring length %/ 
eos = 1; 
while (ccs <= 3) /t start convergence loop on ds and dr$/ 

{ 

cceg = I: 

while (ccg <= 10) /% start convergence loop on length %/ 

{ 


ri = (O8CRHOS(L + 2. 30948P18(rtg) SCP/p)) / (SSF edstP 18 (r+q) S15) ; 

r2 = (6%xsClL]sCRHOs ls (1 + (2trtrtlr)/(lLertptp)))/ (Pier alrtdr) 5 

xbelt = (dtwSLthUtrd(xs5) + xsl7))/(PL¥gtptp)); 

x2Z = (OdwSl seer d(xsCil) + xs013) + xsfl7) + xsf19) + xsf23] + 
x6(29) + xs{3l] + xsf37) + xs€4l)))/(PL&gtp§p); 

xSS = (SIMUSwWELSdst(L+CP))/ (Ls828Pl8(r+q)); 

le = (Plarschstan(Pl/3))/ (3p); 

xse = (27ewtl thUtrtrt(-beskp(p, paPltr/ledtbesip(p, p¥Pltr/le)))/ 

(Pltlespsy); 
xl = xbelt + x22 + xss + x58; 
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Ard = (O8MUtwdltxsll trsf3) / (Pita); 

ees = (OExsCL)BwALOMUR (we tO. S8dr))/ Ur tPislr - (dr /2.0) - (we /8.0))); 
x2 = xrd + Xxrs; 

xm = (OSMUSwRL or dxell)) /(PLsgtptp); 


vag = ZtrBRewsisswt (i) /p: /t air gap voltage %/ 
rel = (xméxmtrd)/(rttrd + (xitx@) &(xdtxa))s 

/t thev equiv resistance %/ 

Kel = (xeddrdtri + xidxl + xidxmd)/(rttrd + (xdtxm) st (xitxea)): 

/% they equiv inductance %/ 
e2emax = sqrt(reigrel + (xel+x2) 8 (xeltx2)); /%r2 at smax %/ 
vt = vag/vrat(r2smax, x2, xl, xa, rl); 

/% terminal voltage at smax %/ 
Vila = (vttxe)/sqrt((xt + xmdt(xt + xm) + (rltri)); 

/% thevenin equivalent voltage %/ 
tt@ax = (L.o8viatvia) / (wat (reltsqrt(reltrelt (xelt+x2) &(xettx2)))); 

/% test maxisua torque %/ 
if (abs((tmax-tt@ax)/teax) <= 0.605) 


{ 
ccg = 15 
break; /% we have convergence %/ 

} 
] 8= tmax/ttmax; /¥ reset 1, reiterate ¥/ 
FFCCQ} 
} /t end convergence loop on length %/ 
if (ccs > 9) 
break; 


Smax = r2/sqrt{rettrel + (xelt+x2) 8 lxeltx2)); 
Slip = saax/3.0; 
/% starting point for slip converge %/ 


cslip = 1; 
while (cslip <= 20) 
{ /% start convergence loop on slip %/ 


tainper = (S8vlatviatr2/slip)/((rettr2/slip)&(rettr2/slip) + 
(xeltx2) $(xeltx2))s 
if (abs ((@inpwr-tainper)/m@inper) <= 0.009) 


break; /¥ we have convergence %/ 

Slip $= teinper/ainper; 

++cS]ip; /$ no gore than 20 tries %/ 

} /% end convergence loop on slip $/ 
i2 = vagssqrt((r2/slip)s(r2/slip) + x28%2)/((r2/slip)&(r2/slip) + x2%x2); 
jr = (3812) / (PTS ir 8dr tr tRSF); /t rotor current density %/ 
if (jr > JSMAX) 

{ 


dr %= jr/JSMAX; 
if (dr > (r/3)) 
dr =r/3; /3 reset to limit %/ 

} 
ila = (r2/stip) 8 (xatx2) - x28r2/slip; 
ilb = x28(xa@tx2) + (r2/slip)t(r2/slip); 
il = vag$sqrt(ilatila + 11b8itb)/(xme(x28x2 + (r2/ship)8ir2/slip))); 
js = (3.0811) /(PIS(rtg) tdstlssSSF); 
if (js > JSMAX) 
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ds %= js/JSMAX; 
if ({js <= JSMAX) && (jr <= JOMAX)) 


break: 
++CCS; 
) /$ end convergence Joop on ds and dr %/ 
if ({ccs ¢= 5) ts (ceg > 10)) 
break; 
++CCr} 
} /t end convergence loop on r %/ 


if ((ccg > 10) tt (ccs > 5) te (cslip > 20) ti (cer > 5)) 
{ /t no design convergence %/ 
for (d=i3 d ¢=10; ++d) 
hCd]C(18) = 10000000.0; 
printf(“burp °); 
break; 


} 


/¥ calculations %/ 
Saax = r2/sqrtirelsrel + (xeltx2)&(xeltx2)); 


Siv = find siv(l, r, g, ds, dcore, 1s); /t stator iron volume %/ 
riv = 18Plars(r - 2kdrtlr); /t rotor iron voluae %/ 
ecv = Q8PU8(rtg)8dsklsk() + 2, 30948PT8 (r tg) 8CP/p); 
/% stator copper volume %/ 
rcy = wretdrt(nrt] + 2%ler); /t rotor copper voluae %/ 
cw = (rcv + scv) DCU; /t total copper weight %/ 
loa = ] + 48(rtq); /$ length-over-all %/ 
doa = 28(rt+gtdstdcore); /t over-all-diaseter %/ 
vol = VOLALLS(loasPI%doatdoa/4); /¥ machine envelope voluee %/ 
wt = MTALLS(cw + DS(BRNOSSriv + siv))s /$ gachine weight in kg %/ 
wt iron = D&(riv + siv); /¥ iron weight only %/ 
ph = SL, BO2253BETASfreqsBRisHCiswt_iron/D; /t hysteresis loss in watts, uses iron weight of 
machine %/ 


pe = (106236. 98NUSBSATSBSATSTI8T18freqtfreqtwt iron) /(RHOSD); 
/¥ eddy current loss in watts, uses iron weight of machine %/ 
i2r = S.08iLtiltrt; /$ stator copper loss in watts %/ 
i2rr = slipssinpwr; /t rotor copper loss in watts %/ 
vt = sart((vagtilSri)s(vagtilgrd) + (ibsxddsidsxtdds  /¥ rated voltage %/ 
va = 3tvttil; /$ VA rating at the terminals %/ 
effcy=(ainpwr)/(ainpwr + ph + pe + i2r + i2rr)s ew=wt + ket(i-effcy) + kv8vol; /t Ef- 
fective weight %/ 


HCidli)=il; hCill2)=freq; hlill3J=w; hCill4)=r; hCiI{SJ=93 
hCil{6J=dcore; hlill7)=ds; hi J(8)=dr; hCidl(9)=1s3 hCiI(tO)=1r; 
HCiNCitJ=vol; mCi dC2)=wt; = ACA ILASI=ph; = Ca t4d=pe; = li 015)=12r; 
hCillidj=va; hCidli7J=effcys hlil18)=en; hCiJCI9J=)5 hCij£20)=1 2rr; 
HOiNC2))=smaxs hlidl(22]=teaxs ACiW2SI=vts; =ACiI(24]=slip; hliI(25)=vaq; 
hCid(26)=r 1; RCIVCO27J=xt; ACA IC2B)=xm; ALL E29) =x 2; hliJ(30)=r2; 
HCiMS1=loa; =hlidl32]=doa; 

/3 this section just changed all the variables in the “hold” array %/ 

mie 

} 


return: 
} 


HO? 





print _out(best, p, @inpwr, ke, ky, rpm) 
int best, p; 
double ainpwr, ke, kv, rpas 

{ 

char outfileli4); 

FILE &fpo, tfopen(); 

int 13 


printf ("\nWhat 1s the nage of the file where you want the output? "); 
scant (*4s", outfile); 


fpo=fopen(outfile, “w*); 


fprintflfpo, "4d", ps 

fprintt(fpo, “\nklf*, minper); 

fprintt(fpa, “\ndlf", ke); 

fprintf(fpo, "\ndlf", kv); 

fprintfifpa, “\ndlf", rpe); 

for (i=ly i ¢= 32: +41) 
fprintf(fpo, “\nde” b{bestI(iJ); 

fprintf(fpa, "\n")s; 

fclose(fpo); 

\ 


4 


double find siv(l, r, q, ds, dcore, 1s) 
double 1, r, q, ds, dcore, 1s; 


f 
i 


double one, two, three, fours 


one = (r+gtdstdcore)s(rtgtdstdcore) - (r+g)s(rtq); 
two = 23F18(r+g)8dstis; 
three = (r+gtdstdcore)&(r+gtdstdcore) - (r+gtds)&(r+qtds); 


1) 


four = 18(PI%one - two) + PI&4e(r+q)&three: 
return (four); 


} 

double vrat(r2smax, x2, xl, xm, rt) {% ratio of vag/vt %/ 
double r2smax, x2, x1, xa, ri; 

{ 


double a, b, c, d, vrat; 
extern double sqrt(); 


r2smax8(xmtxt) + rid(x2+xa); 

= rltr2smax - xmtx2 - xld(x2txe)s 
atr2smax - 8x2; 

d = b3r2smax + atx2: 


7 or ew 


vrat = (xetsqrt(ctc + dtd))/(ata + b&b); 
return(vrat);} 
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Table 43. Listing of induction efficiency program 


tinclude “stdio.h* 
Hinclude “def.h* 


/$ program mage: ieff.c to find efficiency of induction machines %/ 

/% works with a single aachine $/ 

gain() 

{ 

FILE %fopen(), %fp; 

double r, ls, Ir, dcore, ds, dr, g, w, 1, i2Zrr, va, ph, pe, i2r, vol, 
wh, effcy, ew, Siv, riv, wt_iron, find siv(), paeff, pgetf, teff, 
parpa, geff=0.98, pcetf=0.99, gr, dhp, rpm, minpwr, ke, kv, freq, 
sgeff=1.0, wr, wa, ri, r2, xl, x@, vag, vt, via, rel, xel, radical, 
x2, Sa, Sb, Sc, Slip, slip!, slipz, smax, 11, 12, tmax; 

extern double sqrt(); 

int e=0, f, p, 13 

char infile(t4]: 

printt(*\nCalculates efficiency for a single @achine. \n"); 


while (e '= 1) 
= 0; 
printt(“What is the nage of the input file? “°); 
scanf("Zs", infile); 
fp = fopen(intile, “r"); 


fscantifp, "4d", &p); /¥ input nuaber of pole pairs $/ 
fScanfitp, “4l¢", &einowr); 
binpwr t= 746.0; /$ now in watts %/ 


fscanfifp, "41lf", &ke)s 
fscantlfp, "Alf", &kv); 
fscant(tp, "Alf", &rpad; 
fscantltp, "Alf", &il); 
fscantiftp, "Zlf", &freq); 
fscantltp, "Alf", &w); 
fscantitp, “4lf", &r); 
fscantitp, "Alf", &g)s 
fscantifp, "Zlf", &dcore); 
fscantlfp, "Alf", &ds); 
fscantitp, "Alf", &dr)s 
fscantifp, "Alf", &ls); 
fscanfiftp, "“lf", &lr); 
fscantifp, "Alf", &vol); 
fscantiftp, “Zlt", &wt); 
fscantitp, "Alf", &ph); 
fscantltp, “4lf", &pe); 
fscanflfp, "4lf", &i2r); 
fscanfifp, ““lf", &va); 
fscanfifp, "Zlf", keffcy); 
fscanfifp, "Zlf", &ew)s 
fscantitp, “Alf*, &l); 
fscantifp, "“1f", &i2rr)s 
fscantlftp, "Alf", &s@ax); 


peau 





fscantltp, “Alf", &twax): 
fscantitp, "Alt", &vt); 
fscantlfp, ““lf", &slip; 
fscantlfp, “4l¢", &vag)s 
fscant(fp, “Zl#", &ri); 
fecantitp, "Alf", xi); 
fscantlfp, "Zlf", &xa)s 
fscanflfp, "Alf", &%2)3 
tscantltp, "Alf", &r2); 
fcloseltp); 


while (f != 1) 

{ 

printf("\nWhat 15 the sustained speed aachine horsepower? “); 
scanf("Zlf", &dhp); 

dhp t= 746.0; 

printf("what is the sustained speed machine rpa? “); 
scanf("Zlf", &parpa); 

freq = perpatp/60,0; /+ PR frequency %/ 


Via = (vttxa)/sqrt( (xi + xm)S(xi + xm) + (rdtrid); 
/t thevinin equivalent voltage %/ 
rel = (x@sxetri)/(riari + (xitxe)(xitxa))s 
/% they equiv resistance %/ 
(xeb(ritrd + xidxt + xidxa))/(ridri + (xitxe) t(xitxa)); 
/% they equiv inductance $/ 

sa = reltrel + (xel+x2)¢lxeltx2); /% pieces of slip quadratic %/ 
sb = 2treitr2 - (38)tviatviatr2/dhp); 
sc = rztr2; 
radical = sb¥sb - 4%sa%sc; 
if (radical ¢ 9.0) 

abort("\nbot a negative radical in the slip eqn."); 
slipi = (-sb + sart(radical))/(2%sa); 
slip2 = (-sb - sqrt(radical))/(2%sa); 
if ((slipl] ¢ 0.0) && (slip2 ¢ 0,6)) 

abort("\nTwo negative slips."); 

/¥ now we will use the smallest positive slip $/ 

else if ((slipi ¢ 0,0) && (slip2 > 0.0)) 

Slip = slip2s 
else if ((slipi > 0.0) && (slip2 ¢ 0.0)) 

slip = slipl; 


xe] 


+ 


else 
Slip = (slipi > slip2) ? slip2 : slipl; 
/¥ tslip=min(slipi, slip2) ¢/ 
1f (slip > sax) 
abort("\n slip is ore than smax."); 


Siv = find_sivil, r, g, ds, dcore, 1s); /$ stator iron voluae §/ 
riv = 18Plér'lr - 28drsir); /$ rotor iron voluee %/ 
wt_ iran = De(riv + siv); /¥ iron weight only %/ 

ph = Si, B62259BETASfreqeBRitHCiswt iron/D; 


/t hysteresis loss in watts, uses iron weight of sachine %/ 
(106236. 98NUSBSATEBSATST18Tisfreqttregtswt iron) /(RHOSD) ; 
/$ eddy current loss in watts, uses iron weight of sachine %/ 


pe 
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12 = sart((slipsdhp)/(3%r2s(1-slip))); /t load current ¢/ 


tL = 12 + vag/xa; /$ 12 plus g@agnetizing current %/ 
y2r=S. 081181 14rl; /% stator copper loss in watts $/ 
t2rr = sliptdhp; /t rotor copper loss in watts &/ 


paeff = dhp/(dhp + ph + pe + 12r + i2rr); 
printf("\n Sustained efficiency is %1f", pmeff); 


printf("\n\nWhat 1s the endurance speed machine horsepower? "); 
scant("Zlf", S&dhp); 
dhp t= 744.0; 
printf("What is the endurance speed sachine rpa? "); 
scanf("Z1f", &parpa); 
freq = parpats/60.0; /t PA frequency %/ 
sb = 2treltr2 - (Stl tviatviatr2/dhp); 
radical = sb&sb - 4¥safsc; 
if (radical ¢ 0.90) 
abort {*\nbot a negative radical in the slip eqn."); 
slipl = (-sb + sart(radical))/(2%sa); 
slip2 = (-sb - sqrt(radical)}/(28%sa); 
if ((slipl <¢ 0.0) && (slip2 < 0.0)) 
abort("\nTwo negative slips. "); 
else if ((slipl ¢ 9.0) && (slip2 > 0.0)) 
slip = slip2: 
else if ((slipl > 0.0) && (slip2 ¢ 0.0)) 
slip = slip]; 
else 
slip = (slip! > slip2) ? slip2 : slip; 
if (slip > s@ax) 
abort("\n slip is aore than smax."); 


ph = 31. B862258BETASfreqsBRitHC{twt_iron/D; 

pe = (106236. 98NUSBSATSBSATET 1 8Tltfreqtfreqtwt iron) /(RHOMD); 

12 = sqrt((slipsdhp)/(3tr2s(i-slip))); /$ load current %/ 

if = 12 + vag/xa; /ti2 plus agnetizing current %/ 
12r=S. O81 1tiltr ls /t stator copper loss in watts %/ 
i2rr = slipsdhp; /‘t rotor copper loss in watts %/ 


paeff = dhp/(dhp + ph + pe + i2r + izrr); 
printf("\n Endurance efficiency is X%1f", paeff); 


print#("\nSam@e aachine? "); 
scanf("Ad*, &f); 


if (f == 0) 
continue; 
else if (f == 2) 
{ 
e= i; 
break; 
} 
} /t end of f-loop t/ 
} /t end of e-loop %/ 
} /$ end of aain prograa $/ 


double find siv(l, r, g, ds, dcore, 1s) 
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double 1, r, g. ds, dcare, 1s; 
{ 
double one, two, three, four: 


one = (rtgtdstdcore) t{r+gtdstdcore) - (r+q) t(r+q); 
two = 29PIR(rt+q)tdstls; 
three = (r+gtdstdcore)t(r+gtdstdcore) - (r+gtds)t(r+gtds); 


four = 18(PI¥one - two) + PI84s(r+g)tthree; 


return (four); 
} 
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Appendix E. Weight and volume algorithms 


The ability to correctly characterize a new naval ship 
technology depends in part on the ability to calculate the 
weight and volume associated with that technology. Algo- 
rithms for this purpose were taken from a variety of 
sources, including the Advanced Surface Ship Evaluation Tool 
(ASSET) theory manuals. The ASSET algorithms are the result 
of data analysis for naval ships that have been constructed, 
as well as studies for other ship designs. 

The Ship Work Breakdown Structure (SWBS) categorizes 


all ship weights. The general categories are: 


W100 Ship structures 

W200 Propulsion plant 

W300 Electric generation plant 

W400 Command and control equipment 

W500 Auxiliaries and distributed 
systems 

W600 Outfit and furnishings 

W700 Ship armament 

WFOO Ship Loads 


Within W200 are several sub-groups that pertain to 


electric propulsion. They are: 

W235 Electric propulsion devices 

W235.1 Propulsion motors 

W235.2 Propulsion generators 

W235.,.3 Transmission lines and 
propulsion cables 

W235.4 Cooling systems 

W235.5 Switchgear 

W241.1 Locked-train-double-reduction 
reduction gears 

W242 Propulsion clutches and 
couplings 

W243 Shafting 

W244 Propulsion shaft bearings 

W245 Propeller weight 


195 





ASSET allows SWBS groups to be adjusted in weight, 
which allows technology sensitivity analyses, such as this 
thesis, to be performed. Only a few of the W235 sub-groups 
need to be calculated outside of ASSET and adjusted within 
ASSET. These are W235.1, W235.2, W235.4, W235.5, W241.1, 
and W243. The W235.1 and W235.2 weights and volumes were 
calculated as part of the machine design. The rest of the 
needed W235 weights were calculated in "wt.c", a copy of 
which follows. These algorithms are the U. S. Navy 
standard, and were verified against actual ships and com- 
ponents. 

Shafting and transmission line weights are dependent on 
motor and generator positions within a ship. Their weights 
and volumes were calculated from ASSET equations, using the 
layouts of the baseline and variant ships. 

Where no volume equation was found in the ASSET 
documentation, or where the result of such an equation was 
unrealistic, a literature search generally found enough ac- 
tual equipment to permit a relationship to be empirically 
determined. A linear scaling of such volumes provided 


adequate results. 
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Table 44. Listing of off-line weight and volume program 


finclude “stdio.h" 
tdetine K 190.79 /t gear hardness factor %/ 
/¥ prograa nage: wt.c, to find the weights of SWES groups %/ 


main() 

{ 

FILE fp, tfopen(); 

double xq, x@, 2g, z@, pa, nag, na, gra, grag, np, dhp, gn, ds, qg, ne, 
W255, #2595, W254, W2555, w241=0.0, W243, w2979=0.0, xprop, 
wpC, WpC, ppc, wexc, vexc, work, vbrk, pa, v239, 
¥2399=0.0, v239420.0, v2355, v241, v243; 

extern double sqrt(), pow(); 


printf(*\nThis prograa calculates SWBS 200 weights for screen output.”); 
printf("\n\nWhat is the LCG of the propulsion generator (s)? "); 
Beant( 41f°, &%Q): 

printf (*What is the LCb of the propulsion motor(s)? "); 
scant("4lf", &xa); 

printf ("What is the VC6 of the propulsion generator(s)? "); 
scant ("alf", 429); 

printt(*wWhat 1s the VC6 of the propulsion aotor(s)? "); 
scanf(*’lt*, 42a); 

printf ("What is the nuaber of propulsion generator(s)? °); 
scanf(*’lf", 4ng); 

printt#("What is the nugber of propulsion aotor(s)? *); 
scanf("Z1¢", dna); 

printf ("What is the number of gas turbines aboard? "); 
scant("dZit", &ne): 

printt(*What is the rated horsepower of each gas turbine? "); 
scant(*’lf*, dpa); 

printt#("What is the propeller rpa? "): 

scant("%lt°, &np); 

printf (*How auch horsepower is delivered to each propeller? °); 
scant("41¢", &dhp): 

printt(*What is the LC6 of the propeller? *); 

scanf(*Alf", &xprop); 

printt#("What is the gear ratio at the propulsion aotor? "); 
scant("4lf", &gra); 

printf("What is the gear ratio at the propulsion generator? "); 
scant("’lf”, &qrq); 


/t W235.3 Transmission lines %/ 
W2503 = 0,000009% (xa-xgtzg-zat27) ¥(pat746/30000); /¥ LT, enhanced ¢/ 


V2393 = 0.065458 (xa-xgtzg-zat27); /t cubic feet, enhanced %/ 
/t W233.4 Cooling systeas %/ 

w2504 = 0, 26%patng/22405 /t LT t/ 
v2304 = 100. 9%ng; /t swag %/ 


doe 





/¢ W235.9 Switch gear %/ 
W2305 = 0,268(ng + 2hnm); /& enhanced ASSET, LT %/ 


v2899 = 49.08w2355: /¥ switchgear volume, ft*3, enhanced %/ 

ppc = 9.0009395238pa/na; /t power converter rating, MW, enhanced ASSET %/ 

wpc = nats. o5sppc/12; /¥ weight pwr conv, LT, enhanced ASSET %/ 

vpc = naso4sOtppc/i2; /t vol per conv, LT, enhanced ASSET %/ 

pa = pat/46/1000000; /t notor rating in AW %/ 

wexc = (ne + nabSpow((pe/30), 0.3): /t weight of exciters, LT, enhanced %/ 
vexc = /)kwexc; /t vol of exciters, ft*3 %/ 

work = 0. 26%patna: /t weight of braking resistors, LT, enhanced %/ 
vbrk = 37%patna: /t val af resistors, ft*3, enhanced %/ 

WS) = w2355 + w2354 + W255 + pC + Wexc + wbrks 

v259 = v2353 + v2354 + v2359 + vpc + vexc + vbrk: 


/t W24i locked-train double reduction gears %/ 
if (grg '= 1.0) /t there are pg gears, btw gt and pg %/ 
w241 = (L.o78patpow(lgrgti), 3.0)8ng)/ (34008 gragtgrgtargtk); 
if (gra '= 1.0) /t there are pa gears, btw p@ and propeller t/ 
w24i += (L.S78dhptpow((gratl), 3.0)8na)/(nptgratgratk): /t LT %/ 
v24i = w241934.412; /t cubic feet, ratioed fa FFG7 %/ 


/t W243 Shafting %/ 

ds = 2.1528pom( (4.22347 8dhp/np), 0.333); /8 shaft diam, sq inches %/ 

W243 = L.S7088dstdst2t(xprop ~ xa ~ 1.0)/2240; AQ LT t/ 
v243 = 0.01091 Sdstdst(xprop - xa - 6.0); /t ft*3 O/ 


/t 0298.1 LIDR operating fluids, additon to wt in ASSET %/ 
if (grg '= 1.0) 

W298 = 0.27%patng/2240; 
if (gre ‘= 1.0) 


W298 += O.27%patna/2240; /t LT t/ 
Printf("\n\nW239.3 Transmission lines  47.21¢ LT 27.21 ft*3", #2393, v2303); 
print# ("\nW230.4 Cooling systeas Rieeltoep eeraclt ft°3", wesdsy Vest): 
printf ("\nW¥235.9 Switchgear Miecitrlleeavecit ft oS 5 Wedd, V2s00)$ 
printf (*\n Power converters eectsaliehiecit tt 5°, WDE, vpc)s 
printf ("\n Exciters MacltlN ch ecit ft S", WeXC, VEXC)$ 
printf("\n Braking resistors A7.21¢ LT 27.21¢ $t*3", work, vbrk): 
print#(*\n\nW235 Electric propulsion Z7.21¢ LT 47.21f ft*3, less P6s and PMs", w239, v230); 
printf(*\nW241 Reduction gears Rett wonvezit £073", w24l, v24l)s 
print#("\nW243 Shafting biecie Ll ety, 214 ft73", W283, v243)s 
printf("\nW298 Sear operating fluid 47.2) LT\n\n", w298); 

} 
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Appendix F. Advanced Surface Ship Evaluation Tool output 


The output of ASSET is in text and graphic form. The 
total text output for any particular synthesis run is more 
than thirty pages. Following are several graphic outputs of 
ASSET, showing the mechanical and electrical transmission 


ships used in this thesis. 
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Hull isometric view of all ships 


Figure 49. 
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Figure 50. Body plan of all ships 
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Figure 51. 
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